mlpmv.f90

来自「数值计算和数值分析在Fortran下的特殊函数库,是数值计算的必备」· F90 代码 · 共 25 行

MODULE lpmv_func
 
! From the book "Computation of Special Functions"
!      by Shanjie Zhang and Jianming Jin
!   Copyright 1996 by John Wiley & Sons, Inc.
! The authors state:
!   "However, we give permission to the reader who purchases this book
!    to incorporate any of these programs into his or her programs
!    provided that the copyright is acknowledged."
 
IMPLICIT NONE
INTEGER, PARAMETER  :: dp = SELECTED_REAL_KIND(12, 60)

REAL (dp), PARAMETER  :: pi = 3.141592653589793_dp, el = .5772156649015329_dp

CONTAINS
 

SUBROUTINE lpmv(v, m, x, pmv)

!    =======================================================
!    Purpose: Compute the associated Legendre function
!             Pmv(x) with an integer order and an arbitrary
!             nonnegative degree v
!    Input :  x   --- Argument of Pm(x)  ( -1