Std::cyl_bessel_j,std::cyl_bessel_jf,std::cyl_bessel_jl (3) Linux Manual Page

std::cyl_bessel_j,std::cyl_bessel_jf,std::cyl_bessel_jl – std::cyl_bessel_j,std::cyl_bessel_jf,std::cyl_bessel_jl

Synopsis

double cyl_bessel_j( double ν, double x );

float cyl_bessel_jf( float ν, float x );                  (1) (since C++17)

long double cyl_bessel_jl( long double ν, long double x );

Promoted cyl_bessel_j( Arithmetic ν, Arithmetic x );      (2) (since C++17)

1) Computes the cylindrical_Bessel_function_of_the_first_kind of ν and x.
2) A set of overloads or a function template for all combinations of arguments of arithmetic type not covered by (1). If any argument has integral_type, it is cast to double. If any argument is long double, then the return type Promoted is also long double, otherwise the return type is always double.

Parameters

ν- the order of the function
x – the argument of the function

Return value

If no errors occur, value of the cylindrical Bessel function of the first kind of ν and x, that is J
ν(x) = Σ∞
k=0
(-1)k
(x/2)ν+2k
k!Γ(ν+k+1)
(for x≥0), is returned.

Error handling

Errors may be reported as specified in math_errhandling

* If the argument is NaN, NaN is returned and domain error is not reported
* If ν>=128, the behavior is implementation-defined

Notes

Implementations that do not support C++17, but support ISO_29124:2010, provide this function if __STDCPP_MATH_SPEC_FUNCS__ is defined by the implementation to a value at least 201003L and if the user defines __STDCPP_WANT_MATH_SPEC_FUNCS__ before including any standard library headers.
Implementations that do not support ISO 29124:2010 but support TR 19768:2007 (TR1), provide this function in the header tr1/cmath and namespace std::tr1
An implementation of this function is also available_in_boost.math

Example

// Run this code

#include <cmath>

#include <iostream>

int main()

{
    // spot check for ν == 0
    double x = 1.2345;
    std::cout << "J_0(" << x << ") = " << std::cyl_bessel_j(0, x) << '\n';
    // series expansion for J_0
    double fct = 1;
    double sum = 0;
    for (int k = 0; k < 6; fct *= ++k) {
        sum += std::pow(-1, k) * std::pow((x / 2), 2 * k) / std::pow(fct, 2);
        std::cout << "sum = " << sum << '\n';
    }
}

Output:

  J_0(1.2345) = 0.653792
  sum = 1
  sum = 0.619002
  sum = 0.655292
  sum = 0.653756
  sum = 0.653793
  sum = 0.653792

External links

Weisstein,_Eric_W._"Bessel_Function_of_the_First_Kind." From MathWorld–A Wolfram Web Resource.

std::cyl_bessel_j,std::cyl_bessel_jf,std::cyl_bessel_jl (3) Linux Manual Page