mjdzo.f90

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

MODULE jdzo_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 the lower limits of 2 arrays (ZO & ZOC).
! Latest revision - 5 January 2002
! Corrections by Alan Miller (amiller @ bigpond.net.au)
! Variables x, zoc & p were used before values had been assigned to them.

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


SUBROUTINE jdzo(nt, n, m, p, zo)

!       ===========================================================
!       Purpose: Compute the zeros of Bessel functions Jn(x) and
!                Jn'(x), and arrange them in the order of their
!                magnitudes
!       Input :  NT    --- Number of total zeros ( NT