lorenc.m

来自「Show work of ODE45 function. As the exam」· M 代码 · 共 12 行

function  dY=lorenc(t,y)

sigma=10;
b=8/3;
r=16;

dY=zeros(3,1);

dY(1)=sigma*(y(2)-y(1));
dY(2)=r*y(1)-y(2)-y(1)*y(3);
dY(3)=-b*y(3)+y(1)*y(2);