[Help-gsl] regular coulomb wave function array failing

[Andrew Spott]

I'm attempting to find the coulomb wave functions from l = 0 to l = 50 for the 
points between x = 0 and x = 1000 for k's of ~0 (I get a bunch of underflow 
values for k < .0099) to k = 1:

gsl_sf_coulomb_wave_F_array(0, lmax, 1. /k , k * j * dr , values, exponent)

where j goes from 0 to 10001, dr = .1, k goes from  ~ .01 to 1 (in steps of .01 
ish).

However, for certain regions in k * j * dr, I get loss of accuracy errors, and 
I don't know how to make this function iterate for longer to make those errors 
go away.

Any ideas (I've attached the error log out of my debug program which shows the 
beginning and end of the loss of accuracy errors, and were they happen(in terms 
of k * x = k * j * dr, and k).)

-Andrew

errors
Description: Binary data