matrix1 = np.array([[0,1],[1,3]])

print(matrix1)

print(np.dot(matrix1,matrix1)) #matrix or dot product

print(matrix1*matrix1)

stuff = np.zeros(10,float)

t = np.arange(4)

stuff[3:7] = np.sqrt(t+1)

stuff[3:7]