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

#include<iostream>
#include<fstream>
#include<math.h>
using namespace std;
int main(){
 
freopen("inputK.txt", "rt", stdin);
freopen("outputK.txt", "wt", stdout);
    double last=0,s1,s2,eps,err, t=0;
    double *x=new double [6];
    double *p=new double [6];
    double *b=new double [6];
    double **a=new double *[6];
    cin>>eps;
    for(int i=0;i<6;i++)
        a[i]=new double[6];
    for(int i=0;i<6;i++){
        for(int j=0;j<6;j++)
            cin>>a[i][j];
        cin>>b[i];
    }
    cout<<"array of ~"<<endl;
    for(int i=0;i<6;i++){
        p[i]=b[i];
        for (int j=i+1;j<6;j++){
            t+=-a[i][j];
        }
 
        p[i]+=t;
 
        p[i]/=a[i][i];
        cout<<"p["<<i+1<<"]="<<p[i]<<endl;
    }
 
    cout<<endl;
    for(int i=0;i<6;i++) x[i]=0;
    int iter=0;
 
    do{
        err=0;
        for(int i=0;i<6;i++){
            s1=0;
            s2=0;
            for(int j=0;j<i;j++)
                s1+=a[i][j]*x[j];
            for(int j=i;j<6;j++)
                s2+=a[i][j]*x[j];
            last=x[i];
            x[i]-=(1/a[i][i])*(s1+s2-b[i]);
            if (fabs(last-x[i])>err) err=fabs(last-x[i]);
        }
        iter++;
    }
    while(err>=eps);
    cout<<"Number of iteration "<<iter<<endl;
    for(int i=0;i<6;i++)
        cout<<"x["<<i+1<<"]="<<x[i]<<endl;
    for(int i=0;i<6;i++)
        delete [] a[i];
    delete [] a;
    delete [] b;
    delete [] x;
    return 0;
}