import numpy as np
def SistemaNormal(A,y):
    '''Retorna os parametros do sistema normal de um sistema sobredeterminado.'''
    k,n = np.shape(A)
    N = np.array([[sum(A.T[i]*A.T[j]) for j in range(n)] for i in range(n)])
    b = np.array([sum(A.T[i]*y) for i in range(n)])
    return N,b

#teste
A=np.array([[1,2,3,4,5],[2,3,4,5,6],[3,4,5,0,0]]).T
y=np.array([1,0,1,0,1])

N,u = SistemaNormal(A,y)

N

u