#include <cmath>

 double cyl_neumann(double nu, double x)  // C++17

Liefert Wert der Bessel-Funktion 2. Art (auch Weber- oder Neumann-Funktion) $Y_\nu(x) = \frac{J_\nu(x)\cos(\nu\pi)-I_{-\nu} (x)}{\sin (\nu\pi)}$ bzw. $lim_{\nu\to n} Y_\nu(x)$ bei ganzzahligem $\nu$ für $x \geq 0$.

Rückgabewert: $Y_\nu(x)$.

cyl_bessel_i(), cyl_bessel_j(), cyl_bessel_k().

cyl_neumann.cpp

#include <cmath>
#include <iostream>
 
int main()
{
  std::cout << "# x n=0 n=1 n=2\n";
 
  for (int i = 0; i <= 1000; ++i)
  {
    double x = 0.01*i;
    std::cout << x 
      << '\t' << std::cyl_neumann(0, x) 
      << '\t' << std::cyl_neumann(1, x) 
      << '\t' << std::cyl_neumann(2, x) 
      << '\n';
  }
}