Переписать с Паскаля на С++

program PROGRAMMA7;
uses crt, graph;
const n=100; m=100; h=1; dt=0.2;
var  ii,jj,kk,i,j,DriverVar, ModeVar, ErrorCode : integer;
t: array[1..N, 1..M] of real; 
q,a,b,bb :real; naprav,uslovie: boolean;
 
procedure Init; {---- Инициализация графики ----}
begin DriverVar:=Detect; 
      InitGraph(DriverVar,ModeVar,'c:\bp\bgi');
      ErrorCode:=GraphResult;   
      if ErrorCode <> grOK then Halt(1); end;
 
procedure Param_sred; {---Коэффициент температуропроводности---}
begin if j<70 then a:=2 else a:=1; end;
 
procedure Istoch;     {--- Источники тепла ---}
begin if ((i>45)and(i<55))and((j>70)and(j<75)) 
                                 then q:=50 else q:=0; end;
 
procedure Nach_uslov; {-- Начальное распределение температуры --}
begin  For i:=1 to N do For j:=1 to M do
  begin  uslovie:=((j<65)and(j>45)and(i>20)and(i<30))
                or((j<45)and(j>35)and(i>50)and(i<60));
  if uslovie=true then t[i,j]:=450 else t[i,j]:=1;  end; 
end;
 
procedure Raschet; {---- Расчет температуры ----}
begin Istoch; Param_sred;
    t[i,j]:=t[i,j]+a*(t[i,j+1]-2*t[i,j]+t[i,j-1])*dt/(h*h)+q;
   if naprav=true then t[i,j]:=
         t[i,j]+a*(t[i+1,j]-2*t[i,j]+t[i-1,j])*dt/(h*h); 
end;
 
procedure Draw;{---- Вывод на экран ----}
begin if t[i,j]>50 then setcolor(2); 
    if t[i,j]>300 then setcolor(12);
    if (t[i,j]<300)and(t[i,j]>120) then setcolor(10);
    if (t[i,j]<120)and(t[i,j]>70) then setcolor(3);
    if (t[i,j]<70)and(t[i,j]>30) then setcolor(4);
    if (t[i,j]<30)and(t[i,j]>20) then setcolor(5);
    if (t[i,j]<20)and(t[i,j]>10) then setcolor(7);
    if t[i,j]<10 then setcolor(15); 
    rectangle(i*5+50,j*5,i*5+54,j*5+4);   
end;              
 
BEGIN             {---- Основная программа ----}
Init; Nach_uslov;
Repeat kk:=kk+1;
For i:=2 to N-1 do 
   For j:=2 to M-1 do begin naprav:=false; Raschet; end;
For j:=2 to M-1 do 
   For i:=2 to N-1 do begin naprav:= true; Raschet; end;
For i:=2 to N-1 do For jj:=2 to M-1 do  
      begin j:=M+1-jj; naprav:=true; Raschet; end;
For j:=2 to M-1 do For ii:=2 to N-1 do
     begin i:=N+1-ii; naprav:=false; Raschet; end;
if kk/2=round(kk/2) then 
   For i:=2 to N-1 do For j:=2 to M-1 do Draw;
until KeyPressed; CloseGraph;
END.