mrswfo.f90

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

MODULE rswfo_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."

! Latest revision - 30 December 2001
! Corrections by Alan Miller (amiller @ bigpond.net.au)
! Variable sw was used without a value assigned to it.
! There was a call to a routine with the same argument variable repeated.
! This is only permitted when the called routine does not alter either
! argument - i.e. when they are both INTENT(IN) arguments.  The duplicated
! argument was nm2 and the second one has been replaced with nm3.

IMPLICIT NONE
INTEGER, PARAMETER  :: dp = SELECTED_REAL_KIND(12, 60)
 
CONTAINS
 

SUBROUTINE rswfo(m, n, c, x, cv, kf, r1f, r1d, r2f, r2d)

!    ==========================================================
!    Purpose: Compute oblate radial functions of the first
!             and second kinds, and their derivatives
!    Input :  m  --- Mode parameter,  m = 0,1,2,...
!             n  --- Mode parameter,  n = m,m+1,m+2,...
!             c  --- Spheroidal parameter
!             x  --- Argument (x