File: geisel.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 (10 lines) | stat: -rw-r--r-- 219 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
 
# an example due to Geisel with 1/f spectral properties
par a=1.5,b=.5
x'=u
y'=v
u'=(a+b*cos(y))*sin(x)
v'=(a+b*cos(x))*sin(y)
init x=-30,y=-50,u=1,v=.5
@ total=8000,dt=.25,maxstor=100000,bound=100000,meth=qualrk
done