zunik.f
SUBROUTINE ZUNIK (ZRR, ZRI, FNU, IKFLG, IPMTR, TOL, INIT, PHIR,
+ PHII, ZETA1R, ZETA1I, ZETA2R, ZETA2I, SUMR, SUMI, CWRKR, CWRKI)
C***BEGIN PROLOGUE ZUNIK
....
....
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 (CUNIK-A, ZUNIK-A)
C***AUTHOR Amos, D. E., (SNL)
C***DESCRIPTION
C
C ZUNIK COMPUTES PARAMETERS FOR THE UNIFORM ASYMPTOTIC
C EXPANSIONS OF THE I AND K FUNCTIONS ON IKFLG= 1 OR 2
C RESPECTIVELY BY
C
C W(FNU,ZR) = PHI*EXP(ZETA)*SUM
C
C WHERE ZETA=-ZETA1 + ZETA2 OR
C ZETA1 - ZETA2
C
C THE FIRST CALL MUST HAVE INIT=0. SUBSEQUENT CALLS WITH THE
C SAME ZR AND FNU WILL RETURN THE I OR K FUNCTION ON IKFLG=
C 1 OR 2 WITH NO CHANGE IN INIT. CWRK IS A COMPLEX WORK
C ARRAY. IPMTR=0 COMPUTES ALL PARAMETERS. IPMTR=1 COMPUTES PHI,
C ZETA1,ZETA2.
C
C***SEE ALSO ZBESI, ZBESK
C***ROUTINES CALLED D1MACH, ZDIV, ZLOG, ZSQRT
C***REVISION HISTORY (YYMMDD)
C 830501 DATE WRITTEN
C 910415 Prologue converted to Version 4.0 format. (BAB)
C 930122 Added EXTERNAL statement with ZLOG and ZSQRT. (RWC)
C***END PROLOGUE ZUNIK