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 '$gauss 1 '$version)
16
17(putprop '$gauss_a
18	 `((a b c x)
19	   ((unk) "$first" "$gauss_a")
20	   ((unk) "$second" "$gauss_a")
21	   ((unk) "$third" "$gauss_a")
22	   (($dgauss_a) a b c x))
23	 'grad)
24
25(putprop '$dgauss_a
26	 `((a b c x)
27	   ((unk) "$first" "$dgauss_a")
28	   ((unk) "$second" "$dgauss_a")
29	   ((unk) "$third" "$dgauss_a")
30	   ((mtimes) ((mexpt) ((mplus) 1 ((mtimes) -1 x)) -1)
31	    ((mexpt) x -1)
32	    ((mplus) ((mtimes) a (($gauss_a) a b c x) b)
33	     ((mtimes) (($dgauss_a) a b c x)
34	      ((mplus) ((mtimes) -1 c)
35	       ((mtimes) ((mplus) 1 a b) x))))))
36	 'grad)
37
38(putprop '$gauss_b
39	 `((a b c x)
40	   ((unk) "$first" "$gauss_b")
41	   ((unk) "$second" "$gauss_b")
42	   ((unk) "$third" "$gauss_b")
43	   (($dgauss_b) a b c x))
44	 'grad)
45
46(putprop '$dgauss_b
47	 `((a b c x)
48	   ((unk) "$first" "$dgauss_b")
49	   ((unk) "$second" "$dgauss_b")
50	   ((unk) "$third" "$dgauss_b")
51	   ((mtimes) ((mexpt) ((mplus) 1 ((mtimes) -1 x)) -1)
52	    ((mexpt) x -1)
53	    ((mplus) ((mtimes) a (($gauss_b) a b c x) b)
54	     ((mtimes) (($dgauss_b) a b c x)
55	      ((mplus) ((mtimes) -1 c)
56	       ((mtimes) ((mplus) 1 a b) x))))))
57	 'grad)
58