[wumaobo]

目标函数
function dy=rigid(t,y)
global g A as vp vs ap;
dy=zeros(2,1);
dy(1)=g*y(1)*y(2)/A-as*y(1);
dy(2)=-vp*g*y(1)*y(2)/(vs*A)-ap*y(2);
主函数
global g A as vp vs ap ;
g=0.609*10^(-13);A=8*10^(-13);as=-25;ap=4*10^8;vp=1550;vs=1450;
[T,Y]=ode45(@rigid,[0 40],-25);
plot(T,Y)
请问那里出错