# vdpprc.ode
# PRC of the van der Pol oscillator
init x=2
x'=y
y'=-x+y*(1-x^2)+a*pulse(t-tau)
tau'=0
pulse(t)=heav(t)*heav(sigma-t)
par sigma=.2,a=0
par t0=6.65
aux prc=1-t/t0
@ dt=.01
done