%%% Nonlinear Systems Analysis, NTU-EE, 91 Fall
%%% Dynamic Equations of Van der Pol oscillator
%%% Date: Sep. 26, 2002
%%% By Feng-Li Lian, fengli@ntu.edu

function dx = funVanderPol(t,x)

dx = zeros(2,1);   

epsilon = 1.0;

dx(1) = x(2);
dx(2) = - x(1) + epsilon * ( 1 - x(1)^2 ) * x(2);