Work at SourceForge, help us to make it a better place! We have an immediate need for a Support Technician in our San Francisco or Denver office.

## maxima-bugs

 [Maxima-bugs] [ maxima-Bugs-736540 ] gradef for bessel functions From: SourceForge.net - 2003-05-12 16:32:00 Bugs item #736540, was opened at 2003-05-12 11:31 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=736540&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Barton Willis (willisb) Assigned to: Nobody/Anonymous (nobody) Summary: gradef for bessel functions Initial Comment: Consider (C1) display2d : false; Evaluation took 0.00 seconds (0.00 elapsed) (D1) FALSE (C2) diff(bessel_j(x,x),x); Evaluation took 0.00 seconds (0.00 elapsed) (D2) 'DIFF(BESSEL_J[x](x),x,1)-BESSEL_J[x](x) +BESSEL_J[x-1](x) The derivative in the first term of (d2) is should be with respect to the order of the bessel function -- instead it's a _total_ derivative. A good solution isn't easy; in orthopoly, I handle this problem by signaling an error. Thus (from orthopoly) ;; When a user requests the derivative of an a function in this package ;; with respect to the order or some other parameter, return a form ;; ((unk) input from user). We "simplify" this form by printing an error. (defprop unk simp-unk operators) (defun simp-unk (x y z) (declare (ignore y z)) (merror "Maxima doesn't know the derivative of ~:M with respect the ~:M argument" (nth 2 x) (nth 1 x))) (putprop '\$legendre_p '((n x) ((unk) "\$first" "\$legendre_p") ((mtimes) ((mplus) ((mtimes) n ((\$legendre_p) ((mplus) -1 n) x)) ((mtimes) -1 n ((\$legendre_p) n x) x)) ((mexpt) ((mplus) 1 ((mtimes) -1 ((mexpt) x 2))) -1))) 'grad) (C3) build_info(); Maxima version: 5.9.0 Maxima build date: 19:10 2/9/2003 host type: i686-pc-mingw32 lisp-implementation-type: Kyoto Common Lisp lisp-implementation-version: GCL-2-5.0 Evaluation took 0.00 seconds (0.00 elapsed) (D3) (C4) Barton ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=736540&group_id=4933
 [Maxima-bugs] [ maxima-Bugs-736540 ] gradef for bessel functions From: SourceForge.net - 2003-05-25 15:10:43 Bugs item #736540, was opened at 2003-05-12 12:31 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=736540&group_id=4933 Category: None Group: None >Status: Closed >Resolution: Fixed Priority: 5 Submitted By: Barton Willis (willisb) Assigned to: Nobody/Anonymous (nobody) Summary: gradef for bessel functions Initial Comment: Consider (C1) display2d : false; Evaluation took 0.00 seconds (0.00 elapsed) (D1) FALSE (C2) diff(bessel_j(x,x),x); Evaluation took 0.00 seconds (0.00 elapsed) (D2) 'DIFF(BESSEL_J[x](x),x,1)-BESSEL_J[x](x) +BESSEL_J[x-1](x) The derivative in the first term of (d2) is should be with respect to the order of the bessel function -- instead it's a _total_ derivative. A good solution isn't easy; in orthopoly, I handle this problem by signaling an error. Thus (from orthopoly) ;; When a user requests the derivative of an a function in this package ;; with respect to the order or some other parameter, return a form ;; ((unk) input from user). We "simplify" this form by printing an error. (defprop unk simp-unk operators) (defun simp-unk (x y z) (declare (ignore y z)) (merror "Maxima doesn't know the derivative of ~:M with respect the ~:M argument" (nth 2 x) (nth 1 x))) (putprop '\$legendre_p '((n x) ((unk) "\$first" "\$legendre_p") ((mtimes) ((mplus) ((mtimes) n ((\$legendre_p) ((mplus) -1 n) x)) ((mtimes) -1 n ((\$legendre_p) n x) x)) ((mexpt) ((mplus) 1 ((mtimes) -1 ((mexpt) x 2))) -1))) 'grad) (C3) build_info(); Maxima version: 5.9.0 Maxima build date: 19:10 2/9/2003 host type: i686-pc-mingw32 lisp-implementation-type: Kyoto Common Lisp lisp-implementation-version: GCL-2-5.0 Evaluation took 0.00 seconds (0.00 elapsed) (D3) (C4) Barton ---------------------------------------------------------------------- >Comment By: Raymond Toy (rtoy) Date: 2003-05-25 11:10 Message: Logged In: YES user_id=28849 Derivative with respect to order added for Bessel J. It's a bit messy. ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=736540&group_id=4933