实验五常微分方程值数解.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为输出值。 计算流程框图 程序清单 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 = 0.9820 x1 = 0.0400 y1 = 0.9650 x1 = 0.0600 y1 =