1 SUBROUTINE cbinu(Z, FNU, KODE, N, CY, NZ, RL, FNUL, TOL, ELIM,
10 COMPLEX CW, CY, CZERO, Z
11 REAL ALIM, AZ, DFNU, ELIM, FNU, FNUL, RL, TOL
12 INTEGER I, INW, KODE, N, NLAST, NN, NUI, NW, NZ
14 DATA czero / (0.0e0,0.0e0) /
19 dfnu = fnu + float(n-1)
20 IF (az.LE.2.0e0) go
to 10
21 IF (az*az*0.25e0.GT.dfnu+1.0e0) go
to 20
26 CALL
cseri(z, fnu, kode, nn, cy, nw, tol, elim, alim)
31 IF (nw.GE.0) go
to 120
32 dfnu = fnu + float(nn-1)
34 IF (az.LT.rl) go
to 40
35 IF (dfnu.LE.1.0e0) go
to 30
36 IF (az+az.LT.dfnu*dfnu) go
to 50
41 CALL
casyi(z, fnu, kode, nn, cy, nw, rl, tol, elim, alim)
42 IF (nw.LT.0) go
to 130
45 IF (dfnu.LE.1.0e0) go
to 70
50 CALL
cuoik(z, fnu, kode, 1, nn, cy, nw, tol, elim, alim)
51 IF (nw.LT.0) go
to 130
55 dfnu = fnu+float(nn-1)
56 IF (dfnu.GT.fnul) go
to 110
57 IF (az.GT.fnul) go
to 110
59 IF (az.GT.rl) go
to 80
64 CALL
cmlri(z, fnu, kode, nn, cy, nw, tol)
74 CALL
cuoik(z, fnu, kode, 2, 2, cw, nw, tol, elim, alim)
75 IF (nw.GE.0) go
to 100
82 IF (nw.GT.0) go
to 130
83 CALL
cwrsk(z, fnu, kode, nn, cy, nw, cw, tol, elim, alim)
84 IF (nw.LT.0) go
to 130
90 nui = int(fnul-dfnu) + 1
92 CALL
cbuni(z, fnu, kode, nn, cy, nw, nui, nlast, fnul, tol, elim,
94 IF (nw.LT.0) go
to 130
96 IF (nlast.EQ.0) go
to 120
subroutine cseri(Z, FNU, KODE, N, Y, NZ, TOL, ELIM, ALIM)
subroutine cwrsk(ZR, FNU, KODE, N, Y, NZ, CW, TOL, ELIM, ALIM)
subroutine cbuni(Z, FNU, KODE, N, Y, NZ, NUI, NLAST, FNUL, TOL, ELIM, ALIM)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)
subroutine cuoik(Z, FNU, KODE, IKFLG, N, Y, NUF, TOL, ELIM, ALIM)
subroutine casyi(Z, FNU, KODE, N, Y, NZ, RL, TOL, ELIM, ALIM)
subroutine cmlri(Z, FNU, KODE, N, Y, NZ, TOL)
subroutine cbinu(Z, FNU, KODE, N, CY, NZ, RL, FNUL, TOL, ELIM, ALIM)