Решение нелинейных уравнений. Методы бисекций, простой итерации, Ньютона.

function ex2
% Решить уравнение f(x)=0, где где f(x)= x^3 - cos(x) + 1 методом бисекции

% Введём функцию f(x)
f = inline('x.^3 - cos(x) + 1');
root1 = bisec(f, -0.6, -0.4)
root2 = bisec(f, -0.2, 0.2)

% Метод бисекции
function center = bisec(f, left, right)
% Классический метод бисекции с точностью 2 eps
while right - left > eps * 2
center = (right - left) / 2 + left;
if f(center) * f(left) > 0
left = center;
else
right = center;
end
end

>>
root1 = -0.4901

root2 = -1.0537e-008