relaksacyjny

#include <iostream>
#include <fstream>

using namespace std;
const int n=4;
int iter=30;
int i,j,k;
double sum1,sum2;
double w[n][n];
double z[n];
double x[n]={3,6,4,3};
double x1[n];
float a=1.8;
int main()
{
    ifstream plik("M.txt");
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            plik>>w[i][j];
        }
    }
    plik.close();
    ifstream plik1("V.txt");
    for(int i=0;i<n;i++){
        plik1>>z[i];
    }
    plik1.close();
//poczatek
    for(int k=0;k<iter;k++){
        sum1=0;
        for(int j=1;j<n;j++){
            sum1=sum1+w[0][j]*x[j];
        }
        x1[0]=sum1+z[0];
        for(int i=1;i<n-1;i++){
            sum1=0;
            for(int j=0;j<i;j++){
                sum1=sum1+w[i][j]*(x[j]+a*(x1[j]-x[j]));
            }
            sum2=0;
            for(int j=i+1;j<n;j++){
                sum2=sum2+w[i][j]*x[j];
            }
            x1[i]=sum1+sum2+z[i];
        }
        sum2=0;
        for(int j=0;j<n-1;j++){
            sum2=sum2+w[n-1][j]*(x[j]+a*(x1[j]-x[j]));
        }
        x1[n-1]=sum2+z[n-1];
    for(int i=0;i<n;i++){
        x[i]=x1[i];
    }
    cout<<"interacja:"<<k+1<<endl;
for(int i=0;i<n;i++){
    cout<<x1[i]<<endl;

}
}
    return 0;
}