chaoslorenz.m

来自「研究生期间的数值计算的书配套程序源代码 奉献了」· M 代码 · 共 12 行

clear
clc
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]);
figure;
comet3(x(:,1),x(:,2),x(:,3)); 
axis([10 40 -20 20 -20 20]);