qcd.ode 357 B

12345678910111213141516171819202122
  1. # You may run this example by doing
  2. #
  3. # ode < qcd.ode | graph -T X -C
  4. #
  5. # This simulates a coupled system of ordinary differential equations that
  6. # arises in the mathematical theory of quantum chromodynamics (QCD).
  7. #
  8. # Source: Griffiths et al., J. Math. Phys.
  9. f' = fp
  10. fp' = -f*g*g
  11. g' = gp
  12. gp' = g*f*f
  13. f = 0
  14. fp = -1
  15. g = 1
  16. gp = -1
  17. print t, f
  18. step 0, 5