Metodo De La Biseccion

disp('Metodo de la biseccion');
n=6;%numero de iteracciones

a=input('Ingrese el intervalo x1: ');
b=input('Ingrese el intervalo x2: ');
xna=0;
fprintf(' A\t C\t F(a)\t F(c)\t B\t F(b)\terror\');
for i=1:n
f(a)=1+a+exp(a);
f(b)=1+b+exp(b);
xn=(a+b)/2;
f(xn)=1+xn+exp(xn);
error=abs(((xn-xna)/xn)*100);
fprintf('\n %.3f %.3f %.3f %.3f %.3f %.3f %.3f\n',a,b,f(a),f(b),xn,f(xn),error);
if f(a)*f(xn)<=0
b=xn;
f(b)=f(xn);
else
a=xn;
f(a)=f(xn);
end

end
fprintf('\nEl resultado sera %.4f\n',xn);ezplot(Fx);%graficamos la funcion
grid on;


clc;
Fx=input('Ingrese la funcion: ','s');
a=input('Ingrese el intervalo x1 : ');
c=input('Ingrese el intervalo x2 : ');
e=input('Ingrese el error :');

x=a;
Fa=eval(Fx);
x=c;
Fc=eval(Fx);
fprintf('\n %6s %7s %8s %10s %8s %8s %8s \n ','A','C','F(a)','F(c)','B','F(b)','|c-a|');
while abs(c-a)>e
b=(a+c)/2;
x=b;Fb=eval(Fx);

fprintf('\n %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f \n',a,c,Fa,Fc,b,Fb,abs(c-a));
if Fa*Fb<=0
c=b;
Fc=Fb;
else
a=b;
Fa=Fb;
endend
fprintf('\nEl resultado sera %.4f\n',b);
ezplot(Fx);%graficamos la funcion
grid on;

Fx=input('Ingrese la funcion: ','s');
a=input('Ingrese el intervalo x1 : ');
c=input('Ingrese elintervalo x2 : ');
e=input('Ingrese el error : ');
x=a;
Fa=eval(Fx);
x=c;
Fc=eval(Fx);
fprintf('\n %6s %7s %8s %10s %8s %8s %8s \n ','A','C','F(a)','F(c)','B','F(b)','error');

clear all,clc;Fx=input('Ingrese la funcion: ','s');
a=input('Ingrese el intervalo x1 : ');
c=input('Ingrese el intervalo x2 : ');
e=input('Ingrese el error : ');
x=a;
Fa=eval(Fx);
x=c;
Fc=eval(Fx);
fprintf('A\t C\t F(a)\t F(c)\t B\t F(b)\t error\');

while (c-a)>e
b=(a+c)/2;
x=b;
Fb=eval(Fx);

fprintf('\n %.4f %.4f %.4f %.4f %.4f %.4f %.4f...