#include<stdio.h>
#include<stdlib.h>
#include<math.h>
float Gauss(int n, float A[][4], float B[], float X[])
{
	    int i,j,k;
	    float m,d;
	    for(i=1;i<=n-1;i++)
	    {

	    for(k=i+1;k<=n;k++)
	    {
		m=A[k][i]/A[i][i];
		for(j=1;j<=n;j++)
		A[k][j]=A[k][j]-m*A[i][j];
		B[k]=B[k]-m*B[i];
		}

}
       d=1;
       if(A[n][n]==0)
       return 0;
       for(i=1;i<=n;i++)
       d=d*A[i][i];
       for(i=n;i>=1;i--)
       {
       	   X[i]=B[i];
       	   for(j=n,j>=i+1,j--)
       	   X[i]=X[i]-A[i][j]*X[j];
       	   X[i]/A[i][i];
       }

}
return d;
}
int main()
{
	float A[4][4]={{0,0,0,0},{0,6,-1,-1},{0,-1,6,-1},{0,-1,-1,6}};
    float B[4]={0,11.33,32,42},X[4];
	int n=3,i;
	printf("Determiantul sistemului este %f\n",Gauss(n,A,B,X));
	printf("Solutia este ");
	for(i=1;i<=n;i++)
	printf("%f\n",X[i]);
	system("pause");
	return 0;

}