1 SUBROUTINE cbesj(Z, FNU, KODE, N, CY, NZ, IERR)
145 COMPLEX CI, CSGN, CY, Z, ZN
146 REAL AA, ALIM, ARG, DIG, ELIM, FNU, FNUL, HPI, RL, R1, R1M5, R2,
147 * tol, yy,
r1mach, az, fn, bb, ascle, rtol, atol
148 INTEGER I, IERR, INU, INUH, IR, KODE, K1, K2, N, NL, NZ, I1MACH, K
150 DATA hpi /1.57079632679489662e0/
155 IF (fnu.LT.0.0e0) ierr=1
156 IF (kode.LT.1 .OR. kode.GT.2) ierr=1
158 IF (ierr.NE.0)
RETURN
174 k = min0(iabs(k1),iabs(k2))
175 elim = 2.303e0*(float(k)*r1m5-3.0e0)
178 dig = amin1(aa,18.0e0)
180 alim = elim + amax1(-aa,-41.45e0)
181 rl = 1.2e0*dig + 3.0e0
182 fnul = 10.0e0 + 6.0e0*(dig-3.0e0)
183 ci =
cmplx(0.0e0,1.0e0)
190 bb=float(i1mach(9))*0.5e0
193 IF(az.GT.aa) go
to 140
194 IF(fn.GT.aa) go
to 140
205 arg = (fnu-float(inu-ir))*hpi
209 IF (
mod(inuh,2).EQ.1) csgn = -csgn
214 IF (yy.GE.0.0e0) go
to 40
219 CALL
cbinu(zn, fnu, kode, n, cy, nz, rl, fnul, tol, elim, alim)
220 IF (nz.LT.0) go
to 120
231 IF (amax1(
abs(aa),
abs(bb)).GT.ascle) go
to 55
232 zn = zn*
cmplx(rtol,0.0e0)
236 cy(i) = zn*
cmplx(atol,0.0e0)
241 IF(nz.EQ.(-2)) go
to 130
octave_int< T > mod(const octave_int< T > &x, const octave_int< T > &y)
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)
subroutine cbesj(Z, FNU, KODE, N, CY, NZ, IERR)
octave_value cos(void) const
OCTAVE_EXPORT octave_value_list return the value of the option it must match the dimension of the state and the relative tolerance must also be a vector of the same length tem it must match the dimension of the state and the absolute tolerance must also be a vector of the same length The local error test applied at each integration step is xample roup abs(local error in x(i))<
octave_value sqrt(void) const
subroutine cbinu(Z, FNU, KODE, N, CY, NZ, RL, FNUL, TOL, ELIM, ALIM)