[lxcwym]

function dy=plas(x,y)
dy=zeros(5,1);
y(1)=1/(2*y(3)^2);
y(2)=log(y(3));
dy(1)=y(4);
dy(2)=y(5);
dy(4)=y(3)-1-dy(5);

>>[x,y]=ode45(@plas,[-inf,inf],[1,0])

我这样编写，问题出现在哪里呢？谢谢！