clear 
clc

syms x

disp('Algoritmo de Regula Falsin')
f=input(' Función a resolver / comilla: ');
P=input(' Intervalos  [xI xD]: ');
a=P(1);
b=P(2);
tol=input(' Tolerancia: ');
max1=1000;
RR=[1 P 0];

x=a; ya=eval(f);
x=b; yb=eval(f);

if ya*yb>0
      disp('Note: f(a)*f(b)>0')
     
end

max1=1000;
RR=zeros(1,5);

for k=1:max1
    
    x=a; ya=eval(f);
    x=b; yb=eval(f);
    dx=yb*(b-a)/(yb-ya);
    c=b-dx;
    x=c;
    yc=eval(f);
    ac=c-a;
    RR(k,:)=[k a  b  c yc ];
    
    if yc==0
        break
    elseif yb*yc>0
        b=c;
        yb=yc;
    else
        a=c;
        ya=yc;
    end
    dx=min(abs(dx),ac);
    if dx<tol
        break
    end
end


        disp(' ')
        disp('"Resultados Finales" ')
        disp('   "Iteracion"  "xI"    "xD"        "xM"       "f(x)"     ')
        disp(' ')
        
        
        disp(RR)
disp(' La aproximación de la raíz es:')
         disp(c)