#include "header.h" double f(int n,double x) { double ans; ans=gsl_sf_bessel_Jn(n,x); return ans; }