>>> X = scipy.matrix('1 2 3; 4 5 6; 7 8 9')
>>> X*X.T
matrix([[ 14, 32, 50],
[ 32, 77, 122],
[ 50, 122, 194]])
Это дает вам больше, чем вы хотели, но это, несомненно, легко.
Или
>>> X = scipy.array([[1,2,3], [4,5,6], [7,8,9]])
>>> scipy.dot(X, X.T)
array([[ 14, 32, 50],
[ 32, 77, 122],
[ 50, 122, 194]])