常微分方程的常见解法演示文稿.pptVIP

令 令 ) ( 2 1 1 1 - - n n x y K x 预测 处的斜率 如果以 本文档共97页；当前第63页；编辑于星期一\16点20分 取 则 -----------(6) (6)式称为三阶Runge-Kutta方法 本文档共97页；当前第64页；编辑于星期一\16点20分 还可构造四阶(经典)Runge-Kutta方法 四阶(经典)Runge=Kutta方法有4阶精度 本文档共97页；当前第65页；编辑于星期一\16点20分 例 求初始值问题的数值解 利用四阶Runge=Kutta方法计算机编程 给出步长和初始值 循环计算各点上函数的近似值 显示结果 本文档共97页；当前第66页；编辑于星期一\16点20分 printlev1:=0: h:=0.1: x[0]:=0: y[0]:=0.5: f:=(x,y)-1+(y-x)^2; for n from 1 to 10 do x[n]:=h*n; k1:=f(x[n-1],y[n-1]); k2:=f(x[n-1]+h/2,y[n-1]+k1*h/2); k3:=f(x[n-1]+h/2,y[n-1]+k2*h/2); k4:=f(x[n-1]+h,y[n-1]+k3*h); y[n]:=y[n-1]+h*(k1+2*k2+2*k3+k4)/6; u[n]:=x[n]+