#include <stdio.h>
void main ()
{
	float **allocate_real_matrix(int, int, int, int);
	void free_real_matrix(float **, int, int, int);
	void gsssolerb(float **, int, float [], float []);
	int i,j;
	float **a,b[5],aux[12];

	a=allocate_real_matrix(1,4,1,4);
	for (i=1; i<=4; i++) {
		for (j=1; j<=4; j++) a[i][j]=1.0/(i+j-1);
		b[i]=a[i][3];
	}
	aux[0]=aux[2]=aux[6]=1.0e-14;
	aux[4]=8;
	gsssolerb(a,4,aux,b);
	printf("Solution: %e  %e  %e  %e\n",b[1],b[2],b[3],b[4]);
	printf("Sign(Det) =%3.0f\nNumber of eliminations =%3.0f\n"
			"Max(abs(a[i,j])) = %e\nUpper bound growth = %e\n"
			"1-norm of the inverse matrix = %e\n"
			"Upper bound rel. error in the calc. solution = %e",
			aux[1],aux[3],aux[5],aux[7],aux[9],aux[11]);
	free_real_matrix(a,1,4,1);
}