Дайте комментарии к коду

#include "stdafx.h"
#include<iostream>
#include<conio.h>
#include<stdio.h>
 
using namespace std;
 
double a[15][16];
int n,x,y;
 
void input()
{
  cout<<"Razmernost': ";
  cin>>n;
  for(int i=0;i<n;i++)
  {
    for(int j=0;j<n;j++)
    {
      cout<<"a["<<i+1<<"]["<<j+1<<"]=";
      cin>>a[i][j];
    }
    cout<<"b["<<i+1<<"]=";
    cin>>a[i][n];
  }
}
 
void output()
{
  for(int i=0;i<n;i++)
  {
  cout<<"(";
    for(int j=0;j<n;j++)
      printf("% 6.1f",a[i][j]);
  cout<<" |";
  printf("% 6.1f",a[i][n]);
  cout<<")\n";
  }
}
 
void main()
{
input();
cout<<"\n";
output();
for(int i=0;i<n;i++)
{
  //выбор разрешающего эл-та
  for(int j=0;j<n;j++)
  {
    if(a[i][j]!=0)
    {
      y=j;
      break;
    }
    if(j==n-1&&a[i][n]!=0)
    {
      cout<<"Sistema ne sovmestna";
      getch();
      return;
    }
  }
 
 
  for(int u=0;u<n;u++)
  for(int j=0;j<n+1;j++)
    if(j!=y&&u!=i)a[u][j]=(a[i][y]*a[u][j]-a[i][j]*a[u][y])/a[i][y];
  for(int j=0;j<n+1;j++)
    if(j!=y)a[i][j]=a[i][j]/a[i][y];
  a[i][y]=1;
  for(int u=0;u<n;u++)
    if(u!=i)a[u][y]=0;
}
cout<<"\n";
output();
cout<<"\n";
for(int i=0;i<n;i++)
{
  for(int j=0;j<n;j++)
    if(a[i][j]!=0)cout<<"x["<<j+1<<"]=";
  printf("% 15.10f",a[i][n]);
  cout<<"\n";
}
getch();
}