clear all
clc

syms x
disp('Algoritmo de Secante')
f=input(' Función a resolver / comilla : ');

P=input(' Intervalos [a b]: ');
a=P(1);
b=P(2);
tol=input(' Tolerancia (T): ');
max1=3;
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

RR=zeros(1,8);

for k=1:max1
    
    x=a; ya=eval(f);
    x=b; yb=eval(f);
    dif=(b-a);
    c=b-((dif)*(yb))/(yb-ya);
    x=c;
    err=abs((c-b)*100/c);
    RR(k,:)=[k a  b dif yb ya c err ];
    
    if  tol <= err
        a=b;
        b=c;
    else
        break
    end
end

  %format long   
disp('   "Iteracion"  "xa"    "xb"     "xb-xa"     "f(a)"   "f(b)"    "xi+1"       "e%"     ')
        disp(' ')
        disp(RR)

disp(' La aproximación de la raíz es:')
         disp(c)