************80
! OTHPL computes orthogonal polynomials Tn(x), Un(x), Ln(x) or Hn(x).
Discussion:
This procedure computes orthogonal polynomials: Tn(x) or Un(x),
or Ln(x) or Hn(x), and their derivatives.
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:
08 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 ) KT, the function code:
1 for Chebyshev polynomial Tn(x)
2 for Chebyshev polynomial Un(x)
3 for Laguerre polynomial Ln(x)
4 for Hermite polynomial Hn(x)
Input, integer ( kind = 4 ) N, the order.
Input, real ( kind = 8 ) X, the argument.
Output, real ( kind = 8 ) PL(0:N), DPL(0:N), the value and derivative of
the polynomials of order 0 through N at X.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer(kind=4) | :: | kf | ||||
| integer | :: | n | ||||
| real(kind=8) | :: | x | ||||
| real(kind=8) | :: | pl(0:n) | ||||
| real(kind=8) | :: | dpl(0:n) |