#include<iostream>

using namespace std;
const int sizeX = 4;
const int sizeY = 2;
float matrix[sizeX][sizeY];
float X[sizeX][sizeY + 1];
float trans[sizeY + 1][sizeX];
float M[sizeX + 1][sizeY + 1];

void wypelnijMacierz()
{
	for (int i = 0; i < sizeX; i++)
	{
		for (int j = 0; i < sizeY; j++)
		{
			float a;
			cout << "Podaj element macierzy";
			cin >> a;
			matrix[i][j] = a;
		}

	}
}
void add() {

	for (int i = 0; i < sizeX; i++)
	{
		X[i][0] = 1;

		for (int i = 0; i < sizeX; i++)
		{
			for (int j = 0; j < sizeY; j++)
			{
				X[i][j + 1] = matrix[i][j];
			}

		}
	}
}
void transpose()
{
	for (int i = 0; i < sizeY + 1; i++)
	{
		for (int j = 0; j<sizeX; j++)

		{
			trans[i][j] = X[j][i];
		}
	}

}
void GetM()
{
	float tmp = 0;
	for (int x = 0; x <sizeY + 1; x++)
	{
		for (int y = 0; y < sizeY + 1; y++)
		{
			tmp = 0;
			for (int z = 0; z < sizeX; z++)
			{
				tmp += X[z][y] * trans[x][z];
			}	M[x][y] = tmp;
		}
	}


}
void PokazMacierz()
{
	for (int i = 0; i < sizeX; i++)
	{
		for (int j = 0; j < sizeY; j++)
		{
			cout << " " << matrix[i][j] << " ";
		}
	}

}
void PokazMacierzX()
{
	for (int i = 0; i < sizeX; i++)
	{
		for (int j = 0; j <sizeY + 1; j++)
		{
			cout << " " << X[i][j] << " ";

		}
	}
}
void pokazMacierzTrans()
{
	for (int i = 0; i < sizeY + 1; i++)
	{
		for (int j = 0; j < sizeX; j++)
		{
			cout << " " << trans[i][j] << " ";
		}
		cout << endl;
	}


}
void pokazMacierzM()
{
	for (int i = 0; i < sizeY + 1; i++)
	{
		for (int j = 0; j < sizeY + 1; j++)
		{
			cout << " " << M[i][j] << " ";
		}
		cout << endl;
	}

}
int main()
{
	wypelnijMacierz();
	PokazMacierz();
	add();
	PokazMacierzX();
	transpose();
	pokazMacierzTrans();
	GetM();
	pokazMacierzM();

	return 0;
	system("pause");


}