dydt2.m

来自「一个小小的c语言编程」· M 代码 · 共 10 行

function DX=DYDt2(t,X,flag,sigma,b,r)
switch flag
    case ''
        DX=[-sigma*X(1)+sigma*X(2);r*X(1)-X(2)-X(1).*X(3);X(1).*X(2)-b*X(3)];
    case 'init'
        options.RelTol=1e-5;options.AbsTol=1e-5;
    otherwise 
        error(['Unknown flag ''' flag '''.']);
        
end