zunhj.f
SUBROUTINE ZUNHJ (ZR, ZI, FNU, IPMTR, TOL, PHIR, PHII, ARGR, ARGI,
+ ZETA1R, ZETA1I, ZETA2R, ZETA2I, ASUMR, ASUMI, BSUMR, BSUMI)
C***BEGIN PROLOGUE ZUNHJ
....
....
Warning: this routine is not intended to be user-callable.
....
....
C***SUBSIDIARY
C***PURPOSE Subsidiary to ZBESI and ZBESK
C***LIBRARY SLATEC
C***TYPE ALL (CUNHJ-A, ZUNHJ-A)
C***AUTHOR Amos, D. E., (SNL)
C***DESCRIPTION
C
C REFERENCES
C HANDBOOK OF MATHEMATICAL FUNCTIONS BY M. ABRAMOWITZ AND I.A.
C STEGUN, AMS55, NATIONAL BUREAU OF STANDARDS, 1965, CHAPTER 9.
C
C ASYMPTOTICS AND SPECIAL FUNCTIONS BY F.W.J. OLVER, ACADEMIC
C PRESS, N.Y., 1974, PAGE 420
C
C ABSTRACT
C ZUNHJ COMPUTES PARAMETERS FOR BESSEL FUNCTIONS C(FNU,Z) =
C J(FNU,Z), Y(FNU,Z) OR H(I,FNU,Z) I=1,2 FOR LARGE ORDERS FNU
C BY MEANS OF THE UNIFORM ASYMPTOTIC EXPANSION
C
C C(FNU,Z)=C1*PHI*( ASUM*AIRY(ARG) + C2*BSUM*DAIRY(ARG) )
C
C FOR PROPER CHOICES OF C1, C2, AIRY AND DAIRY WHERE AIRY IS
C AN AIRY FUNCTION AND DAIRY IS ITS DERIVATIVE.
C
C (2/3)*FNU*ZETA**1.5 = ZETA1-ZETA2,
C
C ZETA1=0.5*FNU*CLOG((1+W)/(1-W)), ZETA2=FNU*W FOR SCALING
C PURPOSES IN AIRY FUNCTIONS FROM CAIRY OR CBIRY.
C
C MCONJ=SIGN OF AIMAG(Z), BUT IS AMBIGUOUS WHEN Z IS REAL AND
C MUST BE SPECIFIED. IPMTR=0 RETURNS ALL PARAMETERS. IPMTR=
C 1 COMPUTES ALL EXCEPT ASUM AND BSUM.
C
C***SEE ALSO ZBESI, ZBESK
C***ROUTINES CALLED D1MACH, ZABS, ZDIV, ZLOG, ZSQRT
C***REVISION HISTORY (YYMMDD)
C 830501 DATE WRITTEN
C 910415 Prologue converted to Version 4.0 format. (BAB)
C 930122 Added ZLOG and ZSQRT to EXTERNAL statement. (RWC)
C***END PROLOGUE ZUNHJ