# cannon.ode
# with linear friction
x'=vx
vx'=-f*vx
y'=vy
vy'=-f*vy-g
global 0 t  {vx=v0*cos(pi*theta/180);vy=v0*sin(pi*theta/180)}
par theta=30
par f=.01,g=1
par v0=2.5
init x=0,y=0
@ xlo=0,xhi=5
@ ylo=0,yhi=3
@ xp=x,yp=y,bound=100000
done