function dydt = vdpeqMU(t,y,mu) % Evaluate the van der Pol ODEs for mu = 0.2 global mu dydt = [y(2); mu*(1-y(1)^2)*y(2)-y(1)]; % "Complex and Chaotic Nonlinear Dynamics. % Advances in Economics and Finance, % Mathematics and Statistics" % T.Vialar, Springer 2009 % Copyright(c).