matlab解状态方程组

一、四阶龙格库塔法

u 四阶龙格库塔原理：

u y (i+1)=y( i )+k*h  //h 为步长， K 为平均系数

u a=0;b=20;n=1000;h=(b-a)/ n;T = a:h:b;X =zeros(1,n+1);

u Y=zeros(1,n+1);f=@( t,x,y )(1*y);g=@( t,x,y )(-2*x-3*y+1);X(