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.