1 SUBROUTINE zwrsk(ZRR, ZRI, FNU, KODE, N, YR, YI, NZ, CWR, CWI,
12 DOUBLE PRECISION ACT, ACW, ALIM, ASCLE, CINUI, CINUR, CSCLR, CTI,
13 * ctr, cwi, cwr, c1i, c1r, c2i, c2r, elim, fnu, pti, ptr, ract,
14 * sti, str, tol, yi, yr, zri, zrr, xzabs,
d1mach 15 INTEGER I, KODE, N, NW, NZ
23 CALL zbknu(zrr, zri, fnu, kode, 2, cwr, cwi, nw, tol, elim, alim)
25 CALL zrati(zrr, zri, fnu, n, yr, yi, tol)
32 IF (kode.EQ.1)
GO TO 10
42 acw = xzabs(cwr(2),cwi(2))
43 ascle = 1.0d+3*
d1mach(1)/tol
45 IF (acw.GT.ascle)
GO TO 20
50 IF (acw.LT.ascle)
GO TO 30
63 ptr = str*c1r - sti*c1i
64 pti = str*c1i + sti*c1r
67 ctr = zrr*ptr - zri*pti
68 cti = zrr*pti + zri*ptr
75 cinur = ptr*ctr - pti*cti
76 cinui = ptr*cti + pti*ctr
81 ptr = str*cinur - sti*cinui
82 cinui = str*cinui + sti*cinur
double precision function d1mach(i)
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)