matlab解二阶微分方程组ode,MATLAB中用ODE45解两个2阶非线性常微分方程组,Command Window出现警告：写在补充里面>> [t,y]=ode45(@RK4_shiyan01,...

MATLAB中用ODE45解两个2阶非线性常微分方程组,Command Window出现警告：写在补充里面>> [t,y]=ode45(@RK4_shiyan01,[0:0.001:1000],[0,0,0,0])Warning: Failure at t=1.768731e-001. Unable to meet integration tolerances without reducin

2016-01-27

MATLAB中用ODE45解两个2阶非线性常微分方程组,Command Window出现警告：写在补充里面

>> [t,y]=ode45(@RK4_shiyan01,[0:0.001:1000],[0,0,0,0])

Warning: Failure at t=1.768731e-001. Unable to meet integration tolerances without reducing the step

size below the smallest value allowed (4.440892e-016) at time t.

> In ode45 at 371

方程如下：

function dx=RK4_shiyan01(t,x)

a1=4.73001;a2=7.85320;b1=0.005;b2=0.2;b3=10;b4=50;b5=7.6;b6=20;b7=0.05;b8=20;b12=0.8355;b21