home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 9 / 09.iso / e / e004 / 1.img / ANHARM.MCD next >
Encoding:
Text File  |  1987-10-21  |  886 b   |  39 lines
  1. .MCD 20000 0
  2. .CMD PLOTFORMAT logs=0,0 subdivs=1,1 size=5,15 type=l
  3. .CMD FORMAT  rd=d ct=10 im=i et=3 zt=15 pr=3 mass length time charge
  4. .CMD SET ORIGIN 0
  5. .CMD SET TOL 0.001000
  6. .CMD MARGIN 0
  7. .CMD LINELENGTH 78
  8. .CMD SET PRNCOLWIDTH 8
  9. .CMD SET PRNPRECISION 4
  10. .TXT 0 0 1 27 
  11. a1,26,24,25
  12. Anharmonic Oscillations:
  13. .TXT 2 1 1 28 
  14. a1,28,26,27
  15. Over a mesh of N points...
  16. .EQN 0 33 1 10 
  17. N:100
  18. .EQN 0 14 1 12 
  19. i:1;N
  20. .TXT 2 -40 1 43 
  21. a1,43,41,42
  22. ...spanning the range from -α to α, where
  23. .EQN 0 43 1 10 
  24. α:3.1
  25. .EQN 2 -15 3 23 
  26. x[i:α*(2*i/(N+1)-1)
  27. .EQN 1 -32 2 9 
  28. t[0:0
  29. .EQN 0 16 2 9 
  30. x[0:α
  31. .TXT 3 -19 2 77 
  32. a2,77,75,134
  33. We compute the times at which these points are reached by an oscillator 
  34. having a cosine potential energy, and plot their trajectory.
  35. .EQN 3 1 5 34 
  36. t[i:t[(i-1)+α/\|(cos(x[i)-cos(α))
  37. .EQN 0 40 11 35 
  38. α&-α&x[i{1,1,9,30,d}@t[N&0&t[i
  39.