龙格库塔方法的应用
例如求解微分方程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