Gauss Jecobi program

#include<iostream>
using namespace std;
int main()
{

    int n,i,j,it=0;
    cin>>n;
    float mat[n][n+1],x[50]={0},y[50],sum=0;
    for(i=0;i<n;i++)
    {
        for(j=0;j<n+1;j++)
        {
            cin>>mat[i][j];
        }
    }
    do{

    for(i=0;i<n;i++)
    {
       for(j=0;j<n;j++)
        {
            if(i==j)
            {
                continue;
            }
            else
            {
                sum=sum+mat[i][j]*x[j];
            }
        }
        y[i]=(mat[i][n]-sum)/mat[i][i];
        sum=0;
    }
    for(i=0;i<n;i++)
    {
        x[i]=y[i];
    }
    ++it;
    }while(it<=20);
    for(i=0;i<n;i++)
    {
        cout<<x[i]<<endl;
    }
}