lorenz.m

来自「初学者可以参考一下 作为MATLAB的入门教程与实践」· M 代码 · 共 11 行

t_final=100; 
x0=[0; 0; 1e-10];
[t,x] = ode45('lorenzeq',[0, t_final], x0);
figure(1); set(gcf,'position',[7 225 390 270])
plot(t,x),
figure(2); set(gcf,'position',[405 225 390 270])
plot3(x(:,1), x(:,2), x(:,3)); 
axis([10 40 -20 20 -20 20])