微分方程数值解第一次报告——matlab改进欧拉法与四阶龙格-库塔求解.docxVIP

微分方程数值解第一次报告徐松松1304一：实验目的掌握MATLAB语言、C/C++语言编写计算程序的方法、掌握改进欧拉法与四阶龙格-库塔求解一阶常微分方程的初值问题。掌握使用MATLAB程序求解常微分方程问题的方法。二：实验内容分别写出改进欧拉法与四阶龙格-库塔求解的算法，编写程序上机调试出结果，要求所编程序适用于任何一阶常微分方程的数值解问题，即能解决这一类问题，而不是某一个问题。? 实验中以下列数据验证程序的正确性。 (0=x=5) , 步长h=0.25。三：源程序改进后欧拉格法程序源代码：function [] = GJOL(h,x0,y0,X,Y)formatlongh=input(h=); x0=input(x0=); y0=input(y0=); disp(输入的范围是); X=input(X=);Y=input(Y=); n=round((Y-X)/h); i=1;x1=0;yp=0;yc=0; fori=1:1:n x1=x0+h; yp=y0+h*(-x0*(y0)^2);%yp=y0+h*(y0-2*x0/y0);% yc=y0+h*(-x1*(yp)^2);%yc=y0+h*(yp-2*x1/yp);% y1=(yp+yc)/2; x0=x1;y0=y1; y=2/(1+x0^2);%y=sqrt(1+2*x0);%