1 SUBROUTINE zbinu(ZR, ZI, FNU, KODE, N, CYR, CYI, NZ, RL, FNUL,
10 DOUBLE PRECISION ALIM, AZ, CWI, CWR, CYI, CYR, DFNU, ELIM, FNU,
11 * fnul, rl, tol, zeroi, zeror, zi, zr, xzabs
12 INTEGER I, INW, KODE, N, NLAST, NN, NUI, NW, NZ
14 DATA zeror,zeroi / 0.0d0, 0.0d0 /
19 dfnu = fnu + dble(float(n-1))
20 IF (az.LE.2.0d0) go
to 10
21 IF (az*az*0.25d0.GT.dfnu+1.0d0) go
to 20
26 CALL
zseri(zr, zi, fnu, kode, nn, cyr, cyi, nw, tol, elim, alim)
31 IF (nw.GE.0) go
to 120
32 dfnu = fnu + dble(float(nn-1))
34 IF (az.LT.rl) go
to 40
35 IF (dfnu.LE.1.0d0) go
to 30
36 IF (az+az.LT.dfnu*dfnu) go
to 50
41 CALL
zasyi(zr, zi, fnu, kode, nn, cyr, cyi, nw, rl, tol, elim,
43 IF (nw.LT.0) go
to 130
46 IF (dfnu.LE.1.0d0) go
to 70
51 CALL
zuoik(zr, zi, fnu, kode, 1, nn, cyr, cyi, nw, tol, elim,
53 IF (nw.LT.0) go
to 130
57 dfnu = fnu+dble(float(nn-1))
58 IF (dfnu.GT.fnul) go
to 110
59 IF (az.GT.fnul) go
to 110
61 IF (az.GT.rl) go
to 80
66 CALL
zmlri(zr, zi, fnu, kode, nn, cyr, cyi, nw, tol)
76 CALL
zuoik(zr, zi, fnu, kode, 2, 2, cwr, cwi, nw, tol, elim,
78 IF (nw.GE.0) go
to 100
86 IF (nw.GT.0) go
to 130
87 CALL
zwrsk(zr, zi, fnu, kode, nn, cyr, cyi, nw, cwr, cwi, tol,
89 IF (nw.LT.0) go
to 130
95 nui = int(sngl(fnul-dfnu)) + 1
97 CALL
zbuni(zr, zi, fnu, kode, nn, cyr, cyi, nw, nui, nlast, fnul,
99 IF (nw.LT.0) go
to 130
101 IF (nlast.EQ.0) go
to 120
subroutine zbinu(ZR, ZI, FNU, KODE, N, CYR, CYI, NZ, RL, FNUL, TOL, ELIM, ALIM)
subroutine zwrsk(ZRR, ZRI, FNU, KODE, N, YR, YI, NZ, CWR, CWI, TOL, ELIM, ALIM)
subroutine zuoik(ZR, ZI, FNU, KODE, IKFLG, N, YR, YI, NUF, TOL, ELIM, ALIM)
subroutine zbuni(ZR, ZI, FNU, KODE, N, YR, YI, NZ, NUI, NLAST, FNUL, TOL, ELIM, ALIM)
subroutine zasyi(ZR, ZI, FNU, KODE, N, YR, YI, NZ, RL, 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 zseri(ZR, ZI, FNU, KODE, N, YR, YI, NZ, TOL, ELIM, ALIM)
subroutine zmlri(ZR, ZI, FNU, KODE, N, YR, YI, NZ, TOL)