bisect=function(fun, x0, x1 , eps)
{
        iterno = 0
        repeat {
                iterno = iterno + 1
                x2 = (x0+x1)/2
                if (fun(x2)*fun(x0)<0)
                       x1=x2
                else  
                       x0=x2
                if(abs(x1-x0) < eps||abs(fun(x2))<.1e-10)
                        break
                cat("****** Iter. No: ", iterno, "    Current Iterate = ",
                    x2, fill = T)
        }
        return(x2)
}