share/contrib/diffequations/kummer.lisp

   1 ;; Author Barton Willis
   2 ;; University of Nebraska at Kearney
   3 ;; Copyright (C) 2004, Barton Willis
   4
   5 ;; Brief Description: Maxima code for linear homogeneous second order
   6 ;; differential equations.
   7
   8 ;; Maxima odelin is free software; you can redistribute it and/or
   9 ;; modify it under the terms of the GNU General Public License,
  10 ;; http://www.gnu.org/copyleft/gpl.html.
  11
  12 ;; Maxima odelin has NO WARRANTY, not even the implied warranty of
  13 ;; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
  14
  15 ($put '$kummer 1 '$version)
  16
  17 (putprop '$kummer_m
  18          '((a b x)
  19            ((unk) "$order" "$kummer_m")
  20            ((unk) "$order" "$kummer_m")
  21            (($dkummer_m) a b x))
  22          'grad)
  23
  24 (putprop '$dkummer_m
  25          '((a b x)
  26            ((unk) "$order" "$kummer_m")
  27            ((unk) "$order" "$kummer_m")
  28            ((mtimes ) ((mexpt  ) x -1)
  29             ((mplus ) ((mtimes  ) a
  30                        (($kummer_m ) a b x))
  31              ((mtimes  ) -1 (($dkummer_m ) a b x) b)
  32              ((mtimes  ) (($dkummer_m ) a b x) x))))
  33          'grad)
  34
  35 (putprop '$kummer_u
  36          '((a b x)
  37            ((unk) "$order" "$kummer_u")
  38            ((unk) "$order" "$kummer_u")
  39            (($dkummer_u) a b x))
  40          'grad)
  41
  42 (putprop '$dkummer_u
  43          '((a b x)
  44            ((unk) "$order" "$dkummer_u")
  45            ((unk) "$order" "$dkummer_u")
  46            ((mtimes ) ((mexpt  ) x -1)
  47             ((mplus ) ((mtimes  ) a
  48                            (($kummer_u ) a b x))
  49              ((mtimes  ) -1 (($dkummer_u ) a b x) b)
  50              ((mtimes  ) (($dkummer_u ) a b x) x))))
  51          'grad)
  52
  53
  54