# phase_app.ode
# phase model for two coupled oscillators
# using numerically computed H
h(x)=3.34-3.05*cos(x)-.29*cos(2*x)+3.61*sin(x)-.33*sin(2*x)
th1'=w1+a*h(th2-th1)
th2'=w2+a*h(th1-th2)
par w1=1,w2=1,a=.1
@ total=100
@ fold=th1,fold=th2
@ xlo=0,ylo=0,xhi=6.3,yhi=6.3,xp=th1,yp=th2
done