How to make a generic implementation of the bisection method in MATLAB?

Good Night,

I am having the following problem:

• The code developed in Matlab to calculate roots of the equation using bisection method is bugging, perhaps in the calculation process (loops);

• The idea is as follows: The user enters with the function they want, ranges, desired tolerance and the amount of iterations.

clc
clear


%Recebe a Função desejada
disp('Insira a sua função');
f = input('=> ','s');

%Recebe os Intervalos desejados
disp('Insira os valor do Intervalo Xa');
Xa = input('=> '); %Recebe o valor do intervalo A
disp('Insira os valor do Intervalo Xb: ');
Xb = input('=> '); %Recebe o valor do intervalo B

%Recebe a Tolerancia desejada
disp('Insira a Tolerancia desejada');
tolerancia_desejada = input('=>'); %Recebe o valor do erro desejado

%Recebe a quantidade de Operações
disp('Insira o numero de interações');
iteracoes_desejada = input('=> '); %Recebe o valor da quantidade de iterações que o usuário deseja

%Processamento dos Dados
aux = 1;

if((subs(f,Xa))*(subs(f,Xb))>0)
    fprintf('Essa função não existe');
else
    if((subs(f(Xa))*(subs(f,Xb))<0))
        fprintf('Essa função tem raiz');
    end
end

while(aux<iteracoes_desejada)
    media = (Xa+Xb)/2;
    if(subs((f(Xa))*(subs(f(media))<0)))
        Xb = media;
    else
        Xa = media;
    end
    if abs(subs((f(media))<tolerancia_desejada))
        break
    end
end

%Exibindo Resultado da Operação
fprintf('A Raiz é: %f',media);