Click above to get retro games delivered to your door ever month!

bessel()                 Return Result of a Bessel Function of x

 #include <math.h>

 double j0(x);
 double x;                            Floating-point value

 double j1(x);
 double x;                            Floating-point value

 double jn(n,x);
 int n;                               Integer order
 double x;                            Floating-point value

 double y0(x);
 double x;                            Floating-point value

 double y1(x);
 double x;                            Floating-point value

 double yn(n,x);
 int n;                               Integer order
 double x;                            Floating-point value


    The "bessel routine" comprises the six functions j0(),j1(), jn(),
    y0(), y1(), and yn().  The j0(), j1(), and jn() routines return
    Bessel functions of the first kind of orders--0, 1, and n,
    respectively.

    The y0(), y1(), and yn() routines return Bessel functions of the
    second kind of orders, 0, 1, and n respectively.

    Returns:    The result of a Bessel function of 'x'.  For y0(), y1(),
                and yn(), if 'x' is negative, 'errno' is set to EDOM, a
                DOMAIN error message is printed to 'stderr', and the
                value of negative HUGE is returned.

      Notes:    'x' must be positive.

                Error handling can be modified by using the matherr()
                routine.

  -------------------------------- Example ---------------------------------

    The following statements return the Bessel function of 'x'.

          #include <math.h>

          double x, y, z;

          main()
          {
           y = j0(x);
           z = yn(3,x);
           .
           .
          }