clear all
clc
 
disp('Secante:')
disp(' ')
Ea=100;
xi_1=.5;                  %Este es un valor a editar
xi=-.1;                   %Este es un valor a editar
aux=1;
while abs(Ea)>=0.01
    VV=xi-(-11-22*(xi)+17*xi^2-2.5*xi^3)*(xi_1-xi)/((-11-22*(xi_1)+17*(xi_1)^2-2.5*(xi_1)^3)-(-11-22*(xi)+17*xi^2-2.5*xi^3));
    Ea=((xi - VV)/VV)*100;     %Nótese que se asume como el xi-1 anterior para error relativo el valor xi=-0.1 presente en las constantes
    fprintf('Xi+1: %f Xi= %f Ea= %f',VV,xi,Ea);
    xi_1=xi;
    vec1(1,aux)=Ea;
    vec2(1,aux)=aux;
    xi=VV;
    aux=aux+1;
    disp(' ')
end
plot(vec2,vec1);
 
syms VV %Sentencia para graficar la función
    M=-11-22*(VV)+17*VV^2-2.5*VV^3;    
    figure(2);
    ezplot(M,[-1 6]);