Численные методы. Лекции, задания, примеры / Metod_Gausa_4_var_kod

#include <iostream>
#include  <stdio.h>
#include  <math.h>
#include  <stdlib.h>
#include <conio.h>

using namespace std;

int main ()
{

	int i,j;
	double a[4][4],b[4],x[4];
	for (i=1; i<=3; i++)
	{
		cout<<"vvedi b["<<i<<"]  ";
		cin>>b[i]; cout<<"\n";
		for (j=1;j<=3;j++)
		{
			cout<<"vvedi a["<<i<<"]["<<j<<"] "<<"\n";
			cin>>a[i][j];
		}
	}
	
	for (i=2; i<=3; i++)
	{
		b[i]=b[i]-(a[i][1]*b[1])/a[1][1];
		
		for
			(j=2; j<=3; j++)
		{
			a[i][j]=a[i][j]-a[i][1]*a[1][j]/a[1][1];
		}
	}

	a[3][3]=a[3][3]-a[3][2]*a[2][3]/a[2][2];
	b[3]=b[3]-a[3][2]*b[2]/a[2][2];
	x[3]=b[3]/a[3][3];
	x[2]=(b[2]-a[2][3]*x[3])/a[2][2];
	x[1]=(b[1]-a[1][2]*x[2]-a[1][3]*x[3])/a[1][1];
	for (j=1; j<=3; j++)
	{
		cout<<"\n Znach X= "<<x[j];
	}

	getch ();
	return 0;
}