The bisection algorithm in Haskell.

bisection.hs

bisection f a b e = if a < b then bis a b else bis b a
    where err     = if e < 1e-15 then 1e-15 else e
          bis a b = let d = (b - a) / 2; m = (b + a) / 2 in
                    if d < err then
                        m
                    else if f a * f m < 0.0 then
                        bis a m
                    else
                        bis m b