import numpy as np

def matrix_ops(matrix):
    array = np.array(matrix)
    return np.dot(array, array.T)

matrix = [[1,2], [3,4]]
result = matrix_ops(matrix)

print("Matrix Product: ", result)
print("Shape: ", result.shape)