C++ program to find determinant of a 4×4 matrix.
C++ Code:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 |
#include <stdio.h> #include <stdlib.h> #include <time.h> #define d 4 //4x4 double det (double A[d][d], int N) { double c, r=1; for(int i = 0; i < N; i++) { for(int k = i+1; k < N; k++) { c = A[k][i] / A[i][i]; for(int j = i; j < N; j++) A[k][j]= A[k][j] - c*A[i][j]; } } for (int i = 0; i < N; i++) r *=A[i][i]; return r; } int main() { double M[d][d]; for(int i=0; i<d; i++) { for(int j=0; j<d; j++) { M[i][j] = rand() % 20; } } printf("Det(M) = %f\n",det(M,d)); return 0; } |
Output: