Metodo de la Falsa Posicion Matlab

Aqui les traigo el codigo del metodo de la falsa posicion hecho en matlab.

El codigo es el siguiente :

Fx=input('Ingrese la funcion: ','s');
a=input('Ingrese a : ');
c=input('Ingrese c : ');
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','B','C','F(a)','F(b)','F(c)','|c-a|');
while abs(c-a)>e
    b=(c*Fa-a*Fc)/(Fa-Fc);
    x=b;
    Fb=eval(Fx);
    fprintf('\n %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f \n',a,b,c,Fa,Fb,Fc,abs(c-a));
    if abs(Fc)<e
        break;
    else    
        if Fa*Fb<=0
            c=b;
            Fc=Fb;
        else
            a=b;
            Fa=Fb;
        end
    end
end
fprintf('\nEl resultado sera %.4f',b);
ezplot(Fx);%graficamos la funcion
grid on;
    

Espero les haya servido, comenten sobre cualquier duda ...