% LOGARITMO NATURAL DE MATRIZ
% -------------------------------------------------------
% Limpiamos el entorno
clc;
clear all;
close all;

% Informamos sobre inicio de ejecución
cadenainfo1 = "----------------------------------------------";
disp (cadenainfo1);
cadenainfo2 = "Programa de Cálculo de logaritmo natural de matriz I + A: método main";
disp (cadenainfo2);
disp(newline);

% *** CASO APARTADO 2 *****
cadenainfo1 = "Apartado 2: matriz A2";
disp(cadenainfo1);

A2 = 0.1 * eye(3);
tol = 1e-16;

[Asalida2,n2]=F(A2,tol);
disp("La salida de la función es:");
disp(Asalida2);
disp(n2);

% *** CASO APARTADO 3 *****
cadenainfo1 = "Apartado 3: matriz A3";
disp(cadenainfo1);
A3 = randi(10,3);
[Asalida3,n3]=F(A3,tol);
disp("La salida de la función es:");
disp(Asalida3);
disp(n3);

% *** CASO APARTADO 3b *****
cadenainfo1 = "Apartado 3b: matriz A3b";
disp(cadenainfo1);
A3b = rand(3);
A3b = A3b / ( 1.05 *  norm( A3b));
[Asalida3b,n3b]=F(A3b,tol);
disp("La salida de la función es:");
disp(Asalida3b);
disp(n3b);

disp(newline);
cadenainfo2 = "Fín de programa";
disp (cadenainfo2);
cadenainfo1 = "----------------------------------------------";
disp (cadenainfo1);