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
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)