gistfile1.c

void gaussElim(int n, int m, double** A,double* x) {
	int k, i, j,s=0;
	//double x;
	for (k = 0; k <n-1; k++) {

		for (i = k + 1; i < n; i++) {


			for (j = k+1; j < n; j++) {



				A[i][j] = A[i][j] - (A[i][k]/A[k][k] * (A[k][j]));

				//printf("%lf\n",A[i][j]);


				A[i][n] = A[i][n]-((A[i][k] / A[k][k]) * (A[k][n]));

			}
		}
		A[k+1][k]=0;
	}





	 for(i=n-1;i!=0;i--){
		 s=0;


	 for(j=k+1;j < n;j++){

	 s+=A[i][j]*x[j];


	 }
	 x[i]=(A[i][n]-s)/A[i][i];

	 }


}