[Date Index] [Thread Index] [Author Index]

Re: Errors in Jacobian Elliptic Functions

• To: mathgroup at smc.vnet.net
• Subject: [mg7729] Re: Errors in Jacobian Elliptic Functions
• From: Paul Abbott <paul at physics.uwa.edu.au>
• Date: Wed, 2 Jul 1997 14:21:53 -0400 (EDT)
• Organization: University of Western Australia
• Sender: owner-wri-mathgroup at wolfram.com

seanross at worldnet.att.net wrote:
> 
> I have an application which requires extensive use of Jacobian Elliptic
> Functions like JacobiSN, EllipticPi and EllipticF.  I have found that
> these functions have some very small errors near the ends of their
> ranges.  For example, JacobiSN is supposed to be bounded by +1 and -1,
> just like Sine and Cosine.  However, for JacobiSN[x,p], with p very
> close to zero, JacobiSN can occasionally come out slightly greater than
> +1.  

I have not been able to reproduce this behaviour.  However I think I
know what you are encountering -- and the general solution: you probably
need to use arbitrary precision arithmetic.  There is an example in the
Mathematica Journal about computing Legendre functions of large order. 
One gets garbage there unless high precision is used.

Cheers,
	Paul 

____________________________________________________________________ 
Paul Abbott				      Phone: +61-8-9380-2734
Department of Physics                           Fax: +61-8-9380-1014
The University of Western Australia           
Nedlands WA  6907                     mailto:paul at physics.uwa.edu.au 
AUSTRALIA                              http://www.pd.uwa.edu.au/Paul

            God IS a weakly left-handed dice player
____________________________________________________________________