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
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
| #include <iostream>
#include <cstdlib>
#include <math.h>
#include <fstream>
using namespace std;
double relax(double **A,double *b,int n,double ep)
{int i,j,t,nom; double *R0,*R1,*X,*C,**D,S,Max,dx; bool Pr;
R0=new double [n];
R1=new double [n];
X=new double [n];
C=new double [n];
D=new double *[n];
for(i=0;i<n;i++)
D[i]=new double [n];
for(i=0;i<n;i++)
{for(j=0;j<n;j++)
D[i][j]=-A[i][j]/A[i][i];
C[i]=b[i]/A[i][i];}
for(i=0;i<n;i++)
X[i]=0;
for(i=0;i<n;i++)
{S=0;
for(j=0;j<n;j++)
if(i=j) S=S+D[i][j]*X[j];
R1[i]=C[i]-X[i]+S;};
for(;Pr=0;)
{for(i=0;i<n;i++)
R0[i]=R1[i];
Max=R0[0]; nom=0;
for(i=0;i<n;i++)
if(R0[i]>Max) {Max=R0[i];nom=i;};
dx=Max;
for(i=0;i<n;i++)
{if(i=nom) R1[i]=0;
else R1[i]=R0[i]+D[i][nom]*dx;}
X[nom]=X[nom]+dx;
Pr=0;
for(i=0;i<n;i++)
if (R1[i]>ep) Pr=1;}
return *X;};
///////////////////////////////////////////////////////////////////////////
int main()
{setlocale(LC_ALL, "rus");
int i,j,k,n; ifstream g;
double **A,*b,*x,ep=0.001;
g.open("D:\\ab.txt");
g>>n;
x=new double [n];
b=new double [n];
A=new double *[n];
for(i=0;i<n;i++)
A[i]=new double[n];
for(i=0;i<n;i++)
{for(j=0;j<n;j++)
g>>A[i][j];}
for(i=0;i<n;i++)
g>>b[i];
b[3]=30-90+50; b[4]=90-10-20; b[5]=-50+20-40;
g.close();
cout<<"Размерность матрицы:"<<n<<"x"<<n<<endl;
cout<<"Исходная матрица токов:";
for(i=0;i<n;i++)
{cout<<endl;
for(j=0;j<n;j++)
cout<<A[i][j]<<"\t";}
cout<<endl<<endl;
cout<<"Массив свободных членов:"<<endl;
for(i=0;i<n;i++)
cout<<b[i]<<"\t";
cout<<endl<<endl;
for(i=0;i<n;i++)
x[i]=relax(A,b,n,ep);
cout<<"Результат"<<endl;
for(i=0;i<n;i++)
cout<<x[i]<<"\t";
cout<<endl;
system("PAUSE");
return 0;} |
(
