************80
! CIK01: modified Bessel I0(z), I1(z), K0(z) and K1(z) for complex argument.
Discussion:
This procedure computes the modified Bessel functions I0(z), I1(z),
K0(z), K1(z), and their derivatives for a complex argument.
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:
31 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, complex ( kind = 8 ) Z, the argument.
Output, complex ( kind = 8 ) CBI0, CDI0, CBI1, CDI1, CBK0, CDK0, CBK1,
CDK1, the values of I0(z), I0'(z), I1(z), I1'(z), K0(z), K0'(z), K1(z),
and K1'(z).
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
complex(kind=8) | :: | z | ||||
complex(kind=8) | :: | cbi0 | ||||
complex(kind=8) | :: | cdi0 | ||||
complex(kind=8) | :: | cbi1 | ||||
complex(kind=8) | :: | cdi1 | ||||
complex(kind=8) | :: | cbk0 | ||||
complex(kind=8) | :: | cdk0 | ||||
complex(kind=8) | :: | cbk1 | ||||
complex(kind=8) | :: | cdk1 |