#include "header.h"

/* GSL Bessel Function */
double
f( int n, double x)
{
    return gsl_sf_bessel_Jn( n, x);
}

/* TODO: Always returns zero for some reason!
 *
 * I have checked the function is reading in the correct
 * values and that the result isn't being changed outside of
 * this function.
 */