function ydot=vdpol(t,y)
global MUTEMP
mu=MUTEMP;
ydot = [y(2); mu*(1-y(1)^2)*y(2)-y(1)];