龙格库塔方法的应用.doc

龙格库塔方法的应用
例如求解微分方程y'=y-(2*x)/y在区间[0 1]上的数值解，初始条件y(1)=1
在matlab中的运行结果如下
h =

0.2000
x =

0
y =

1

微分方程y’=y-(2*x)/y在点
0

处的数值解为
1

微分方程y’=y-(2*x)/y在点
0.2000

处的数值解为
1.1832

微分方程y’=y-(2*x)/y在点
0.4000

处的数值解为
1.3417

微分方程y’=y-(2*x)/y在点
0.6000

处的数值解为
1.4833

微分方程y’=y-(2*x)/y在点
0.8000

处的数值解为
1.6125

微分方程y’=y-(2*x)/y在点
1

处的数值解为
1.7321