Численные методы. Лекции, задания, примеры / 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;
}