first lab / division (1)

function [a, b, x0] = divison(a, b, func, n)
% Проверка значений функции на границах интервала
f_a = func(a);
f_b = func(b);

% Убедимся, что функция имеет разные знаки на концах интервала
if f_a * f_b > 0
error('Функция должна иметь разные знаки на концах интервала.');
end

for i = 1:n
x0 = (a + b) / 2; % Находим середину отрезка
f_x0 = func(x0); % Вычисляем значение функции в x0

% Проверяем, какой отрезок сужать
if f_x0 * f_a < 0
b = x0; % Сужаем отрезок
f_b = f_x0; % Обновляем значение функции в b
else
a = x0; % Сужаем отрезок
f_a = f_x0; % Обновляем значение функции в a
end

% Проверка на точность
if abs(b - a) < 1e-6
fprintf('Корень найден: x0 = %.6f\n', x0);
return; % Завершение функции
end
end
error_bisection =abs(func(a))
disp(['Невязка метода деления отрезков = ', num2str(error_bisection)])

fprintf('Приблизительное значение корня: x0 = %.6f\n', x0);
end