clear all ;
clc ;

format long
syms x
disp('Algoritmo de Punto Fijo')
f=input(' Función a resolver f(x) / comilla: ');
g=input(' Función despejada g(x) / comilla : ');
P=input(' Primer Punto xi: ');
tol=input(' Tolerancia: ');
max1=100;


RR=[0 P 0 0  0 ];

for k=2:max1
    x=P(k-1);
    P(k)=eval(g);
    P1(k)=eval(f);
    err=abs(P(k)-P(k-1));
%     err=abs((P(k)-P(k-1))/abs(P(k)));
    relerr=err/(abs(P(k))+eps);
    p=P(k-1);
    p1=P(k);
    p2=abs(P1(k));
    RR(k-1,:)=[k-2   p  p1  p2 err];
    
  
        
    if k==max1
        disp(' ')
        disp('"Cambiar el punto de inicio o la ecuación despejada "')
        disp(' ')
    end
     
%     if (err<tol) | (relerr<tol)
        
    if  (relerr<tol)
        disp(' ')
        disp('"Resultados Finales" ')
        disp('             i              x_i                g(x_1)          abs(f(x_i))      error  x_i+1 - x_i ')
        disp(RR)
        disp(' La aproximación de la raíz es:')
         disp(p1)
        break
         
    end

end