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])