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

#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;
}