Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- def my_det(X):
- '''
- Parameters
- ----------
- X : array_like
- Returns
- -------
- det : float
- Determinant of `a`.
- '''
- if(len(X)!=len(X[0])):
- return 0
- change = 0
- det = 1
- for i in range(len(X)):
- j=i
- while X[i][j] == 0:
- j+=1
- if j>=len(X):
- return 0
- if i!=j:
- for k in range(len(X)):
- X[i][k], X[j][k]=X[j][k], X[i][k]
- change+=1
- j+=1
- while(j!=len(X)):
- for k in range(len(X)):
- X[j][k]-=X[i][k]*(X[j][i]/X[i][i])
- j+=1
- for i in range(len(X)):
- det *= X[i][i]
- if(change % 2 == 1):
- det*=-1
- return det
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement