实验五常微分程数值解.docVIP

实验五 常微分方程数值解 欧拉法 算法说明 对于xi，i=0,1,2，…，n，取步长h为定值时，有h=xi+1-xi，EURLER法的计算公式为：yi+1=yi+h*f(xi,yi)；i=0,1,2，…，n。 程序中主要符号说明 a为x的下界， b为x的上界， h为步长， n为循环次数（即为x的数值点数减一）， x0、y0为循环初始值， x1、y1为输出值。 计算流程框图 开始 开始 读入数据x0，y0，h，n FOR I=1 TO N X0+H=X1,Y0+H*F(X0,Y0)=Y1 输出X1，Y1 X1=X0 Y1=Y0 NEXT 结束 程序清单 syms a b h n x0 y0 x1 y1; a=0; b=0.1; n=5; x0=0; y0=1; h=(b-a)/n; for j=1:n x1=x0+h; y1=y0-h*0.9*y0/(1+2*x0); x1 y1 x0=x1; y0=y1; end 算例及输入数据说明 算例：求解初值问题 当x=0,0.02,0.04,…,0.10时的数值解。 输入数据说明： a=0; a为x 的下界 b=0.1; b为x的上界 n=5; n为循环次数，即为x的数值点数减一 x0=0; x（0）的值 y0=1; y(0)的值 程序运行结果及结果分析 运行结果： x1 = 0.0200 y1