File: invpend.ode

package info (click to toggle)
xppaut 8.0%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 7,332 kB
  • sloc: ansic: 74,690; makefile: 127; sh: 92
file content (12 lines) | stat: -rwxr-xr-x 265 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
 
# the inverted pendulum driven periodically
# start at 0 and give it a push
par a=.3,omega=10,k=.1
th'=v
v'=-k*v-(1+a*omega*omega*cos(omega*t))*sin(th)
init th=0,v=3
aux mth=mod(th,2*pi)
@ meth=qualrk,total=30
@ xp=t,yp=mth,xlo=0,xhi=30,ylo=0,yhi=6.3
@ lt=-2
done