Boa noite,
Estou tendo o seguinte problema:
O código desenvolvido em Matlab para calcular raízes da equação usando método da bisecção está bugando, talvez no processo de calculo (loops);
A ideia é a seguinte: O usuário entrar com a função que deseja, intervalos, tolerância desejada e a quantidade de iterações.
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);
aux
em nenhum lugar, o que possivelmente está resultando em um loop infinito.