1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
| #include <conio.h>
#include <stdio.h>
#include <windows.h>
#include <iostream>
#include <math.h>
using namespace std;
int main()
{
setlocale(LC_ALL,"Russian");
int itr, k;
const int m=3;
long double temp;
long double a[m][m]={ {5, 1, 2},
{1, 4, 1},
{1, 1, 3},
};
long double b[m] = {10, 7, 8};
long double c[m]={0,0,0};
itr=0;
long double c_old[m] = {0,0,0};
do {
itr++;
for (k=0; k<m; k++)
{
c[0]=(b[0]-a[0][1]*c[1]-a[0][2]*c[2])/a[0][0];
c[1]=(b[1]-a[1][0]*c[0]-a[1][2]*c[2])/a[1][1];
c[2]=(b[2]-a[2][0]*c[0]-a[2][1]*c[1])/a[2][2];
}
//x1:=(b1-a[1,2]*x2-a[1,3]*x3)/a[1,1];
//x2:=(b2-a[2,1]*x1-a[2,3]*x3)/a[2,2];
//x3:=(b3-a[3,1]*x1-a[3,2]*x2)/a[3,3];
for(int j=0; j <3 ; j++)
printf("x%d = %.10lf\n", j+1,c[j]);
cout<<endl;
temp=0;
for (k=0; k<m; k++) temp+=fabs(c_old[k]-c[k]);
for (k=0; k<m; k++) c_old[k]=c[k];
}
while (temp>0.001);
cout << "Кiлькiсть iтерацiй:" << itr << endl;
for(int j=0; j <3 ; j++) printf("x%d= %.10lf\n", j+1,c[j]);
system("PAUSE");
return 0;
} |