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

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

using namespace std;

int main ()
{
double a[4][4],b[4],E,x1,x2,x3,y1,y2,y3,z1,z2,z3,k;
int m,L,i,j;
for (i=1; i<=3; i++)
{
for (j=1; j<=3;j++)
{
cout<<" vvedi a["<<i<<"]["<<j<<"] ";
cin>>a[i][j];
cout<<"\n";
}
}

cout<<"vvedi B1 "; cin>>b[1]; cout<<"\n";
cout<<"vvedi B2 "; cin>>b[2]; cout<<"\n";
cout<<"vvedi B3 "; cin>>b[3]; cout<<"\n";

cout<<"vvedi E"; cin>>E; cout<<"\n";
cout<<"vvedi m"; cin>>m; cout<<"\n";

x1=b[1];
x2=b[2];
x3=b[3];
L=0;

m1: y1=x1;
x1=(b[1]-a[1][2]*x2-a[1][3]*x3)/a[1][1];

y2=x2;
x2=(b[2]-a[2][1]*x1-a[2][3]*x3)/a[2][2];

y3=x3;
x3=(b[3]-a[3][1]*x1-a[3][2]*x2)/a[3][3];

z1=fabs(x1-y1);
z2=fabs(x2-y2);
z3=fabs(x3-y3);
cout<<"z1= "<<z1;cout<<"\n";
cout<<"z2= "<<z2;cout<<"\n";
cout<<"z3= "<<z3;cout<<"\n";
k=z1;
if (z2>k) k=z2;
if (z3>k) k=z3;
if (k<E) {cout<<"X1= "<<x1<<" X2= "<<x2<<" X3= "<<x3; goto m2;}

else {
L=L+1;
if (L<m) goto m1;else cout<<" SOS !!!!!!!!!!!";
}


m2:

getch ();
return 0;
}