%-------------------------------------GUÍA 1: EJERCICIO 5------------------------------------%

clear all
close all

%-------------------------------INICIALIZACIÓN DE LAS VARIABLES------------------------------%

N=17.35; 		% Número al que se quiere calcular la raíz
n=3; 			% Grado de la raíz
error_relativo=0.01; 	% Cota del error relativo (Criterio de corte)
error_actual=1; 	% Inicialización de la variable con un valor distinto de 0
indice=1; 		% Inicialización de la variable con un valor distinto de 0

% Se elige un valor semilla para empezar a calcular (es recomendable elegir un valor cercano a
% la solución para comenzar) y se guarda en la posición 1 del vector.

sol(indice)=3;

%-------------------------RESOLUCIÓN UTILIZANDO LA INSTRUCCION WHILE-------------------------%

while error_actual>=error_relativo    
    sol(indice+1)=sol(indice)+((N-(sol(indice))^n)/(n*(sol(indice))^(n-1)));
    error_actual=abs((sol(indice+1)-sol(indice))/sol(indice));
    indice=indice+1; % Se aumenta el índice en una unidad
end

% Se muestran los resultados por consola
disp('UTILIZANDO SENTENCIA WHILE')
disp('Numero a calcular la raiz:'); disp(N)
disp('Grado de la raiz:'); disp(n)
disp('Resultado:'); disp(sol(indice))
disp('Error:'); disp(error_actual)

%--------------------------------------------------------------------------------------------%

clear all  % Se eliminan las variables creadas. Se van a tener que inicializar nuevamente

%-------------------------------INICIALIZACIÓN DE LAS VARIABLES------------------------------%

N=17.35;
n=3;
error_relativo=0.01;
error_actual=1;
indice=1; 
sol(indice)=3;     

%---------------------RESOLUCIÓN UTILIZANDO INSTRUCCIONES FOR, IF Y BREAK--------------------%

for indice=1:500 	% Se elige a priori un número grande porque no se sabe cuantas 
			% iteraciones se necesitan
    sol(indice+1)=sol(indice)+((N-(sol(indice))^n)/(n*(sol(indice))^(n-1)));
    error_actual=abs((sol(indice+1)-sol(indice))/sol(indice));
    if error_actual<error_relativo
        break;		% Corta el ciclo en el momento que se cumple la condición de la cota
    end
end

disp('UTILIZANDO SENTENCIA FOR, IF Y BREAK')
disp('Numero a calcular la raiz:'); disp(N)
disp('Grado de la raiz:'); disp(n)
disp('Resultado:'); disp(sol(indice))
disp('Error:'); disp(error_actual)

%--------------------------------------------------------------------------------------------%