lorenz.m

1）先利用文件Lorenz.m求解出Lorenz系统

function fty=lorenz(t,u)
det=10;
b=8/3;
r=28;
fty=[det*(u(2)-u(1))
    r*u(1)-u(2)-u(1)*u(3)
    u(1)*u(2)-b*u(3)];