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

aW1wb3J0IG51bXB5IGFzIG5wCgptQSA9IG5wLmFycmF5KFtbMSwgMiwgM10sWzQsIDUsIDZdLFs3LCA4LCA5XV0pCm1CID0gbnAuYXJyYXkoW1sxNCwgMzIsICA1MF0sWzMyLCA3NywgIDEyMl0sWzUwLCAxMjIsIDE5NF1dKQoKI3RyYW5zcG9zdGEKbVQgPSBtQS50cmFuc3Bvc2UoKQoKI211bHRpcGxpY2EKbVIgPSBucC5kb3QobUEsIG1UKQoKI2NvbXBhcmEgbWF0cml6ZXMKaWYgbnAuYXJyYXlfZXF1YWwobVIsIG1CKToKICAgIHByaW50ICJjb3JyZXRvXG4iCmVsc2U6CiAgICBwcmludCAiaW5jb3JyZXRvXG4iCgojYXBlbmFzIHBhcmEgbW9zdHJhciBvcyB2YWxvcmVzIGRhcyBtYXRyaXplcwpwcmludCAiTWF0cml6IEEiCnByaW50IG1BCiAgICAKcHJpbnQgIlxuTWF0cml6IEEnIgpwcmludCBtVAoKcHJpbnQgIlxuTWF0cml6QSAqIE1hdHJpekEnIgpwcmludCBtUgoKcHJpbnQgIlxuTWF0cml6IEIiCnByaW50IG1C