bug in function airy_i in file hypergeom.c

[Peter <petsie@dordos.net>]

Dear developers,

I am absolutely new to the internals of the code. Therefore I can just
describe the misbehavior of airy_i:

the test file for hypergeometric functions contains the instruction
"airy(-100)[2]" which should return

airybi(-100) ~= 0.02427... but it returns 0.01213...

This might be caused by a malfunction of lines 74 - 80 in hypergeom.c


I get the wrong result from above if i assume -100>0 and set
a=b=sqrt(abs(x))

wit the setting for a and b from lines 78,79 I get the correct result.


Strange enough: results of airy(2) and airy(-1) are correct.


Kind regards,

Peter Pein