mmtu12.f90

数值计算和数值分析在Fortran下的特殊函数库,是数值计算的必备

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

! Zhang & Jin's errata applied to initialize array FC in routine FCOEF,
! and substantial reductions in routines REFINE & CVA2.
! Latest revision - 16 December 2002
! Corrections by Alan Miller (amiller @ bigpond.net.au)

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

SUBROUTINE mtu12(kf, kc, m, q, x, f1r, d1r, f2r, d2r)

!       ==============================================================
!       Purpose: Compute modified Mathieu functions of the first and
!                second kinds, Mcm(1)(2)(x,q) and Msm(1)(2)(x,q),
!                and their derivatives
!       Input:   KF --- Function code
!                       KF=1 for computing Mcm(x,q)
!                       KF=2 for computing Msm(x,q)
!                KC --- Function Code
!                       KC=1 for computing the first kind
!                       KC=2 for computing the second kind
!                            or Msm(2)(x,q) and Msm(2)'(x,q)
!                       KC=3 for computing both the first and second kinds
!                m  --- Order of Mathieu functions
!                q  --- Parameter of Mathieu functions ( q