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)];
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 42 -20 20 -20 25]);
错在哪里啊？