перевод на С

//---------------------------------------------------------------------------
 
#include <vcl.h>
#pragma hdrstop
 
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
int timer, timer_1, r, fi, x;
double T[9][359][1200];
double dT[9][359][1200];
double a, lam1=0.00038, lam2=0.00003026, N;
//---------------------------------------------------------------------------
 
__fastcall TForm1::TForm1(TComponent* Owner)
    : TForm(Owner)
{
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::Edit1KeyPress(TObject *Sender, System::WideChar &Key)
{
if (Key>='0' && Key<='9' || Key=='\b')
  {
  if (Edit1->Text=="")
    {
    if (Key=='0') {Key=0;}
    timer=0;
    }
  else timer=StrToInt(Edit1->Text);
  }
else Key=0;
}
//---------------------------------------------------------------------------
 
   void __fastcall TForm1::Button1Click(TObject *Sender)
{
for (r=1; r<=8; r++)
 for (fi=1; fi<=358; fi++)
  for (x=1; x<=1199; x++)
     {
      T[r][fi][x]=30;
      dT[r][fi][x]=30;
     }
 
for (timer_1=1; timer_1<=timer; timer_1++)
   {
   for (r=1; r<=8; r++)
    for (fi=1; fi<=358; fi++)
     for (x=1; x<=1199; x++)
        {
        if ((r<5)||(r>8)) {a=22.16;}
        else a=0.3769;
        dT[r][fi][x]=(T[r-1][fi][x]-2*T[r][fi][x]+T[r+1][fi][x])/1;
        dT[r][fi][x]=dT[r][fi][x]+(1/4.5*((T[r][fi][x]-T[r-1][fi][x])/1));
        dT[r][fi][x]=dT[r][fi][x]+(1/20.25*((T[r][fi-1][x]-2*T[r][fi][x]+T[r][fi+1][x])/1));
        dT[r][fi][x]=dT[r][fi][x]+((T[r][fi][x-1]-2*T[r][fi][x]+T[r][fi][x+1])/1);
        dT[r][fi][x]=a*0.01*dT[r][fi][x];
        }
 
   for (r=0; r<=9; r++)
    for (fi=0; fi<=359; fi++)
        {
        T[r][fi][0]=1400;
        T[r][fi][1200]= T[r][fi][1199];
        }
 
   for (fi=0; fi<=359; fi++)
    for (x=1; x<=1199; x++)
       {
       T[0][fi][x]=T[1][fi][x];
       T[5][fi][x]=(lam1*T[6][fi][x]-lam2*T[4][fi][x])/(lam1+lam2);
       T[9][fi][x]=T[4][fi][x];
       T[8][fi][x]=(lam1*T[7][fi][x]-lam2*T[9][fi][x])/(lam1+lam2);
       }
 
   for (r=1; r<=8; r++)
    for (fi=1; fi<=358; fi++)
     for (x=1; x<=1199; x++)
        {
        T[r][fi][x]=T[r][fi][x]+dT[r][fi][x];
        }
 
   N=T[6][357][1198];
   Series1->AddXY(timer_1*0.01, N);
   ProgressBar1->Position=100*timer_1/timer;
   }
Label11->Caption= "T[7][357][1198] = " + FloatToStr(N);
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::Edit1Change(TObject *Sender)
{
timer=StrToInt(Edit1->Text)/0.01;
}
//---------------------------------------------------------------------------
 
void __fastcall TForm1::FormShow(TObject *Sender)
{
timer=StrToInt(Edit1->Text)/0.01;
}
//---------------------------------------------------------------------------