Текст программы

#include "stdafx.h"

#include <iostream>

#include <iomanip>

#include <math.h>

#include <stdio.h>

using namespace std;

long double eps=0.00001;

long double A[3][3]={{3.5, -1.7, 2.8}, {2.7, 3.3, 1.3}, {1.4, 2.5, -3.0}};

long double B[3]={1.7, 2.1, -1.3};

long double D[3], C[3][3];

long double n[3], p[3]={1.7, 2.1, -1.3};

bool check(long double X1[3], long double X2[3])

{

for (int i=0; i<3; i++)

if (abs(X1[i]-X2[i])>=eps) return true;

return false;

}

void Zeydel ()

{

n[0]=C[0][1]*p[1]+C[0][2]*p[2]+D[0];

n[1]=C[1][0]*n[0]+C[1][2]*p[2]+D[1];

n[2]=C[2][0]*n[0]+C[2][1]*n[1]+D[2];

}

int _tmain(int argc, _TCHAR* argv[])

{

setlocale(LC_ALL, "rus");

for (int i=0; i<3; i++)

{

D[i]=B[i]/A[i][i];

for (int j=0; j<3; j++)

{

C[i][j]=-A[i][j]/A[i][i];

}

C[i][i]=0;

}

cout << "Итерация" << "\tx1" << "\t\t\tx2" << "\t\t\tx3" << endl;

Zeydel();

int i=0;

while (check(n,p))

{

p[0]=n[0]; p[1]=n[1]; p[2]=n[2];

Zeydel();

i++;

cout<< i << "\t " << setprecision(15) << n[0] << "\t" << n[1] << "\t" <<n[2] << endl ;

};

cout<<"\n\nРешение: x1 = " << n[0] << endl << "\t x2 = " << n[1] << endl << "\t x3 = " << n[2] << "\n\n";

system ("pause");

return 0;

}

Вывод программы

Для сравнения, воспользуемся данными сайта Wolfram alpha:

Это означает, что мы нашли правильное решения для данной точности (Epsilon = 0.001)