clearallClcr=设定值；f=@(t,x)[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(1)*x(2)+r*x(2)-x(3)];t_final=100;x0=[0;0;1e-10][t,x]=ode45(f,[0,t_final],x0)figure(2)plot3(x(:,1),x(:,2),x(:,3));axis([1042-2020-2020])

　　clearall

　　Clc

　　r=设定值；

　　f=@(t,x)[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(1)*x(2)+r*x(2)-x(3)];

　　t_final=100;

　　x0=[0;0;1e-10]

　　[t,x]=ode45(f,[0,t_final],x0)

　　figure(2)

　　plot3(x(:,1),x(:,2),x(:,3));

　　axis([1042-2020-2020])