While performing the naive gaussian elimination in C++, i can get a upper triangular matrix, but it doesn't work when the system of equation is more than 5*5 . any idea why this might be happening? i am using double as data type for the coefficients.
double factor =0;
factor = (doubleCoefficient[i][k])/(doubleCoefficient[k][k]);
doubleCoefficient[i][j] -= (doubleCoefficient[k][j]*factor);