1 SUBROUTINE zbunk(ZR, ZI, FNU, KODE, MR, N, YR, YI, NZ, TOL, ELIM,
13 DOUBLE PRECISION ALIM, AX, AY, ELIM, FNU, TOL, YI, YR, ZI, ZR
14 INTEGER KODE, MR, N, NZ
17 ax = dabs(zr)*1.7321d0
19 IF (ay.GT.ax) go
to 10
24 CALL
zunk1(zr, zi, fnu, kode, mr, n, yr, yi, nz, tol, elim, alim)
32 CALL
zunk2(zr, zi, fnu, kode, mr, n, yr, yi, nz, tol, elim, alim)
subroutine zbunk(ZR, ZI, FNU, KODE, MR, N, YR, YI, NZ, 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)
subroutine zunk1(ZR, ZI, FNU, KODE, MR, N, YR, YI, NZ, TOL, ELIM, ALIM)
subroutine zunk2(ZR, ZI, FNU, KODE, MR, N, YR, YI, NZ, TOL, ELIM, ALIM)