The new Bessel code (Revision 1.58) is available on Sourceforge.net. Here is the
link directly to the code:
http://maxima.cvs.sourceforge.net/maxima/maxima/src/bessel.lisp?revision=1.58&vi
ew=markup
There a lot of relationsships between the Bessel functions of negative and
positve arg or order. Because we have only numerical Slatec routines for
positive order and positive or complex arg, these relationsships have been used
to calculate the Bessel functions for negative order and arg. You can found most
of them as comments in the source file. The formulas are taken from A&S.
Dieter Kaiser
-----Urspr?ngliche Nachricht-----
Von: rvh2007 at comcast.net [mailto:rvh2007 at comcast.net]
Gesendet: Freitag, 13. Juni 2008 23:42
An: Dieter Kaiser
Cc: maxima at math.utexas.edu
Betreff: Re: AW: [Maxima] Bessel plotting problem
Could you send me where the new bessel code is by sending me the link so I can
get it? I thought I did this before but I think I got the wrong one.
r
------------Original Message------------
From: "Dieter Kaiser" <drdieterkaiser at web.de>
To: "'Richard Hennessy'" <rvh2007 at comcast.net>
Cc: maxima at math.utexas.edu
Date: Tue, Jun-10-2008 6:02 PM
Subject: AW: [Maxima] Bessel plotting problem
The Bessel functions are defined for all real and complex numbers of the order
and the argument. That's not the problem. But Maxima has limitations. The NEW
CVS Bessel code gives numerical values for all negative and positive orders and
real or complex arguments. The evaluation for complex order is still missing.
The error of the slatec routines below you can generated with the OLD Bessel
code.
Dieter Kaiser
-----Urspr?ngliche Nachricht-----
Von: maxima-bounces at math.utexas.edu [mailto:maxima-bounces at math.utexas.edu] Im
Auftrag von Richard Hennessy
Gesendet: Dienstag, 10. Juni 2008 23:37
An: fateman at cs.berkeley.edu
Cc: Maxima List
Betreff: Re: [Maxima] Bessel plotting problem
I am going to take my conclusion back. Just because Maxima gives an error
message while plotting does not mean the function is only defined for positive
x.
Learn, live.
R
------------Original Message------------
From: "Richard Hennessy"<rvh2007 at comcast.net>
To: fateman at cs.berkeley.edu
Date: Sat, Jun-7-2008 11:03 AM
Subject: RE: [Maxima] Bessel plotting problem
I think I made a mistake. What I wanted to plot was the real valued function
(f+g+conjugate(f+g))/2 where f and g are linearly independent (complex)
solutions to my differential eq. I tried using realpart() but that did not
work. I should have tried the following.
plot2d(-(sqrt(3)/2-%i/2)*%i*'(bessel_y(1/6,-x^3/3))+'(bessel_j(1/6,-x^3/3))*sqrt
(x),[x,-20,.001], plot_format,gnuplot]), orthopoly_returns_intervals : false;
That worked.
The domain of the bessel_j or bessel_y function must be positive x only since
when I try
plot2d(-(sqrt(3)/2-%i/2)*%i*'(bessel_y(1/6,-x^3/3))+'(bessel_j(1/6,-x^3/3))*sqrt
(x),[x,-15,15],[plot_format,gnuplot]), orthopoly_returns_intervals : false;
I get
-->
***MESSAGE FROM ROUTINE DBESY IN LIBRARY SLATEC. ***POTENTIALLY RECOVERABLE
ERROR, PROG ABORTED, TRACEBACK REQUESTED * X LESS THAN OR EQUAL TO ZERO *
ERROR NUMBER = 2 * ***END OF MESSAGE ***JOB ABORT DUE TO UNRECOVERED ERROR.0
ERROR MESSAGE SUMMARY LIBRARY SUBROUTINE MESSAGE START NERR
LEVEL COUNT SLATEC DBESY X LESS THAN OR EQUAL 2 1
1
When x is positive, the plot shows only the second and third quadrants with
actual points. I think the first a fourth are not in the domain since I get
this ugly error message when it is done.
I would call this a bug, I am not sure if it is. Is there a reason for the
restriction x > 0?
Rich
------------Original Message------------
From: "Richard Fateman" <fateman at cs.berkeley.edu>
To: "'Richard Hennessy'" <rvh2007 at comcast.net>
Date: Sat, Jun-7-2008 10:11 AM
Subject: RE: [Maxima] Bessel plotting problem
I haven't tried it, but how about removing the ' and
doing plot2d( realpart( ....)..)
or plot2d ('(realpart( ...) ..))..
RJF
_______________________________________________
Maxima mailing list
Maxima at math.utexas.edu
http://www.math.utexas.edu/mailman/listinfo/maxima