lorenz.m

来自「1）先利用文件Lorenz.m求解出Lorenz系统」· M 代码 · 共 7 行

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