clear all
clc

disp('Punto Fijo:')
disp(' ')
Ea=100;
Ev=5;
while abs(Ea)>=0.01
    VV=sqrt((1.7*Ev+2.5)/0.9);
    Ea=((Ev - VV)/VV)*100;
    fprintf('Xi+1: %f Xi= %f Ea= %f',VV,Ev,Ea);
    Ev=VV;
    disp(' ')
end

disp(' ')
disp('Newton:')
disp(' ')
clear all
Ea=100;
Ev=5;
while abs(Ea)>=0.01
    VV=Ev-(-0.9*Ev^2+1.7*Ev+2.5)/(-1.8*Ev+1.7);
    Ea=((Ev - VV)/VV)*100;
    fprintf('Xi+1: %f Xi= %f Ea= %f',VV,Ev,Ea);
    Ev=VV;
    disp(' ')
end