import numpy as np

mA = np.array([[1, 2, 3],[4, 5, 6],[7, 8, 9]])
mB = np.array([[14, 32,  50],[32, 77,  122],[50, 122, 194]])

#transposta
mT = mA.transpose()

#multiplica
mR = np.dot(mA, mT)

#compara matrizes
if np.array_equal(mR, mB):
    print "correto\n"
else:
    print "incorreto\n"

#apenas para mostrar os valores das matrizes
print "Matriz A"
print mA
    
print "\nMatriz A'"
print mT

print "\nMatrizA * MatrizA'"
print mR

print "\nMatriz B"
print mB