use python

Use the rref function in sympy to find the inverse of the transformation and store it as P_inv. (Hint: you can use numpy. concatenate to join two 33 matrices into one 36 matrix.) Recall that, if we do Gauss-Jordan process on [PI3] where I3 is the 33 identity matrix, then at the end we get [I3P1] Note: Call P.evalf (2) to round 1.5000 to 1.5 in your final answer, if P is a sympy.matrix