我初学 各位大侠帮忙看一下拉 不尽感激

vvyin

<>function xdot = lorenzeq(t,x)
xdot=[-8/3*x(1)+x(2)*x(3);
      -10*x(2)+10*x(3);
      -x(1)*x(2)+28*x(2)-x(3)];</P>
<>t_final=100; x0=[0;0;1e-10];
   [t,x]=ode45('lorenzeq',[0,t_final],x0);
   plot(t,x),
   figure; plot3(x(:,1),x(:,2),x(:,3)); axis([10 40 -20 20 -20 20]);</P>
<>这个程序错在那里的呀 我是初学 实在不知道了 各位帮忙看下好不？？
</P>[em01]

vvyin

<>题目是要画出那个方程组的解图</P>