#include "Header.h"

//computes bessel function of order n at double value x
double f(int n, double x)
{
	return gsl_sf_bessel_Jn (n, x);
}