Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- swit=10
- j=0
- det=0
- import numpy
- matrix=numpy.random.rand(swit,swit)
- print matrix
- a= matrix
- print numpy.linalg.det(matrix)
- def get_minor(a, Iind, Jind):
- b = []
- for i in xrange(len(a)):
- if i == Iind:
- continue
- b.append([])
- for j in xrange(len(a)):
- if j == Jind:
- continue
- b[-1].append(a[i][j])
- return b
- def get_det(a):
- if len(a) == 2:
- return a[0][0] * a[1][1] - a[0][1] * a[1][0]
- res = 0.0
- mn = 1.0
- for i in xrange(len(a)):
- val = a[0][i]
- minor = get_minor(a, 0, i)
- res += mn * val * get_det(minor)
- mn = -mn
- return res
- print get_det(matrix)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement