1 SUBROUTINE cbuni(Z, FNU, KODE, N, Y, NZ, NUI, NLAST, FNUL, TOL,
14 COMPLEX CSCL, CSCR, CY, RZ, ST, S1, S2, Y, Z
15 REAL ALIM, AX, AY, DFNU, ELIM, FNU, FNUI, FNUL, GNU, TOL, XX, YY,
16 * ascle, bry, str, sti, stm,
r1mach 17 INTEGER I, IFLAG, IFORM, K, KODE, N, NL, NLAST, NUI, NW, NZ
25 IF (ay.GT.ax) iform = 2
26 IF (nui.EQ.0)
GO TO 60
28 dfnu = fnu + float(n-1)
30 IF (iform.EQ.2)
GO TO 10
35 CALL cuni1(z, gnu, kode, 2, cy, nw, nlast, fnul, tol, elim, alim)
43 CALL cuni2(z, gnu, kode, 2, cy, nw, nlast, fnul, tol, elim, alim)
51 bry(1) = 1.0e+3*
r1mach(1)/tol
57 cscl =
cmplx(ax,0.0e0)
58 IF (ay.GT.bry(1))
GO TO 21
62 cscl =
cmplx(ax,0.0e0)
65 IF (ay.LT.bry(2))
GO TO 25
69 cscl =
cmplx(ax,0.0e0)
72 cscr =
cmplx(ay,0.0e0)
75 rz =
cmplx(2.0e0,0.0e0)/z
78 s2 =
cmplx(dfnu+fnui,0.0e0)*rz*s2 + s1
81 IF (iflag.GE.3)
GO TO 30
88 IF (stm.LE.ascle)
GO TO 30
95 cscl =
cmplx(ax,0.0e0)
96 cscr =
cmplx(ay,0.0e0)
107 s2 =
cmplx(fnu+fnui,0.0e0)*rz*s2 + s1
113 IF (iflag.GE.3)
GO TO 40
119 IF (stm.LE.ascle)
GO TO 40
126 cscl =
cmplx(ax,0.0e0)
127 cscr =
cmplx(ay,0.0e0)
137 IF (iform.EQ.2)
GO TO 70
142 CALL cuni1(z, fnu, kode, n, y, nw, nlast, fnul, tol, elim, alim)
150 CALL cuni2(z, fnu, kode, n, y, nw, nlast, fnul, tol, elim, alim)
152 IF (nw.LT.0)
GO TO 50
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)