#include <iostream>

using namespace std;


int main() {
	double **A = new double *[3];


	double x = 0;
	double y = 0;
	double z = 0;
	
	float tmp_x, tmp_y, tmp_z;

	for(int i = 0; i < 3; i++){
		A[i] = new double [3];
	}
		A[0][0] = 5;
		A[0][1] = -2;
		A[0][2] = 3;
		A[1][0] = 2;
		A[1][1] = 4;
		A[1][2] = 2;
		A[2][0] = 2;
		A[2][1] = -1;
		A[2][2] = -4;

		double *B = new double[3];
		B[0] = 10;
		B[1] = 0;
		B[2] = 0;
	

	for(int j = 0; j < 10; j++){
		tmp_x = (B[0] - A[0][1] * y - A[0][2] * z)/A[0][0]; 
		tmp_y = (B[1] - A[1][0] * x - A[1][2] * z)/A[1][1];
		tmp_z = (B[2] - A[2][0] * x - A[2][1] * y)/A[2][2];

		x = tmp_x;
		y = tmp_y;
		z = tmp_z;

		cout << "X = " << x << " Y = " << y << " Z: " << z << endl;
	}



	cout << A[0][0] * x + A[0][1] * y + A[0][2] * z << endl;
	cout << A[1][0] * x + A[1][1] * y + A[1][2] * z << endl;
	cout << A[2][0] * x + A[2][1] * y + A[2][2] * z << endl;

	system("PAUSE");
}