funa191.m

来自「matlab算法集 matlab算法集」· M 代码 · 共 11 行

      function dx = funa191 (t,x)
      %--------------------------------------------------------------- 
      % Van der Pol oscillator 
      %--------------------------------------------------------------- 
         mu = 2;
         dx = zeros(2,1);      % column vector
         dx(1) =  x(2);
         dx(2) = -x(1) + mu*(1 - x(1)*x(1))*x(2);