lorenz.m

来自「lorenz吸引子的画法」· M 代码 · 共 5 行

function dy = lorenz(t,y)
dy = zeros(3,1); 
dy(1) = 10*(-y(1)+y(2));
dy(2) = 28*y(1)-y(2)-y(1)*y(3);
dy(3) = y(1)*y(2)-8*y(3)/3;