Сможете дописать? Высчитать элементы квадратной матрицы и пр.

//---------------------------------------------------------------------------
 
#include <vcl.h>
#pragma hdrstop
#include <math.h>
#include "Unit1.h"
 
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//глобальные параметры
const int n=4;
double a[5][5];
double x[5];
//подпрограмма для элементов матрицы. Нужно исправить
void  elem_matr( )
{ int i,j;
    for(i=0;i<=n;i++)
    for(j=0;j<=n;j++)
    a[i][j]=pow(3,-j+1)*(i-((j+1)/(3)))*((i+3)/(j+6)); //тут вот ошибка гдето. почти все на ноль выводит
    }
//подпрограмма элементов вектора. Вообще не знаю как написать по такому алгоритму
 void elem_vectora( )
 {int i,j;
 
 
 
 
 
 
 }
 
 
 //подпрограмма функции G. Выбивает ошибку в программе. Потому, что я краб.
 double fun_G( )
 { int i,k; float S,P;
 for(i=0, P=1; i<=4; i++) {
  for(k=0, S=1; k<=i; k++) S+=x[k];
  P *= (x[i] + S);
}
return (P)/(x[i]) ;
}
 
//---------------------------------------------------------------------------
 
 
void __fastcall TForm1::Button1Click(TObject *Sender)
{
 int i,j;
 //вызов подпрограммы для матрицы
 elem_matr( );
  for(i=0;i<=n;i++)
  {
  
  StringGrid1->Cells[0][i+1]=IntToStr(i+1)+" -й ряд";
  }
  for(j=0;j<=n;j++)
  {
  
  StringGrid1->Cells[j+1][0]=IntToStr(j+1)+" -й столбик";
  }
  
  for(i=0;i<=n;i++)
  for(j=0;j<=n;j++)
  StringGrid1->Cells[j+1][i+1]=FormatFloat("0.000",a[i][j]);
}
//---------------------------------------------------------------------------
void __fastcall TForm1::Button2Click(TObject *Sender)
{
int i,j;
//вызов подпрограммы вектора
elem_vectora( );
for(j=0;j<=n;j++)
{
StringGrid2->Cells[0][j]=IntToStr(j+1)+" -й элемент";
StringGrid2->Cells[1][j]=FormatFloat("0.00",x[j]);
}
}
//---------------------------------------------------------------------------
void __fastcall TForm1::Button3Click(TObject *Sender)
 
//вычисление G
{
int i; double G;
G=fun_G( );
Edit1->Text=FormatFloat("0.000",G);
}
//---------------------------------------------------------------------------
void __fastcall TForm1::Button4Click(TObject *Sender)
{
 int i,j;
    for(i=0;i<=n;i++)
    for(j=0;j<=n;j++)
     StringGrid1->Cells[j+1][i+1]=" ";
     for(i=0;i<=n;i++)
     StringGrid2->Cells[1][i+1]=" ";
     Edit1->Text=" ";
}
//---------------------------------------------------------------------------
void __fastcall TForm1::Button5Click(TObject *Sender)
{
Close();
}
//---------------------------------------------------------------------------