Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- double deter(double A[6][6], int n)
- {
- double b[6][6],d = 0;
- int k = 1;
- if(n==1) d = A[1][1];
- else if(n==2) d = A[1][1]*A[2][2]-A[1][2]*A[2][1];
- else if(n==3) d = A[1][1]*A[2][2]*A[3][3]+A[2][1]*A[3][2]*A[1][3]+A[1][2]*A[2][3]*A[3][1]-A[3][1]*A[2][2]*A[1][3]-A[2][1]*A[1][2]*A[3][3]-A[1][1]*A[2][3]*A[3][2];
- else {
- for(int g = 1; g<=n;g++){
- for(int i = 2; i<=n; i++){
- for(int j = 1; j<=n-1;j++){
- if(k==j) continue;
- else{
- b[i][j] = A[i][j];
- }
- }
- }
- d += pow(-1,k-1)*A[1][g]*deter(b,n-1);
- k++;
- }
- }
- return d;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement