************************************************************************
C math library
function ZEROIN - obtain a function zero within the given range
Output
Zeroin returns an estimate for the root with accuracy
4*EPSILON*abs(x) + tol
Algorithm
G.Forsythe, M.Malcolm, C.Moler, Computer methods for mathematical
computations.