lqmn Subroutine

public subroutine lqmn(mm, m, n, x, qm, qd)

************80

! LQMN computes associated Legendre functions Qmn(x) and derivatives.

Discussion:

This routine computes the associated Legendre functions of the
second kind, Qmn(x) and Qmn'(x).

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:

13 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 ) MM, determines the leading dimension 
of QM and QD.

Input, integer ( kind = 4 ) M, the order of Qmn(x).

Input, integer ( kind = 4 ) N, the degree of Qmn(x).

Output, real ( kind = 8 ) QM(0:MM,0:N), QD(0:MM,0:N), contains the values
of Qmn(x) and Qmn'(x).

Arguments

Type IntentOptional Attributes Name
integer(kind=4) :: mm
integer(kind=4) :: m
integer(kind=4) :: n
real(kind=8) :: x
real(kind=8) :: qm(0:mm,0:n)
real(kind=8) :: qd(0:mm,0:n)