************80
! CYZO computes zeros of complex Bessel functions Y0(z) and Y1(z) and Y1'(z).
Parameters:
Ths procedure computes the complex zeros of Y0(z), Y1(z) and Y1'(z),
and their associated values at the zeros using the modified Newton's
iteration method.
Licensing:
This routine is copyrighted by Shanjie Zhang and Jianming Jin. However,
they give permission to incorporate this routine into a user program
provided that the copyright is acknowledged.
Modified:
22 July 2012
Author:
Shanjie Zhang, Jianming Jin
Reference:
Shanjie Zhang, Jianming Jin,
Computation of Special Functions,
Wiley, 1996,
ISBN: 0-471-11963-6,
LC: QA351.C45.
Parameters:
Input, integer ( kind = 4 ) NT, the number of zeros.
Input, integer ( kind = 4 ) KF, the function choice.
0 for Y0(z) and Y1(z0);
1 for Y1(z) and Y0(z1);
2 for Y1'(z) and Y1(z1').
Input, integer ( kind = 4 ) KC, complex/real choice.
0, for complex roots;
1, for real roots.
Output, real ( kind = 8 ) ZO(NT), ZV(NT), the zeros of Y0(z) or Y1(z)
or Y1'(z), and the value of Y0'(z) or Y1'(z) or Y1(z) at the L-th zero.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer(kind=4) | :: | nt | ||||
| integer(kind=4) | :: | kf | ||||
| integer(kind=4) | :: | kc | ||||
| complex(kind=8) | :: | zo(nt) | ||||
| complex(kind=8) | :: | zv(nt) |