import numpy as np
from numpy.random import rand
from numpy.linalg import solve, inv
a = np.array([[1, 2, 3], [3, 4, 6.7], [5, 9.0, 5]])
print a.transpose()
print inv(a)
b =  np.array([3, 2, 1])
print solve(a, b)  # solve the equation ax = b

c = rand(3, 3)  # create a 3x3 random matrix
print c
print np.dot(a, c)  # matrix multiplication
