Метод половинного деления

clc;
clear;
% Основные условия
f=@(x) cos(2*x)+x-5; %это искомая функция
a=5;
b=6;
e=0.0001;
 
% Осноной цикл 
while abs(b-a)>e
    c=(b+a)/2;
    if sign(f(c))==sign(f(a))
        a=c;
    else
        b=c;
    end
end
disp(['Ответ x=' num2str(c,8)]);
disp('для проверки точный ответ =');
syms x;
disp(solve(f(x)));

if sign( cos(2*c)+c-5 )==sign( cos(2*a)+a-5 )

disp(solve('cos(2*x)+x-5'));

f1 = tan(x.*0.44 + 0.3)-x.^2;
f2 = x.^3 - x.^2*0.1 + x.*0.4 + 1.2;
a1 = -2;
b1 = 1;
a2 = -3;
b2 = 3;
e = 0.001; %точность (1 тысячная)
to4ka_a_a = a1; %точка А для отрезка (а; b) уравнения а)
to4ka_b_a = b1; %точка В для отрезка (a; b) уравнения a)
to4ka_a_b = a2; %точка А для отрезка (а; b) уравнения б)
to4ka_b_b = b2; %точка В для отрезка (а; b) уравнения б)
 
while(abs(b1-a1)>e) %пока модуль (В-А) больше точности "е"
    c = (b1+a1)/2; %делим отрезок точкой С пополам
    if sign(tan(c.*0.44 + 0.3)-c.^2)==sign(a1.^3 - a1.^2*0.1 + a1.*0.4 + 1.2) %функция sign() определяет знак функции
        a1 = c;
    else b1 = c;
    end
end
disp(['Ответ к уравнению А: приближенный корень х = ' num2str(c,8)]); %функция num2str
%преобразует число в строку
disp('Проверка точного ответа: ');
syms x; %команда syms создает символьную переменную х
disp(solve(f1(x)));
 
while(abs(b2-a2)>e)
    c = (b2+a2)/2;
    if sign(tan(c.*0.44 + 0.3)-c.^2)==sign(a2.^3 - a2.^2*0.1 + a2.*0.4 + 1.2) 
        a2 = c;
    else b2 = c;
    end
end
disp(['Ответ к уравнению B: приближенный корень х = ' num2str(c,8)]); %функция num2str
%преобразует число в строку
disp('Проверка точного ответа: ');
syms x; %команда syms создает символьную переменную х
disp(solve(f1(x)));

f1 =@(x) tan(x.*0.44 + 0.3)-x.^2;
f2 =@(x) x.^3 - x.^2*0.1 + x.*0.4 + 1.2;