f1=inline(['[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(1)*x(2)+28*x(2)_x(3)]'],'t','x');
t_final=100;x0=[0;0;1e-10];[t,x]=ode45(f1,[0,t_final],x0);
plot(t,x),figure;
plot3(x(:,1),x(:,2),x(:,3));
也不知错在哪里啊