function xdot=lorenz1(t,x,flag,beta,rho,sigma)
xdot=[-beta*x(1)+x(2)*x(3);
-rho*x(2)+rho*x(3);
-x(1)*x(2)+sigma*x(2)-x(3)];
t_final=100;x0=[0;0;1e-10];
b1=8/3;r1=10;s1=28;
[t,x]=ode45('lorenz1',[0,t_final],x0,[],b1,r1,s1);
plot(t,x)
figure;plot3(x(:,1),x(:,2),x(:,3));
axis([10 42 -20 20 -20 25]);
谢谢您们指出错在哪里