lorenc.m

Show work of ODE45 function. As the example use lorenz equation. Example show ensemble Puankare and

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);