In [5]:

    

import random
from sympy import *
init_printing()


    

In [ ]:

    

Matrix??


    

In [6]:

    

A = Matrix([
[3,     3],
[2, S(3)/2]
])

A


    

    
Out[6]:





$$\left[\begin{matrix}3 & 3\\2 & \frac{3}{2}\end{matrix}\right]$$

In [8]:

    

A.inv()*Matrix([6,5])


    

    
Out[8]:





$$\left[\begin{matrix}4\\-2\end{matrix}\right]$$

In [7]:

    

A.det()


    

    
Out[7]:





-3/2

In [8]:

    

A.rref()


    

    
Out[8]:





(Matrix([
 [1, 0],
 [0, 1]]), [0, 1])

In [9]:

    

A.inv()


    

    
Out[9]:





Matrix([
[ -1,  2],
[4/3, -2]])

In [11]:

    

alpha = symbols('alpha')


    

In [12]:

    

A = Matrix([
[3,     3],
[2, alpha]
])


    

In [13]:

    

A.det()


    

    
Out[13]:





3*alpha - 6

In [14]:

    

# if alpha = 2 then det(A)=0  => A is not invertible
# otherwise A is invertible.


    

In [ ]:

    




    

In [ ]:

    




    

In [ ]:

    




    

In [ ]:

    




    

In [ ]: