[zr,zi,fnu,kode,mr,n,yr,yi,nz,tol,elim,alim]=zbunk(zr,zi,fnu,kode,mr,n,yr,yi,nz,tol,elim,alim);

function [zr,zi,fnu,kode,mr,n,yr,yi,nz,tol,elim,alim]=zbunk(zr,zi,fnu,kode,mr,n,yr,yi,nz,tol,elim,alim);
%***BEGIN PROLOGUE  ZBUNK
%***SUBSIDIARY
%***PURPOSE  Subsidiary to ZBESH and ZBESK
%***LIBRARY   SLATEC
%***TYPE      ALL (CBUNI-A, ZBUNI-A)
%***AUTHOR  Amos, D. E., (SNL)
%***DESCRIPTION
%
%     ZBUNK COMPUTES THE K BESSEL FUNCTION FOR FNU.GT.FNUL.
%     ACCORDING TO THE UNIFORM ASYMPTOTIC EXPANSION FOR K(FNU,Z)
%     IN ZUNK1 AND THE EXPANSION FOR H(2,FNU,Z) IN ZUNK2
%
%***SEE ALSO  ZBESH, ZBESK
%***ROUTINES CALLED  ZUNK1, ZUNK2
%***REVISION HISTORY  (YYMMDD)
%   830501  DATE WRITTEN
%   910415  Prologue converted to Version 4.0 format.  (BAB)
%***end PROLOGUE  ZBUNK
%     COMPLEX Y,Z
persistent ax ay ; 

if isempty(ax), ax=0; end;
if isempty(ay), ay=0; end;
%***FIRST EXECUTABLE STATEMENT  ZBUNK
nz = 0;
ax = abs(zr).*1.7321d0;
ay = abs(zi);
if( ay>ax )
%-----------------------------------------------------------------------
%     ASYMPTOTIC EXPANSION FOR H(2,FNU,Z*EXP(M*HPI)) FOR LARGE FNU
%     APPLIED IN PI/3.LT.ABS(ARG(Z)).LE.PI/2 WHERE M=+I OR -I
%     AND HPI=PI/2
%-----------------------------------------------------------------------
[zr,zi,fnu,kode,mr,n,yr,yi,nz,tol,elim,alim]=zunk2(zr,zi,fnu,kode,mr,n,yr,yi,nz,tol,elim,alim);
else;
%-----------------------------------------------------------------------
%     ASYMPTOTIC EXPANSION FOR K(FNU,Z) FOR LARGE FNU APPLIED IN
%     -PI/3.LE.ARG(Z).LE.PI/3
%-----------------------------------------------------------------------
[zr,zi,fnu,kode,mr,n,yr,yi,nz,tol,elim,alim]=zunk1(zr,zi,fnu,kode,mr,n,yr,yi,nz,tol,elim,alim);
end;
end
%DECK ZDIV