program na inverznu maticu v C++

#include <vcl.h>
#pragma hdrstop

#include "Unit1.h"
#include <fstream.h>
#include <stdlib.h>
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma link "CSPIN"
#pragma resource "*.dfm"
TForm1 *Form1;
const  poc=4;
int M[poc][poc];
int i,j;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
void __fastcall TForm1::Button1Click(TObject *Sender)
{
int x= StrToInt(CSpinEdit1->Text);
 Memo1->Lines->Clear();
 for (int i=0; i<x; i++)
      for (int j=0; j<x; j++) M[i][j]=random(100);
 for (int i=0; i<x; i++)
 { AnsiString Riadok= " ";
   for (int j=0; j<x; j++)
   { Riadok= Riadok + AnsiString (M[i][j]) + '\t';
   }
 Memo1->Lines->Add(Riadok);
 }
}

float determinant;
determinant=Matica[0][0]*Matica[1][1]*Matica[2][2]+Matica[1][0]*Matica[2][1]*Matica[0][2]+Matica[2][0]*Matica[0][1]*Matica[1][2]-Matica[2][0]*Matica[1][1]*Matica[0][2]-Matica[1][0]*Matica[0][1]*Matica[2][2]-Matica[0][0]*Matica[2][1]*Matica[1][2];
float a,b,c,d,e,f,g,h,i;
a=pow(-1,2)*(Matica[1][1]*Matica[2][2]-Matica[2][1]*Matica[1][2]);
b=pow(-1,3)*(Matica[1][0]*Matica[2][2]-Matica[2][0]*Matica[1][2]);
c=pow(-1,4)*(Matica[1][0]*Matica[2][1]-Matica[2][0]*Matica[1][1]);
d=pow(-1,3)*(Matica[0][1]*Matica[2][2]-Matica[2][1]*Matica[0][2]);
e=pow(-1,4)*(Matica[0][0]*Matica[2][2]-Matica[2][0]*Matica[0][2]);
f=pow(-1,5)*(Matica[0][0]*Matica[2][1]-Matica[2][0]*Matica[0][1]);
g=pow(-1,4)*(Matica[0][1]*Matica[1][2]-Matica[1][1]*Matica[0][2]);
h=pow(-1,5)*(Matica[0][0]*Matica[1][2]-Matica[1][0]*Matica[0][2]);
i=pow(-1,6)*(Matica[0][0]*Matica[1][1]-Matica[1][0]*Matica[0][1]);

Matica3[0][0]=(1/determinant)*a;
Matica3[1][0]=(1/determinant)*b;
Matica3[2][0]=(1/determinant)*c;
Matica3[0][1]=(1/determinant)*d;
Matica3[1][1]=(1/determinant)*e;
Matica3[2][1]=(1/determinant)*f;
Matica3[0][2]=(1/determinant)*g;
Matica3[1][2]=(1/determinant)*h;
Matica3[2][2]=(1/determinant)*i;
AnsiString riadok1,riadok2,riadok3;
riadok1=FloatToStr(Matica3[0][0])+'\t'+FloatToStr(Matica3[0][1])+'\t'+FloatToStr(Matica3[0][2]);

riadok2=FloatToStr(Matica3[1][0])+'\t'+FloatToStr(Matica3[1][1])+'\t'+FloatToStr(Matica3[1][2]);
riadok3=FloatToStr(Matica3[2][0])+'\t'+FloatToStr(Matica3[2][1])+'\t'+FloatToStr(Matica3[2][2]);
Memo3->Clear();
Memo3->Lines->Add(AnsiString(R)+" "+AnsiString(S));
Memo3->Lines->Add(riadok1);
Memo3->Lines->Add(riadok2);
Memo3->Lines->Add(riadok3);