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

Bug in StruveH

• To: mathgroup at yoda.ncsa.uiuc.edu
• Subject: Bug in StruveH
• From: pjanhune at finsun.csc.fi
• Date: Wed, 23 Jan 91 09:37:48 +0200

I dont know how to contact R.Maeder, so I put the report here.

[This message has been emailed to Maeder for future response. -smc]

The package StruveH.m from the book by Roman Maeder contains the
following bug:

When called with an integer or rational second argument, the function
StruveH returns an exact-looking result, which is WRONG!! It's only
the first term of the series expansion. For floating-point arguments
or complex numbers the function behaves correctly. One cure is to insert
the line

If[prec == Infinity, prec = Precision[1.]];

just before the line

zf = (z/2)^(nu+1); g1 = Gamma[k1]; g2 = Gamma[k2];

in function StruveH.