#include "header.h"

double f(int n, double x)
{
  double tempbessel=gsl_sf_bessel_Jn (n,x);		//use a temperary variable to store the output for the bessel function
  return tempbessel;					//return the value for the bessel function
}