Bessel Function – C++

				
					#include <iostream>
#include <cmath>

int main() {
    double x = 2.5;
    double j0_result = std::cyl_bessel_j(0, x); // Bessel function of the first kind of order zero
    double j1_result = std::cyl_bessel_j(1, x); // Bessel function of the first kind of order one
    double yn_result = std::cyl_neumann(2, x); // Bessel function of the second kind of order two

    std::cout << "J0(" << x << ") = " << j0_result << std::endl;
    std::cout << "J1(" << x << ") = " << j1_result << std::endl;
    std::cout << "Y2(" << x << ") = " << yn_result << std::endl;

    return 0;
}