SONE :: Secant method

function x = secant_method(a, b, max_iter, max_error, my_function)

    for i = 1:max_iter
        x = b - feval(my_function, b) * ((b-a)/(feval(my_function, b)-feval(my_function, a)));

        if abs(feval(my_function, x)) < max_error
            return;
        end

        a = b;
        b = x;
    end

end