ex1_3.m

源代码集合

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 40 -20 20 -20 20]);