1 SUBROUTINE cbesk(Z, FNU, KODE, N, CY, NZ, IERR)
149 REAL AA, ALIM, ALN, ARG, AZ, DIG, ELIM, FN, FNU, FNUL, RL, R1M5,
150 * tol, ufl, xx, yy,
r1mach, bb
151 INTEGER IERR, K, KODE, K1, K2, MR, N, NN, NUF, NW, NZ, I1MACH
158 IF (yy.EQ.0.0e0 .AND. xx.EQ.0.0e0) ierr=1
159 IF (fnu.LT.0.0e0) ierr=1
160 IF (kode.LT.1 .OR. kode.GT.2) ierr=1
162 IF (ierr.NE.0)
RETURN
179 k = min0(iabs(k1),iabs(k2))
180 elim = 2.303e0*(float(k)*r1m5-3.0e0)
183 dig = amin1(aa,18.0e0)
185 alim = elim + amax1(-aa,-41.45e0)
186 fnul = 10.0e0 + 6.0e0*(dig-3.0e0)
187 rl = 1.2e0*dig + 3.0e0
189 fn = fnu + float(nn-1)
194 bb=float(i1mach(9))*0.5e0
196 IF(az.GT.aa) go
to 210
197 IF(fn.GT.aa) go
to 210
206 IF (az.LT.ufl) go
to 180
207 IF (fnu.GT.fnul) go
to 80
208 IF (fn.LE.1.0e0) go
to 60
209 IF (fn.GT.2.0e0) go
to 50
210 IF (az.GT.tol) go
to 60
213 IF (aln.GT.elim) go
to 180
216 CALL
cuoik(z, fnu, kode, 2, nn, cy, nuf, tol, elim, alim)
217 IF (nuf.LT.0) go
to 180
224 IF (nn.EQ.0) go
to 100
226 IF (xx.LT.0.0e0) go
to 70
230 CALL
cbknu(z, fnu, kode, nn, cy, nw, tol, elim, alim)
231 IF (nw.LT.0) go
to 200
239 IF (nz.NE.0) go
to 180
241 IF (yy.LT.0.0e0) mr = -1
242 CALL
cacon(z, fnu, kode, mr, nn, cy, nw, rl, fnul, tol, elim,
244 IF (nw.LT.0) go
to 200
252 IF (xx.GE.0.0e0) go
to 90
254 IF (yy.LT.0.0e0) mr = -1
256 CALL
cbunk(z, fnu, kode, mr, nn, cy, nw, tol, elim, alim)
257 IF (nw.LT.0) go
to 200
261 IF (xx.LT.0.0e0) go
to 180
268 IF(nw.EQ.(-1)) go
to 180
subroutine cacon(Z, FNU, KODE, MR, N, Y, NZ, RL, FNUL, TOL, ELIM, ALIM)
subroutine cbunk(Z, FNU, KODE, MR, N, Y, NZ, TOL, ELIM, ALIM)
subroutine cbesk(Z, FNU, KODE, N, CY, NZ, IERR)
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 cuoik(Z, FNU, KODE, IKFLG, N, Y, NUF, TOL, ELIM, ALIM)
subroutine cbknu(Z, FNU, KODE, N, Y, NZ, TOL, ELIM, ALIM)
octave_value sqrt(void) const