Метод Гаусса Visual С++

    #include <iostream>
#include <iomanip>
#include <conio.h>
#include <math.h>
using namespace std;
int F;
int main()
 {
int n=3,i,j,k;
  float a[9][9],b[6],x[6],h;
 cout<<"Enter Array:"<<endl;
 for (i=0;i<n;i++)
 {
     for (j=0;j<n;j++)
     {
     cin>>a[i][j];
     }
 }
 cout<<endl<<"Enter C:"<<endl;
 for (i=0;i<n;i++)
 {
 cin>>b[i];
 }
  for (i=0;i<n;i++)
 {
 cout<<a[i][0]<<"x1+"<<a[i][1]<<"x2+"<<a[i][2]<<"x3= "<<b[i]<<endl;
 }
 for (i=0;i<n;i++) //возможно неправильно считает
 {
     for (j=i+1;j<n;j++)
     {
         a[j][i]=(n-4)*a[j][i]/a[i][i];
         for (k=i+1;k<n;k++)
         {
         a[j][k]+=a[j][i]*a[i][k];
         b[j]+=a[j][i]*b[i];
         }
     }  
 }
     x[2]=b[2]/a[2][2];
 for (i=n-1;i>=0;i--) //возможно здесь
{
 h=b[i];
 for (j=i+1;j<n;j++)
    h-=x[j]*a[i][j];
    x[i]=h/a[i][i];
}
 
 cout<<"Solution: "<<endl;
 for (i=0;i<n;i++)
 {
     std::cout << std::fixed << std::setprecision(3)<<"X["<<i<<"]="<<x[i]<<std::endl;
 }
 for(i=0;i<n;i++) //просто проверка
 {
     F+=a[0][i]*x[i];
 }
 cout<<endl<<"Otvet:"<<F<<endl;
 system("pause");
 return 0;
}

 x[2]=b[2]/a[2][2];

x[6]=b[2]/a[2][2];

x[2]=b[2]/a[2][2]