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