ode.scm 1.7 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152
  1. #| -*-Scheme-*-
  2. Copyright (C) 1986, 1987, 1988, 1989, 1990, 1991, 1992, 1993, 1994,
  3. 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005,
  4. 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2013 Massachusetts
  5. Institute of Technology
  6. This file is part of MIT/GNU Scheme.
  7. MIT/GNU Scheme is free software; you can redistribute it and/or modify
  8. it under the terms of the GNU General Public License as published by
  9. the Free Software Foundation; either version 2 of the License, or (at
  10. your option) any later version.
  11. MIT/GNU Scheme is distributed in the hope that it will be useful, but
  12. WITHOUT ANY WARRANTY; without even the implied warranty of
  13. MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  14. General Public License for more details.
  15. You should have received a copy of the GNU General Public License
  16. along with MIT/GNU Scheme; if not, write to the Free Software
  17. Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301,
  18. USA.
  19. |#
  20. #|
  21. ;;; Let (sigma t) be the state of a system at time t. Let the
  22. ;;; (first-order) system of differential equations governing the
  23. ;;; evolution of this state be:
  24. ;;; ((D sigma) t) = (R (sigma t))
  25. ;;; or (D sigma) = (compose R sigma)
  26. ;;; i.e. R is a system derivative.
  27. ;;; Let F be any function of state, then a differential equation for
  28. ;;; the evolution of F, as it is dragged along the integral curve
  29. ;;; sigma is:
  30. ;;; (D (compose F sigma)) = (* (compose (D F) sigma) (D sigma))
  31. ;;; = (compose (* (D F) R) sigma)
  32. ;;; Let's call this operation Lie-D (the Lie derivative for
  33. ;;; coordinates). We define
  34. |#
  35. (define (Lie-D R)
  36. (define (the-LD F) (* (D F) R))
  37. (make-operator the-LD `(Lie-D ,R)))