home *** CD-ROM | disk | FTP | other *** search
/ ProfitPress Mega CDROM2 …eeware (MSDOS)(1992)(Eng) / ProfitPress-MegaCDROM2.B6I / MISC / EDUCATIO / ODE25.ZIP / GEISEL.ODE < prev    next >
Encoding:
Text File  |  1990-01-21  |  430 b   |  21 lines
  1. # GEISEL.ODE
  2. # T. Geisel, Phys Rev Lett 59, 2503 (1988). A diffusive chaotic system.
  3. # An exquisite example!!
  4. # Plot x vs. y and keep continuing the solution with F7.
  5. # A particle moves in a 2D periodic potential.
  6. # Author claims spectrum is 1/f.
  7.  
  8. $m 5
  9. $n 500
  10. $t 200
  11.  
  12. # Particle position is (x,y)
  13. x   = 0
  14. y   = 0
  15. vx  = 1
  16. vy  = 0.5
  17. x'  = vx
  18. y'  = vy
  19. vx' = (1.5+0.5*cos(y))*sin(x)
  20. vy' = (1.5+0.5*cos(x))*sin(y)
  21.