GNU Octave  6.2.0
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
cbunk.f
Go to the documentation of this file.
1  SUBROUTINE cbunk(Z, FNU, KODE, MR, N, Y, NZ, TOL, ELIM, ALIM)
2 C***BEGIN PROLOGUE CBUNK
3 C***REFER TO CBESK,CBESH
4 C
5 C CBUNK COMPUTES THE K BESSEL FUNCTION FOR FNU.GT.FNUL.
6 C ACCORDING TO THE UNIFORM ASYMPTOTIC EXPANSION FOR K(FNU,Z)
7 C IN CUNK1 AND THE EXPANSION FOR H(2,FNU,Z) IN CUNK2
8 C
9 C***ROUTINES CALLED CUNK1,CUNK2
10 C***END PROLOGUE CBUNK
11  COMPLEX Y, Z
12  REAL ALIM, AX, AY, ELIM, FNU, TOL, XX, YY
13  INTEGER KODE, MR, N, NZ
14  dimension y(n)
15  nz = 0
16  xx = real(z)
17  yy = aimag(z)
18  ax = abs(xx)*1.7321e0
19  ay = abs(yy)
20  IF (ay.GT.ax) GO TO 10
21 C-----------------------------------------------------------------------
22 C ASYMPTOTIC EXPANSION FOR K(FNU,Z) FOR LARGE FNU APPLIED IN
23 C -PI/3.LE.ARG(Z).LE.PI/3
24 C-----------------------------------------------------------------------
25  CALL cunk1(z, fnu, kode, mr, n, y, nz, tol, elim, alim)
26  GO TO 20
27  10 CONTINUE
28 C-----------------------------------------------------------------------
29 C ASYMPTOTIC EXPANSION FOR H(2,FNU,Z*EXP(M*HPI)) FOR LARGE FNU
30 C APPLIED IN PI/3.LT.ABS(ARG(Z)).LE.PI/2 WHERE M=+I OR -I
31 C AND HPI=PI/2
32 C-----------------------------------------------------------------------
33  CALL cunk2(z, fnu, kode, mr, n, y, nz, tol, elim, alim)
34  20 CONTINUE
35  RETURN
36  END
subroutine cbunk(Z, FNU, KODE, MR, N, Y, NZ, TOL, ELIM, ALIM)
Definition: cbunk.f:2
subroutine cunk1(Z, FNU, KODE, MR, N, Y, NZ, TOL, ELIM, ALIM)
Definition: cunk1.f:2
subroutine cunk2(Z, FNU, KODE, MR, N, Y, NZ, TOL, ELIM, ALIM)
Definition: cunk2.f:2
ColumnVector real(const ComplexColumnVector &a)
Definition: dColVector.cc:137
static T abs(T x)
Definition: pr-output.cc:1678