#include <gsl/gsl_sf.h>
#include <stdio.h>

int main()
{
	gsl_sf_result f, fp, g, gp;
	double exp_f, exp_g;
	int i;
	i = gsl_sf_coulomb_wave_FG_e(0.,1.2693881947287221e-07, 37, 1,
		&f,&fp, & g,&gp,  &exp_f, &exp_g);

	printf(" f=(%g, %g)\n"
		" fp=(%g, %g)\n"
		" g=(%g, %g)\n"
		" gp=(%g, %g)\n"
		" exp_f=%g exp_g=%g\n",
		f.val, f.err,
		fp.val, fp.err, 
		g.val, g.err, 
		gp.val, gp.err, 
		exp_f, exp_g
		);
	if (i) printf("%i\n", i);
	return i;
}