d2/d7a/cbunk_8f_source.html

       SUBROUTINE cbunk(Z, FNU, KODE, MR, N, Y, NZ, TOL, ELIM, ALIM)

 C***BEGIN PROLOGUE  CBUNK

 C***REFER TO  CBESK,CBESH

 C

 C     CBUNK COMPUTES THE K BESSEL FUNCTION FOR FNU.GT.FNUL.

 C     ACCORDING TO THE UNIFORM ASYMPTOTIC EXPANSION FOR K(FNU,Z)

 C     IN CUNK1 AND THE EXPANSION FOR H(2,FNU,Z) IN CUNK2

 C

 C***ROUTINES CALLED  CUNK1,CUNK2

 C***END PROLOGUE  CBUNK

       COMPLEX Y, Z

       REAL ALIM, AX, AY, ELIM, FNU, TOL, XX, YY

       INTEGER KODE, MR, N, NZ

       dimension y(n)

       nz = 0

       xx = REAL(z)

       yy = aimag(z)

       ax = abs(xx)*1.7321e0

       ay = abs(yy)

       IF (ay.GT.ax) go to 10

 C-----------------------------------------------------------------------

 C     ASYMPTOTIC EXPANSION FOR K(FNU,Z) FOR LARGE FNU APPLIED IN

 C     -PI/3.LE.ARG(Z).LE.PI/3

 C-----------------------------------------------------------------------

       CALL cunk1(z, fnu, kode, mr, n, y, nz, tol, elim, alim)

       go to 20

    10 CONTINUE

 C-----------------------------------------------------------------------

 C     ASYMPTOTIC EXPANSION FOR H(2,FNU,Z*EXP(M*HPI)) FOR LARGE FNU

 C     APPLIED IN PI/3.LT.ABS(ARG(Z)).LE.PI/2 WHERE M=+I OR -I

 C     AND HPI=PI/2

 C-----------------------------------------------------------------------

       CALL cunk2(z, fnu, kode, mr, n, y, nz, tol, elim, alim)

    20 CONTINUE

       RETURN

       END

cunk1
subroutine cunk1(Z, FNU, KODE, MR, N, Y, NZ, TOL, ELIM, ALIM)
Definition: cunk1.f:1

cbunk
subroutine cbunk(Z, FNU, KODE, MR, N, Y, NZ, TOL, ELIM, ALIM)
Definition: cbunk.f:1

cunk2
subroutine cunk2(Z, FNU, KODE, MR, N, Y, NZ, TOL, ELIM, ALIM)
Definition: cunk2.f:1

to
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
Definition: Quad-opts.cc:233

dimension
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)

abs
OCTAVE_EXPORT octave_value_list return the value of the option it must match the dimension of the state and the relative tolerance must also be a vector of the same length tem it must match the dimension of the state and the absolute tolerance must also be a vector of the same length The local error test applied at each integration step is xample roup abs(local error in x(i))<