1 SUBROUTINE cacai(Z, FNU, KODE, MR, N, Y, NZ, RL, TOL, ELIM, ALIM)
18 COMPLEX CSGN, CSPN, C1, C2, Y, Z, ZN, CY
19 REAL ALIM, ARG, ASCLE, AZ, CPN, DFNU, ELIM, FMR, FNU, PI, RL,
20 * sgn, spn, tol, yy,
r1mach
21 INTEGER INU, IUF, KODE, MR, N, NN, NW, NZ
23 DATA pi / 3.14159265358979324e0 /
28 dfnu = fnu + float(n-1)
29 IF (az.LE.2.0e0) go
to 10
30 IF (az*az*0.25e0.GT.dfnu+1.0e0) go
to 20
35 CALL
cseri(zn, fnu, kode, nn, y, nw, tol, elim, alim)
38 IF (az.LT.rl) go
to 30
42 CALL
casyi(zn, fnu, kode, nn, y, nw, rl, tol, elim, alim)
49 CALL
cmlri(zn, fnu, kode, nn, y, nw, tol)
55 CALL
cbknu(zn, fnu, kode, 1, cy, nw, tol, elim, alim)
59 csgn =
cmplx(0.0e0,sgn)
60 IF (kode.EQ.1) go
to 50
64 csgn = csgn*
cmplx(cpn,spn)
71 arg = (fnu-float(inu))*sgn
75 IF (
mod(inu,2).EQ.1) cspn = -cspn
78 IF (kode.EQ.1) go
to 60
81 CALL
cs1s2(zn, c1, c2, nw, ascle, alim, iuf)
84 y(1) = cspn*c1 + csgn*c2
subroutine cseri(Z, FNU, KODE, N, Y, NZ, TOL, ELIM, ALIM)
octave_int< T > mod(const octave_int< T > &x, const octave_int< T > &y)
octave_value sin(void) const
subroutine cs1s2(ZR, S1, S2, NZ, ASCLE, ALIM, IUF)
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 casyi(Z, FNU, KODE, N, Y, NZ, RL, TOL, ELIM, ALIM)
subroutine cbknu(Z, FNU, KODE, N, Y, NZ, TOL, ELIM, ALIM)
subroutine cmlri(Z, FNU, KODE, N, Y, NZ, TOL)
subroutine cacai(Z, FNU, KODE, MR, N, Y, NZ, RL, TOL, ELIM, ALIM)