Main.m

%% Universidade Estadual de Campinas
%% IM325 - Metodos Matematicos para Sistemas Mecanicos
%% Autor : Wanderson Vinicius de Oliveira Monteiro

%% Limpeza de memoria 
clc; close all;clear
n = 20;  % Ordem da matriz de Hilbert
%% Exata
% H = sym(zeros(n));  % cria matriz simbólica 
% for i = 1:n
%     for j = 1:n
%         H(i,j) = sym(1) / sym(i + j - 1);
%     end
% end

%% Ponto flutuante
H = zeros(n);  % Cria uma matriz de zeros de ordem n

for i = 1:n
    for j = 1:n
        H(i,j) = 1 / (i + j - 1);  % Formula de Hilbert
        %% Referencia: https://en.wikipedia.org/wiki/Hilbert_matrix
    end
end

% A_inv_Jordan   = Gauss_Jordan_Inversa(H,"flutuante")

[A_inv_eqn,H1]  = Inversa_Hilbert(n,"flutuante")