File: invpend.ode

package info (click to toggle)
xppaut 5.85-3.3
  • links: PTS
  • area: main
  • in suites: etch, etch-m68k
  • size: 11,568 kB
  • ctags: 7,729
  • sloc: ansic: 80,804; cpp: 965; makefile: 271
file content (9 lines) | stat: -rw-r--r-- 195 bytes parent folder | download | duplicates (5) 
1
2
3
4
5
6
7
8
9
 
# 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
@ meth=qualrk,total=30
done