#include<stdio.h>
#include<tgmath.h>
#include<gsl/gsl_sf_bessel.h>

int main(void){

	double x, a=0, dx=0.1, b=10;
	for(x=a;x<=b;x+=dx)
		printf( "%lf %lf %lf\n", x, gsl_sf_bessel_J0(x), gsl_sf_bessel_J1(x) );

	x=0.1;
	fprintf(stderr,"%g %g %g\n",x,gsl_sf_bessel_Jn(2,x),x*x/8);
	fprintf(stderr,"%g %g %g\n",x,gsl_sf_bessel_Yn(2,x),-4/x/x/M_PI);

return 0;
}