import numpy as np
import matplotlib.pyplot as plt

A=np.array([[1,3,2],
       [4,8,2],
       [3,9,7],
       [22,11,17],
       [55,33,66]])


A=A-np.mean(A,0)
print(np.dot(A.T,A)/(A.shape[0]-1))
co=np.cov(A.T)
print(co)
#%%
[D,UI]=np.linalg.eigh(co)
[aa,bb,cc]=np.linalg.svd(co)
print (UI)
print(aa)