I found 2 missprints in gsl-ref.ps regarding Coulomb Wave Functions gsl-ref.ps must be, see gsl_sf_coulomb.h int gsl_sf_coulomb_wave_FG (double.... gsl_sf_coulomb_wave_FG_e int gsl_sf_coulomb_wave_CL (double .... gsl_sf_coulomb_wave_CL_e --Vlad