Метод Зейделя с++

[Fackey]

Здравствуйте ) Я хотел бы обратиться к вам за помощью по поводу " Решения слау методом Зейделя". Мб у вас есть рабочий исходник этой программы , т.к в поиске я ничего не нашел (включая гугл) . Буду очень благодарен за рабочую программу .


Вот есть одна , но она неработает , мб вы что-нибудь сможете с ней сделать? (

#include <math.h>
#include <stdio.h>
#include <iostream.h>

const int nn=100;
typedef double arr[nn][nn];
typedef double vect[nn];

void input(int*, double *);
bool exact(const double *px, const double *pp, const int *n);

arr A;
double e;

int main()
{

int n,i,j;
double sum;
vect x,p,B;

input(&n,B); // input n,A,B.
for (i=0; i<n; i++){ // initial approximatization
x[i]=1;
p[i]=1;
}

do {
for (i=0; i<n; i++){
sum = 0;
for (j=0; j<n; j++)
if (i!=j) sum+=x[j]*A[i][j];
// printf("sum= %.2f",sum);
// cout << "i = " << i;
p[i]=x[i];
x[i]=(B[i] - sum)/A[i][i];
/* for (j=0; j<n; j++)
printf("%.2f ",x[j]);
cout << exact(x,p,&n);*/

}
} while (!exact(x,p, &n)) ;


for (j=0; j<n; j++)
printf("X%d = %.3f\n",j+1,x[j]);

puts("Verification:");
for (i=0; i<n; i++){
sum = 0;
for (j=0; j<n; j++){
sum+=x[j]*A[i][j];
if ((j==0) && (A[i][j]*x[j]<0)) cout <<"-";
printf("%.2f*%.2f",fabs(A[i][j]),fabs(x[j]));
if ((j!=n-1)&&(A[i][j+1]*x[j]>=0)) cout <<" + ";
else if ((j!=n-1)&&(A[i][j+1]*x[j]<0)) cout <<" - ";
else {cout << " = ";
printf("%.3f ;\n", sum);
}
}
}

}

void input(int *m, double *PtrB)
{
int i,j;

cout << "Enter e.\n";
cin >> e;
cout << "\nEnter n.\n";
cin >> *m;
cout << "\nEnter matrix A.\n";
for (i=0; i<*m; i++)
for (j=0; j<*m; j++)
cin >> A[i][j];
cout << "\nEnter matrix B.\n";
for (i=0; i<*m; i++)
cin >> *(PtrB+i);

/*for (i=0; i<*m; i++){
cout << endl;
for (j=0; j<*m; j++)
printf("%.2f ",A[i][j]);}*/

}

bool exact(const double *px, const double *pp, const int *n)
{
int i;
bool f=true;

for (i=0; i<*n; i++)
if ( fabs( *(px+i) - *(pp+i) ) > e){
f=false;
break;
}
return f;
}

Буду очень благодарен за помощь .

[NiCola999]

Задолбали выкладывать код как попало! Заносите его в специальный тег CODE, чтобы были видны отступы и как следует отформатируйте код

[Vago]

Сообщение от Fackey

в поиске я ничего не нашел (включая гугл) .

http://www.google.com/search?rls=en&...C%2B%2B+source
Лично мне понравилась шестая ссылка.

[Fackey]

Vago спасибо огромное . Ты действительно мне очень помог )