************80
! FCOEF: expansion coefficients for Mathieu and modified Mathieu functions.
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:
01 August 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 ) KD, the case code.
1, for cem(x,q) ( m = 0,2,4,...)
2, for cem(x,q) ( m = 1,3,5,...)
3, for sem(x,q) ( m = 1,3,5,...)
4, for sem(x,q) ( m = 2,4,6,...)
Input, integer ( kind = 4 ) M, the order of the Mathieu function.
Input, real ( kind = 8 ) Q, the parameter of the Mathieu functions.
Input, real ( kind = 8 ) A, the characteristic value of the Mathieu
functions for given m and q.
Output, real ( kind = 8 ) FC(*), the expansion coefficients of Mathieu
functions ( k = 1,2,...,KM ). FC(1),FC(2),FC(3),... correspond to
A0,A2,A4,... for KD = 1 case,
A1,A3,A5,... for KD = 2 case,
B1,B3,B5,... for KD = 3 case,
B2,B4,B6,... for KD = 4 case.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer(kind=4) | :: | kd | ||||
| integer(kind=4) | :: | m | ||||
| real(kind=8) | :: | q | ||||
| real(kind=8) | :: | a | ||||
| real(kind=8) | :: | fc(251) |