#include <gsl/gsl_sf_bessel.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>

int main(void) {
  int i;
  double x, y;
  FILE *out;
 
  out=fopen("results.dat","w");
 
  for(i=0;i<400;i++) {
    x=pow(10,i/10.0);
    y=gsl_sf_bessel_jl(0,x);
    fprintf(out,"%e % e\n",x,y);
  }

  fclose(out);
  return 0;
}
