[Date Index] [Thread Index] [Author Index]
Bug in StruveH
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.