In [17]:

    

import random
from sympy import *


    

In [ ]:

    

Matrix??


    

In [19]:

    

B = Matrix(3,3, lambda i,j: S(random.randrange(-1,4))/random.randrange(1,4) )
B


    

    
Out[19]:





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

In [20]:

    

B.det()


    

    
Out[20]:





25/27

In [8]:

    




    

In [ ]:

    




    

In [ ]:

    

beta = symbols('beta')


    

In [ ]:

    

A


    

In [ ]:

    

A.det()


    

In [ ]:

    

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


    

In [ ]:

    

S(random.randrange(1,4))/random.randrange(2,4)


    

In [ ]:

    

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


    

In [ ]:

    

A


    

In [ ]: