#dae_ex5.ode
w'=v_
0= v_*(1-v_*v_)-w
solv v_=1
aux v=v_
@ NEWT_ITER=1000,NEWT_TOL=1e-3,JAC_EPS=1e-5,METH=qualrk
@ total=10,dt=.01
@ xlo=-1.5,xhi=1.5,ylo=-.5,yhi=.5,xp=v,yp=w
done