% myvdpol.m, m-function by ralph, 14 oct 2003
function dydt = myvdpol(t,y,Mu);
% exemplary dyn scheme (2D flow with one param)
dydt = [y(2); (Mu*1-y(1)^2)*y(2)-y(1)];
end;