#include <stdio.h>
#include <stdlib.h>

void read (float **ptr, int n);
void read2 (float **ptr, int n);

int main (void)
{
float co1;
float co2;
float co3;
//float **ptr=NULL;
float temp;
int n = 3;
int column, column1, row1, column11, row11;
int row;
int i, ii;


float **ptr = (float **)malloc(n * sizeof(float *));


for (i=0; i<n;i++)
{

	ptr[i] = (float *) malloc ((n+1) * sizeof(float)); //array of pointers

}


printf ("Please input the matrix A of coefficients column by column.");

read(ptr, n);




for (column1=0; column1<n; column1++)
{
	for (row1=0; row1<n; row1++)
	{

	printf ("%g", ptr[row1][column1]);

	}
}



printf ("\nPlease input solutions column:\n");

read2(ptr,n);

for (column11=(n+1); column11<(n+2); column11++)
{
	for (row11=0; row11<n; row11++)
	{

	printf ("%g", ptr[row11][column11]);

	}
}


for (ii=0; ii<n; ii++)
{
	free(ptr[ii]);
}
free(ptr);



return (0);


}


void read (float **ptr, int n)
{
	int column, row;

	for (column=0; column<n; column++)
	{
		for (row=0; row<n; row++)
		{
			scanf ("%g", &ptr[row][column]);
		}
	}
}

void read2 (float **ptr, int n)
{
	int column, row;
	n++;
	for (column=n; column<(n+1); column++)
	{
		for (row=0; row<(n-1); row++)
		{
			scanf ("%g", &ptr[row][column]);
		}
	}
}