Advertisement
Jokker115

Matrix determinant

Apr 20th, 2017
82
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.85 KB | None | 0 0
  1. import random
  2. #import numpy
  3.  
  4. def fill_matrix (size):
  5.     for i in xrange(size):
  6.         for j in xrange(size):
  7.             matrix[i][j] = random.randint(1, size)
  8.     return matrix
  9.  
  10. def minor_matrix (matrix, size, i, j):
  11.     di = 0
  12.     for ki in xrange(size-1):
  13.         if (ki == i):
  14.             di = 1
  15.         dj = 0
  16.         for kj in xrange(size-1):
  17.             if (kj == j):
  18.                 dj = 1
  19.             tmp[ki][kj] = matrix[ki + di][kj + dj]
  20.     return tmp
  21.                
  22. def det_matrix (size):
  23.     if (size > 1):
  24.         for i in xrange(size):
  25.             for j in xrange(size):
  26.                 M = minor_matrix(matrix, 100, i, j)
  27.                 A = (-1) ** (i + j) * M
  28.                 det = matrix[i][j] * A
  29.         return det_matrix(size) + det_matrix(size - 1)
  30.     elif (size == 1):
  31.         return det
  32.  
  33.  
  34. matrix = [[0] * 100 for i in xrange(100)]
  35. matrix = fill_matrix(100)
  36. #matrix2 = numpy.array(matrix)
  37. #print matrix2
  38. det = det_matrix (100)
  39. print det
  40. input ()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement