Bisection

Computes the root of a continuous function f in the interval [a,b] by bisectioning.

For this algorithm the function values of endpoints must differ in sign.

Usage: [c,fc,err] = bisection(f,a,b,delta)

Input arguments:
f function
a left endpoint
b right endpoint
delta error tolerance
Output arguments:
c solution
Optional Output arguments:
fc function value of solution
err error estimate