Нужно перевести данную программу из Turbo Pascal в Delphi

Program Statist2
Uses Crt,Graph; Const k=10;
Var i,j,n,m1,k1,k2,k3,l,r,v,j1,j2,
    dx1,dx2,dy1,dy2,X0,Y0,YF0     :Integer;
    ErrorNumber,DriverVar,
    ModeVar,Regim                 :Integer;
    r1,r2,Xmin,Xmax,fmax,lambda,mx,Dx,sigma,hi_kvad,p_kvad,
    s1,s2,s3,s4,s5,h,dh,alfa,beta,x1,x2,y1,y2         :Real;
    x                             :Array[1..100] of Real;
    xt,ft,Ftx                     :Array[0..100] of Real;
    p_hi                          :Array[1..7] of Real;
    r_hi                          :Array[1..20,1..7] of Real;
    p,p_r,f                       :Array[1..k] of Real;
    m,fg                          :Array[1..k] of Integer;
    Fx                            :Array[1..(k+1)] of Real;
    Frx                           :Array[0..k] of Real;
    y                             :Array[1..k,1..50] of Real;
    St,St1,St2,St3,St4,St5        :String;
    PR                            :Char;
    fd,fr,fhi                     :Text;
BEGIN ClrScr;{ БЛОК ФАЙЛОВОЙ СИСТЕМЫ ВВОДА СВ X }
Assign(fd,'302-1S');
    Reset(fd); Read(fd,n); For i:=1 to n do Read(fd,x[i]);
    Close(fd);{БЛОК ФАЙЛОВОЙ СИСТЕМЫ ВВОДА "ХИ-КВАДРАТ" }
    Assign(fhi,'tab_hi.pas');
    Reset(fhi); For j:=1 to 7 do Read(fhi,p_hi[j]);
    For i:=1 to 20 do
    For j:=1 to 7 do
    Read(fhi,r_hi[i,j]);
    Close(fhi);{ РАСЧЕТ СТАТИСТИЧЕСКОГО СРЕДНЕГО }
    s1:=0; For i:=1 to n do s1:=s1+x[i]; mx:=s1/n;
    { РАСЧЕТ СТАТИСТИЧЕСКОЙ ДИСПЕРСИИ }
    s2:=0; For i:=1 to n do s2:=s2+Sqr(x[i]-mx); Dx:=s2/n;
    { РАСЧЕТ СТАТИСТИЧЕСКОГО СКО }
    sigma:=Sqrt(Dx);{БЛОК СОРТИРОВКИ ЗНАЧЕНИЙ СВ X ПО ВОЗРАСТАНИЮ}
     For k1:=n downto 1 do
             Begin m1:=1; For i:=1 to k1 do
                 If x[i] >= x[m1] Then m1:=i; r1:=x[k1];
                 x[k1]:=x[m1]; x[m1]:=r1
         End;
    Xmin:=x[1]; Str(Xmin,St4); Xmax:=x[n]; Str(Xmax,St5);
    h:=(Xmax-Xmin)/k;{БЛОК СОРТИРОВКИ СВ ПО РАЗРЯДАМ И РАСЧЕТ
            СТАТИСТИЧЕСКОЙ ФУНКЦИИ РАСПРЕДЕЛЕНИЯ }
    For j:=1 to k do
    Begin
        k2:=0; For i:=1 to n do
        if (x[i]>=(Xmin+(j-1)*h)) and (x[i]<(Xmin+j*h)) then
           Begin k2:=k2+1; y[j,k2]:=x[i]
       End;
       m[j]:=k2
        End;{ФАЙЛОВАЯ СИСТЕМА ЗАПИСИ РЕЗУЛЬТАТОВ РАСЧЕТОВ}
Assign(fr,'302rz1');Rewrite(fr);
Writeln(fr,'Аппаратура П-302, 1-е сечение');
 
 
    Writeln(fr,'Выборка N= 53 ');
    Writeln(fr,'Статистич. среднее M[y]=',mx:5:4);
    Writeln(fr,'Статистич. дисперсия D[y]=',Dx:5:4);
    Writeln(fr,'Статистич. СКО sigma=',sigma:5:4);
        Writeln(fr);Writeln(fr);
    Writeln(fr,'Nразр m[j]   P[y]    f[y]    F[y]');
{БЛОК РАСЧЕТА ЧАСТОТ РАЗРЯДОВ И СТАТИСТИЧ. ФУНКЦИИ РАСПРЕДЕЛЕНИЯ}
    s3:=0;Fx[1]:=0.0; For j:=1 to k do
    Begin
        p[j]:=m[j]/n; f[j]:=p[j]/h; s3:=s3+p[j]; Fx[j+1]:=s3;
    Writeln(fr,'j=',j:2,' | ',m[j]:2,' |',
               p[j]:6:3,' | ', f[j]:6:4,' |', Fx[j+1]:6:3)
    End;
    Writeln(fr);
        fmax:=f[1]; For j:=1 to k do If f[j] >= fmax then
        fmax:=f[j];{ИНИЦИАЛИЗАЦИЯ ГРАФИЧЕСКОГО РЕЖИМА }
    DriverVar:=Detect; ModeVar:=Regim;
    InitGraph(DriverVar,ModeVar,'');
    ErrorNumber:=GraphResult; If ErrorNumber <> grOk Then
       Writeln(GraphErrorMsg(ErrorNumber)) Else
Begin { УСТАНОВКА ЦВЕТА ФОНА }
    SetBkColor(15);{ПОСТРОЕНИЕ РАМКИ ЭКРАНА }
    SetColor(LightMagenta);{РАСЧЕТ МАСШТАБОВ }
    dx1:=Round((GetMaxX-490)/k);
    dx2:=Round((GetMaxX-490)/(Xmax-Xmin));
    dy1:=Round((GetMaxY-250)/fmax); dy2:=Round((GetMaxY-250)/10);
    {----- ПОСТРОЕНИЕ ОСИ ОРДИНАТ:
    F(y) и f(y);F5*0(y) и f5*0(y)}
    SetColor(8); Line(35,185,35,GetMaxY-45);
    OutTextXY(7,167,'*'); OutTextXY(0,175,'F (y)');
    OutTextXY(0,186,'F(y)'); OutTextXY(53,167,'*');
    OutTextXY(45,175,'f (y)'); OutTextXY(45,186,'f(y)');
    For j:=0 to 10 do
        Begin
        PutPixel(35,GetMaxY-45-j*dy2,Red);
        PutPixel(36,GetMaxY-45-j*dy2,Red);
        Str(j,St); If ( j>0 ) and ( j<10 ) Then
        OutTextXY(7,GetMaxY-50-j*dy2,'0.'+St);
        OutTextXY(22,GetMaxY-50,'0')
        End;{ ПОСТРОЕНИЕ ОСИ АБСЦИСС: СВ y }
    Line(35,GetMaxY-45,GetMaxX-435,GetMaxY-45);
    OutTextXY(GetMaxX-428,GetMaxY-48,'y'); For j:=1 to (k+1) do
        Begin Str(j,St);
        End;
              GetMaxY-35,'my');}
    Line(30+dx1+Round((mx-Xmin)*dx2),GetMaxY-40,
         30+dx1+Round((mx-Xmin)*dx2),GetMaxY-45);
    {----- ПОСТРОЕНИЕ ГИСТОГРАММЫ РАСПРЕДЕЛЕНИЯ ------------}
    For j:=1 to k do
        Begin fg[j]:=Round(f[j]*dy1);
        Rectangle(30+j*dx1,GetMaxY-45-fg[j],
                  30+(j+1)*dx1,GetMaxY-45);
        SetFillStyle(0,j+1); Bar(30+j*dx1+1,GetMaxY-45-fg[j]+1,
            30+(j+1)*dx1-1,GetMaxY-45-1);
        Line(45,GetMaxY-45,45,GetMaxY-40);
        Line(195,GetMaxY-45,195,GetMaxY-40);
        End;
    OutTextXY(0,155,'Аппаратура П302,1-е сечение');
    SetColor(Blue); PR:=ReadKey;
    {----- ГРАФИК СТАТИСТИЧЕСКОЙ ФУНКЦИИ РАСПРЕДЕЛЕНИЯ -----}
 
 
    SetColor(Green); For j:=1 to k do
        Line(30+j*dx1,GetMaxY-45-Round(Fx[j]*10*dy2),
             30+(j+1)*dx1,GetMaxY-45-Round(Fx[j+1]*10*dy2));
    Delay(1000);{ГРАФИКИ ТЕОРЕТИЧЕСКОЙ ПЛОТНОСТИ РАСПРЕДЕЛЕНИЯ
           И ФУНКЦИИ РАСПРЕДЕЛЕНИЯ -------------------------}
    X0:=0;Case PR of
        '1':Begin lambda:=1.0/sigma;
            WriteLn(fr,'Экспоненциальное распределение');
            WriteLn(fr,'lambda=',lambda:8:6); WriteLn(fr);
            Y0:=Round(lambda*dy1) End;
        '2':Begin alfa:=mx-Sqrt(3*Dx);
            beta:=mx+Sqrt(3*Dx);
            WriteLn(fr,'Распределение равномерной плотности');
            Writeln(fr,'alfa=',alfa:6:2,'  ',
            'beta=','  ',beta:6:2);
            Y0:=Round(1.0/(beta-alfa)*dy1) End;
        '3':Begin WriteLn(fr,'Нормальное распределение');
           { WriteLn(fr,'sigma=',sigma:5:2);}
            Y0:=Round(1.0/Sqrt(2*Pi)/sigma*
                Exp(-Sqr(Xmin-mx)/2/Sqr(sigma))*dy1) End;
    End;
    YF0:=0; s4:=0; j:=0; For i:=0 to 100 do
        Begin
        Case PR of
            '1':Begin dh:=(Xmax-Xmin)/100;r:=k-2;
                xt[i]:=Xmin+i*dh;
                ft[i]:=lambda*Exp(-lambda*(xt[i]-Xmin))
                End;
            '2':Begin dh:=(beta-alfa)/100;r:=k-3;
                xt[i]:=alfa+i*dh;
                ft[i]:=1.0/(beta-alfa)
                End;
 
            '3':Begin dh:=(Xmax-Xmin)/100;r:=k-3;
                xt[i]:=Xmin+i*dh;
                ft[i]:=1.0/Sqrt(2*Pi)/sigma*
                       Exp(-Sqr(xt[i]-mx)/2/Sqr(sigma))
                End
        End;
        s4:=s4+ft[i]*dh; Ftx[i]:=s4;
    { РАСЧЕТ ФУНКЦИИ РАСПРЕДЕЛЕНИЯ НА ГРАНИЦАХ РАЗРЯДОВ}
        If i=Round(j*100/k) then
           Begin
           Frx[j]:=Ftx[i]; j:=j+1
           End;{ГРАФИК ТЕОРЕТИЧЕСКОЙ ПЛОТНОСТИ РАСПРЕДЕЛЕНИЯ }
        SetLineStyle(SolidLn,0,NormWidth);
        SetColor(Blue);
        Line(30+dx1+X0,GetMaxY-45-Y0,30+dx1+Round((xt[i]-Xmin)*dx2),
             GetMaxY-45-Round(ft[i]*dy1));
    { ГРАФИК ТЕОРЕТИЧЕСКОЙ ФУНКЦИИ РАСПРЕДЕЛЕНИЯ}
        SetColor(Red);Line(30+dx1+X0,GetMaxY-45-YF0,
             30+dx1+Round((xt[i]-Xmin)*dx2),
             GetMaxY-45-Round(Ftx[i]*10*dy2));
        X0:=Round((xt[i]-Xmin)*dx2); Y0:=Round(ft[i]*dy1);
        YF0:=Round(Ftx[i]*10*dy2)
        End;
    {ОПРЕДЕЛЕНИЕ СООТВЕТСТВИЯ ТЕОРЕТИЧЕСКОГО И СТАТИСТИ-
         ЧЕСКОГО РАСПРЕДЕЛЕНИЙ.КРИТЕРИЙ К.ПИРСОНА }
         s5:=0; For j:=1 to k do
        Begin
 
 
        SetColor(DarkGray);
        p_r[j]:=Frx[j]-Frx[j-1];
        s5:=s5+Sqr(m[j]-n*p_r[j])/n/p_r[j]; hi_kvad:=s5
        End;
    If (hi_kvad > r_hi[r,1]) and (hi_kvad < r_hi[r,7]) then
        Begin
        OutTextXY(53,GetMaxY-10,'Рабочая гипотеза'); j1:=1;
        While (hi_kvad > r_hi[r,j1]) do
              Begin
                 x1:=p_hi[j1]; y1:=r_hi[r,j1]; j1:=j1+1
              End;
        j2:=7;
        While (hi_kvad < r_hi[r,j2]) do
              Begin
                 x2:=p_hi[j2]; y2:=r_hi[r,j2]; j2:=j2-1
              End;
        p_kvad:=x1+(x2-x1)/(y2-y1)*(hi_kvad-y1);
                 Writeln(fr,'r=',r,'  ',
                   'hi_kvad=',hi_kvad:2:2,' ',
                   'p_kvad=',p_kvad:2:2);
                 Writeln(fr)
       End
       Else
       Begin
       OutTextXY(51,GetMaxY-10,'Гипотеза неверна');
                  Writeln(fr,'r=',r,'  ',
                 'hi_kvad=',hi_kvad:2:2,'  ',
                 'p_kvad=0.0001');
                 Writeln(fr)
        End;
    Close(fr);{ ВЫВОД ОБОБЩАЮЩЕЙ ИНФОРМАЦИИ}
    Case PR of
        '1':OutTextXY(150,GetMaxY-20,
                      'Экспоненциальное распределение');
        '2':Begin
            OutTextXY(20+dx1+Round((alfa-Xmin)*dx2)-16,
                      GetMaxY-35,'alfa');
            OutTextXY(20+dx1+Round((beta-Xmin)*dx2)-16,
                      GetMaxY-35,'beta');
            Line(20+dx1+Round((alfa-Xmin)*dx2),GetMaxY-45,
                 20+dx1+Round((alfa-Xmin)*dx2),
                            GetMaxY-45-Round(ft[i]*dy1));
            Line(20+dx1+Round((beta-Xmin)*dx2),GetMaxY-45,
                 20+dx1+Round((beta-Xmin)*dx2),
                            GetMaxY-45-Round(ft[i]*dy1));
            OutTextXY(150,GetMaxY-20,
                   'Распределение равномерной плотности')
            End;
        '3':OutTextXY(20,GetMaxY-20,
                   'Нормальное распределение');
    End;
    Str(n,St1); Str(mx:3:3,St2); OutTextXY(106,GetMaxY-35,St2);
    Str(Xmin:3:2,St4); OutTextXY(27,GetMaxY-35,St4);
    Str(Xmax:3:2,St5); OutTextXY(177,GetMaxY-35,St5);
    While not KeyPressed do;
    { ЗАКРЫТИЕ ГРАФИЧЕСКОГО РЕЖИМА}
    CloseGraph
End
END.