[Maxima-commits] CVS: maxima/src bessel.lisp, 1.62, 1.63 expintegral.lisp, 1.9, 1.10

[David B. <bil...@us...>]

Update of /cvsroot/maxima/maxima/src
In directory 23jxhf1.ch3.sourceforge.com:/tmp/cvs-serv21737

Modified Files:
	bessel.lisp expintegral.lisp 
Log Message:
Move definition of $expint from bessel.lisp to expintegral.lisp.


Index: bessel.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/src/bessel.lisp,v
retrieving revision 1.62
retrieving revision 1.63
diff -u -d -r1.62 -r1.63
--- bessel.lisp	12 Dec 2008 11:28:23 -0000	1.62
+++ bessel.lisp	12 Dec 2008 11:33:22 -0000	1.63
@@ -672,18 +672,6 @@
     (print-invert-case (implode (mstring `(($random_normal) ,@ (cdr form))))))
   '$done)
 
-;; I think this is the function E1(x).  At least some simple numerical
-;; tests show that this expint matches the function de1 from SLATEC
-
-;; Exponential integral E1(x).  The Cauchy principal value is used for
-;; negative x.
-(defun $expint (x)
-  (cond ((numberp x)
-	 (values (slatec:de1 (float x))))
-	(t
-	 (list '($expint simp) x))))
-
-
 ;; Define the Bessel funtion J[n](z)
 
 (defprop %bessel_j simp-bessel-j operators)

Index: expintegral.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/src/expintegral.lisp,v
retrieving revision 1.9
retrieving revision 1.10
diff -u -d -r1.9 -r1.10
--- expintegral.lisp	5 Dec 2008 23:15:09 -0000	1.9
+++ expintegral.lisp	12 Dec 2008 11:33:22 -0000	1.10
@@ -16,6 +16,8 @@
 ;;;  $expintegral_shi (z) - Exponential Integral Shi(z)
 ;;;  $expintegral_chi (z) - Exponential Integral Chi(z)
 ;;;
+;;;  $expint (x)          - Exponential Integral E1(x) (depreciated)
+;;;
 ;;;  Global variables for the Maxima User:
 ;;;
 ;;;  $expintrep    - Change the representation of the Exponential Integral to
@@ -1793,3 +1795,13 @@
         (mul -1 ($log z))))))
 
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;; Moved from bessel.lisp 2008-12-11.  Consider deleting it.
+;;
+;; Exponential integral E1(x).  The Cauchy principal value is used for
+;; negative x.
+(defun $expint (x)
+  (cond ((numberp x)
+	 (values (slatec:de1 (float x))))
+	(t
+	 (list '($expint simp) x))))
+