1 SUBROUTINE cbuni(Z, FNU, KODE, N, Y, NZ, NUI, NLAST, FNUL, TOL,
14 COMPLEX CSCL, CSCR, CY, RZ, ST, S1, S2, Y, Z
15 REAL ALIM, AX, AY, DFNU, ELIM, FNU, FNUI, FNUL, GNU, TOL, XX, YY,
16 * ascle, bry,
str, sti, stm, r1mach
17 INTEGER I, IFLAG, IFORM, K, KODE, N, NL, NLAST, NUI, NW, NZ
25 IF (ay.GT.ax) iform = 2
26 IF (nui.EQ.0) go
to 60
28 dfnu = fnu + float(n-1)
30 IF (iform.EQ.2) go
to 10
35 CALL
cuni1(z, gnu, kode, 2, cy, nw, nlast, fnul, tol, elim, alim)
43 CALL
cuni2(z, gnu, kode, 2, cy, nw, nlast, fnul, tol, elim, alim)
51 bry(1) = 1.0
e+3*r1mach(1)/tol
57 cscl =
cmplx(ax,0.0e0)
58 IF (ay.GT.bry(1)) go
to 21
62 cscl =
cmplx(ax,0.0e0)
65 IF (ay.LT.bry(2)) go
to 25
69 cscl =
cmplx(ax,0.0e0)
72 cscr =
cmplx(ay,0.0e0)
75 rz =
cmplx(2.0e0,0.0e0)/z
78 s2 =
cmplx(dfnu+fnui,0.0e0)*rz*s2 + s1
81 IF (iflag.GE.3) go
to 30
88 IF (stm.LE.ascle) go
to 30
95 cscl =
cmplx(ax,0.0e0)
96 cscr =
cmplx(ay,0.0e0)
107 s2 =
cmplx(fnu+fnui,0.0e0)*rz*s2 + s1
113 IF (iflag.GE.3) go
to 40
119 IF (stm.LE.ascle) go
to 40
126 cscl =
cmplx(ax,0.0e0)
127 cscr =
cmplx(ay,0.0e0)
137 IF (iform.EQ.2) go
to 70
142 CALL
cuni1(z, fnu, kode, n, y, nw, nlast, fnul, tol, elim, alim)
150 CALL
cuni2(z, fnu, kode, n, y, nw, nlast, fnul, tol, elim, alim)
152 IF (nw.LT.0) go
to 50
subroutine cuni1(Z, FNU, KODE, N, Y, NZ, NLAST, FNUL, TOL, ELIM, ALIM)
subroutine cuni2(Z, FNU, KODE, N, Y, NZ, NLAST, FNUL, TOL, ELIM, ALIM)
subroutine cbuni(Z, FNU, KODE, N, Y, NZ, NUI, NLAST, FNUL, TOL, ELIM, ALIM)
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_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))<