# Import sympy package:

import sympy as sym

# Defining symbols:

a, b = sym.symbols('a, b')

# Defining a symbolic matrix:

M1 = sym.Matrix([[a, 0], [0, b]])
print('M1 =', M1)

# Matrix multiplication:

M2 = M1*M1
print('M2 =', M2)

# Matrix inversion:

inv_M1 = M1**(-1)
print('inv_M1 =', inv_M1)

# Determinant:

det_M1 = M1.det()
print('det_M1 =', det_M1)

# Calculating num value of M1 using lambdify():

fun_M = sym.lambdify([a, b], M1, 'numpy')
val_M = fun_M(a=2, b=3)
print('val_M =', val_M)