	subroutine mcon(eta,infiny,smalno,base)

c magic numbers

	doubleprecision eta,infiny,smalno,base
	integer m,n,t

	base = 2.0d0
	t = 55
	m = 126
	n = -126
	eta = base**(1-t)
	infiny = base*(1.0d0-base**(-t))*base**(m-1)
	smalno = (base**(n+3))/base**3
	return
	end
on and multiplication

	doublecomplex q(nn)
	doubleprecision ms,mp,are,mre,e

	e = dcabs(q(1))*mre / (are+mre)
	do 10 i = 1,nn
		e = e*ms + dcabs(q(i))
10	continue
	errev = e*(are+mre) - mp*mre
	return
	end
the zero is stored and the
c polynomial