x'=v
v'=g-(k*x+(f+mu)*v)/m
m'=f
global 1 x-1 {m=max(m-h*m*v,m0)}
par f=.4,g=.32,h=4,m0=.01,mu=.6,k=1
init x=0,v=0,m=1
@ total=200,bound=1000000
done