#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Wed Nov 23 10:54:07 2022

@author: ojacques
"""

import numpy as np
import numpy.linalg as la
A=np.array([[0.59955556, 0.61444444],
       [0.61444444, 0.73555556]])
print('Dada uma matriz quadrada A:')
print(A)
print('\nExtraimos os seus autovalores D em diagonal:')
D,V = la.eigh(A)
D=np.diag(D)
print(D)
print('\nE os seus autovetores V correspondentes:')
print(V)
print('\n---------------------------')
print('VERIFICA-SE QUE  AV=VD')
print('---------------------------')
print('\nImprimindo AV:')
print(A@V)
print('\nImprimindo VD:')
print(V@D)
print('Porcentagem explicativa do primeiro autovetor:')
print(f'{(D[0,0]/np.sum(D))*100:.2f}%')
print('Porcentagem explicativa do segundo autovetor:')
print(f'{(D[1,1]/np.sum(D))*100:.2f}%')