00001 // This is core/vnl/vnl_bessel.h 00002 #ifndef vnl_bessel_h_ 00003 #define vnl_bessel_h_ 00004 //: 00005 // \file 00006 // \brief Bessel functions of the first kind 00007 // \author Tim Cootes 00008 00009 #include <vnl/vnl_vector.h> 00010 00011 //: Returns J_n(x), the value of the Bessel function of order n at x 00012 // Bessel function of the first kind of order n 00013 double vnl_bessel(unsigned n, double x); 00014 00015 //: Returns J_0(x), the value of the Bessel function of order 0 at x 00016 // Bessel function of the first kind of order zero 00017 double vnl_bessel0(double x); 00018 00019 //: Compute Bessel functions of first kind up to order n_max 00020 // On exit, J[i] = J_i(x) for i=0..n_max 00021 void vnl_bessel(unsigned n_max, double x, vnl_vector<double>& J); 00022 00023 #endif // vnl_bessel_h_