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