#include #include // number of variables #define N 2 double A[N][N]; double x[N]; double b[N]; /* An simple implementation of Gauss elimination. It solves the equation Ax = b (A, b are given, A:N*N matrix, x:N-vector, b:N-vector). It transformes A into upper triangle matrix U a b c a b c a b c d e f -> 0 j k -> 0 j k g h i 0 l m 0 0 n An equation Ux = b' is easily can be solved. If it encounts A[i][i]==0, it doesn't work. */ void gauss_elimination_simple(){ int i, j, k; double tmp; // elimination for(i=0; i=0; i--){ x[i] = b[i]; for(j=i+1; j=0; i--){ x[i] = b[i]; for(j=i+1; j