double BesselJ0(double x);
double BesselJm(double x, int m);