Программистам. Нужен код программы (простой, на С/С++) , решающий СЛАУ методом Зейделя. Срочно.

Нужен код программы (простой, на С/С++) , решающий СЛАУ методом Зейделя. Срочно.

Гаусс-Зейдель.
Дана матрица 4х4.
Что-то вроде этого:
#include
#include
#include
#include

double A[10][10],T[10][10],X[10],Y[10],C[10][10],D[10];
int n;
double E;

int vvodsort()
{
int i,j,k,max;
cout<<"Vvedi pogreshnosti ";
cin>>E;
cout<<"\nVvedi n ";
cin>>n;
cout<<"\nVvedi matritu koeffitientov A \n";
for(i=0;i>A[j];
T[j]=A[j];
}
cout<<"\nVvedi matritu svob koeffitientov B\n";
for(i=0;i>A[j];

for(i=0;i)>abs(A))max=j;
}

for(k=0;k[k];
int y;
y=A[n];
A[n]=A[max][n];
A[max][n]=y;
for(k=0;k[k]=A[max][k];
A[max][k]=Y[k];
}
}
}

int proverka()
{
int i,j,k,z=0;double sum;

for(i=0;i[j]);
if(abs(A)))){z=1;break;}
}
if(z==0)return 1; else return 0;
}

int zegeli()
{
int i,j;
for(i=0;i=X=A[n]/A;
D=A[n]/A;
}

for(i=0;i[j]=0;
else C[j]=A[j]/(-A);
}

int z=0,c=0;
while(z!=n)
{
for(i=0;i=X;
double t=0;
for(i=0;i[j]*X[j];
X=t+D;
t=0;
}
z=0;
for(i=0;i-Y<<" ";
cout<<"\nKolicestvo iteratii = "<[j]<<"X"<[j+1]<0);
else if(j!=2)cout<<"+";
t=t+A[j]*X[j];
}
cout<<"="<

t=t+A[j]*X[j];
}
cout<<"="<

t=t+A[j]*X[j];
}
cout<<"="<