// Win32Lab.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"
#include <iostream>
#include <vector>
#include <cmath>
using namespace std;
int _tmain(int argc, _TCHAR* argv[])
{
	double a[4][5] = {  {  0.0309, -0.0179, -0.0353, 0.048,   0.082 },   		
						{ -0.0179, 0.0114, 0.0203, -0.0298, -0.0481 },
						{ -0.0353, 0.0203, 0.0426, -0.0546, -0.1018 },
						{  0.048, -0.0298, -0.0564, 0.0849,   0.1452 } };

	int size = 4;
	//cin >> size;

	vector <vector <long double> > matrix;

	// ??????? ????? ????? ?????? (size) x (size + 1),
	// c ?????? ??????? ????????? ??????    
	matrix.resize(size);
	for (int i = 0; i < size; i++)
	{
		matrix[i].resize(size + 1);

		for (int j = 0; j < size + 1; j++)
		{
			matrix[i][j] = a[i][j];
			//cin >> matrix[i][j];
		}
	}

	// ????????? ??????????? ???????? ???????
	long double eps;
	// cin >> eps;
	eps = 0.001;
	// ?????? ?????? ???????? ??????????? ?? ?????????? ????????,
	// ?????? ???????? ????? ????? ????? ? ???????, ?.?. size,
	// ?????? ???????? ?????? ?????????? ????????? ??? ??????
	vector <long double> previousVariableValues(size, 0.0);

	// ????? ????????? ???????????? ??????? ?? ??? ???, 
	// ???? ?? ????? ?????????? ??????????? ????????    
	while (true)
	{
		// ?????? ?????? ???????? ??????????? ?? ??????? ????       
		vector <long double> currentVariableValues(size);

		// ????????? ???????? ??????????? ?? ??????? ????????
		// ? ???????????? ? ?????????????? ?????????
		for (int i = 0; i < size; i++)
		{
			// ?????????????? i-?? ??????????? ????????? 
			// ?????????? ????? i-?? ?????? ???????
			currentVariableValues[i] = matrix[i][size];

			// ???????? ????? ?? ???? ???????? ?? i-?? ???????????
			for (int j = 0; j < size; j++)
			{
				if (i != j)
				{
					currentVariableValues[i] -= matrix[i][j] * previousVariableValues[j];
				}
			}

			// ????? ?? ??????????? ??? i-?? ???????????
			currentVariableValues[i] /= matrix[i][i];
		}

		// ????????? ??????? ??????????? ???????????? ?????????? ????????
		long double error = 0.0;

		for (int i = 0; i < size; i++)
		{
			error += fabs(currentVariableValues[i] - previousVariableValues[i]);
		}

		// ???? ??????????? ???????? ??????????, ?? ????????? ???????
		if (error < eps)
		{
			break;
		}

		// ????????? ? ????????? ????????, ??? 
		// ??? ??????? ???????? ??????????? 
		// ?????????? ?????????? ?? ?????????? ????????
		previousVariableValues = currentVariableValues;
	}

	// ??????? ????????? ???????? ??????????? ? 8 ??????? ????????
	for (int i = 0; i < size; i++)
	{
		printf("%.8llf ", previousVariableValues[i]);
	}

	return 0;
}