core/vnl/vnl_bessel.cxx File Reference

Bessel functions of the first kind. More...

#include "vnl_bessel.h"
#include <vcl_algorithm.h>

Go to the source code of this file.

Functions
void	vnl_bessel (unsigned n_max, double x, vnl_vector< double > &J)
	Compute Bessel functions of first kind up to order n_max.
double	vnl_bessel0 (double x)
	Returns J_0(x), the value of the Bessel function of order 0 at x.
double	vnl_bessel (unsigned n, double x)
	Returns J_n(x), the value of the Bessel function of order n at x.

---

Detailed Description

Bessel functions of the first kind.

Author:

Tim Cootes

Definition in file vnl_bessel.cxx.

---

Function Documentation

void vnl_bessel	(	unsigned	n_max,
		double	x,
		vnl_vector< double > &	J
	)

Compute Bessel functions of first kind up to order n_max.

On exit, J[i] = J_i(x) for i=0..n_max

Uses recurrence relation: J_(n-1)(x)+J_(n+1)=(2n/x)J_n(x) Thus J_n(x) = (2(n+1)/x)J_(n+1)(x) - J_(n+2)(x) Start with arbitrary guess for high n and work backwards. Normalise suitably.

Definition at line 17 of file vnl_bessel.cxx.

double vnl_bessel	(	unsigned	n,
		double	x
	)

Returns J_n(x), the value of the Bessel function of order n at x.

Bessel function of the first kind of order zero.

Uses recurrence relation: J_(n-1)(x)+J_(n+1)=(2n/x)J_n(x)

Definition at line 69 of file vnl_bessel.cxx.

double vnl_bessel0

(

double

x

)

Returns J_0(x), the value of the Bessel function of order 0 at x.

Bessel function of the first kind of order zero.

Uses recurrence relation: J_(n-1)(x)+J_(n+1)=(2n/x)J_n(x)

Definition at line 45 of file vnl_bessel.cxx.