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
