5楼: Originally posted by laosam280 at 2014-02-27 11:38:50
你这个方程右端因为有一个很小的系数，使得方程的解在很长的时间内都稳定在初值附近。你的初值给的是0.1918，要想解减少到哪怕0.09，x会上好几十万这么大。我刚刚试了一下。
代码都给你了，你自己去matlab上运行就 ...

2楼: Originally posted by laosam280 at 2014-02-27 11:16:58
function SolveRK4
clear all; clc;
y0 = 0.11;
T = 50;
h = 1/20;
=ode45(@myfun,,y0);

plot(t,y,'linewidth',3);
xlabel('x');
ylabel('y(x)')

function dt = myfun(t,y)
dt = ;

untitled.jpg
...

2楼: Originally posted by laosam280 at 2014-02-27 11:16:58
function SolveRK4
clear all; clc;
y0 = 0.11;
T = 50;
h = 1/20;
=ode45(@myfun,,y0);

plot(t,y,'linewidth',3);
xlabel('x');
ylabel('y(x)')

function dt = myfun(t,y)
dt = ;

untitled.jpg
...