ff4.m

来自「求解方程 包括求解线性方程非线性方程」· M 代码 · 共 3 行

    [T,Y]=ode15s('vdp1000',[0 3000],[2 0]); 
     plot(T,Y(:,1),'-')