1 SUBROUTINE cwrsk(ZR, FNU, KODE, N, Y, NZ, CW, TOL, ELIM, ALIM)
10 COMPLEX CINU, CSCL, CT, CW, C1, C2, RCT, ST, Y, ZR
11 REAL ACT, ACW, ALIM, ASCLE, ELIM, FNU, S1, S2, TOL, YY
12 INTEGER I, KODE, N, NW, NZ
20 CALL
cbknu(zr, fnu, kode, 2, cw, nw, tol, elim, alim)
22 CALL
crati(zr, fnu, n, y, tol)
27 cinu =
cmplx(1.0e0,0.0e0)
28 IF (kode.EQ.1) go
to 10
41 ascle = 1.0
e+3*r1mach(1)/tol
42 cscl =
cmplx(1.0e0,0.0e0)
43 IF (acw.GT.ascle) go
to 20
44 cscl =
cmplx(1.0e0/tol,0.0e0)
48 IF (acw.LT.ascle) go
to 30
49 cscl =
cmplx(tol,0.0e0)
60 rct =
cmplx(1.0e0/act,0.0e0)
subroutine cwrsk(ZR, FNU, KODE, N, Y, NZ, CW, TOL, ELIM, ALIM)
subroutine crati(Z, FNU, N, CY, TOL)
octave_value sin(void) const
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)
octave_value cos(void) const
subroutine cbknu(Z, FNU, KODE, N, Y, NZ, TOL, ELIM, ALIM)