Gauss

#include <iostream>
#include <fstream>

using namespace std;
const int n=4;
int iter=50;
int i,j,k;
double sum1,sum2;
double w[n][n];
double z[n];
double x[n]={1,1,1,1};
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];
        }
        x[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];
            }
            sum2=0;
            for(int j=i+1;j<n;j++){
                sum2=sum2+w[i][j]*x[j];
            }
            x[i]=sum1+sum2+z[i];
        }
        sum2=0;
        for(int j=0;j<n-1;j++){
            sum2=sum2+w[n-1][j]*x[j];
        }
        x[n-1]=sum2+z[n-1];
        cout<<"iteracja numer:"<<k+1<<endl;
        for(int i=0;i<n;i++){

    cout<<x[i]<<endl;
}
    }

    return 0;
}