#include "header.h" // func(int n,double x) is a function that evaluates the bessel function of order n given a double value x double func(int n,double x) { double ans=gsl_sf_bessel_Jn (n,x); // bessel function in GSL library return ans; }