# cubic nullcline
param tau=1
g(x) = x^3/3-x
du/dt = v - g(u)
dv/dt = -u/tau
init u=1 v=0
@ xp=u yp=v xlo=-3 xhi=3 ylo=-3 yhi=3
done