Перевод кода с Си подобного языка - PascalABC.NET

#include < stdio.h >
#include < stdlib.h >
#define _USE_MATH_DEFINES
#include < math.h >

//*************************************************************************************************************
// аппроксимация гамма-функции в интервале от 1 до 2
// отношением полиномов 8 степени
double gammaapprox(double x){
        double p[]={-1.71618513886549492533811e+0,
                        2.47656508055759199108314e+1,
                       -3.79804256470945635097577e+2,
                        6.29331155312818442661052e+2,
                        8.66966202790413211295064e+2,
                       -3.14512729688483675254357e+4,
                       -3.61444134186911729807069e+4,
                        6.64561438202405440627855e+4};
 
        double q[]={-3.08402300119738975254353e+1,
                        3.15350626979604161529144e+2,
                       -1.01515636749021914166146e+3, 
                       -3.10777167157231109440444e+3,
                        2.25381184209801510330112e+4,
                        4.75584627752788110767815e+3,
                       -1.34659959864969306392456e+5,
                       -1.15132259675553483497211e+5};
        double z = x - 1.0;
        double a = 0.0;
        double b = 1.0;
        for(int i = 0; i<8; i++){
                a =(a+p[i])*z;
                b = b*z+q[i];
        }
        return (a/b+1.0);
}

//*************************************************************************************************************
// Гамма-функция вещественного агрумента
// возвращает значение гамма-функции аргумента z
double gamma(double z){
       if((z>0.0)&&(z<1.0))   return gamma(z+1.0)/z;        // рекурентное соотношение для 0
       if(z>2)                return (z-1)*gamma(z-1);      // рекурентное соотношение для z>2 
       if(z<=0)               return M_PI/(sin(M_PI*z)*gamma(1-z)); // рекурентное соотношение для z<=0 
       return gammaapprox(z);   // 1<=z<=2 использовать аппроксимацию
}

//*************************************************************************************************************
 
// Основная программа для рассчета значения гамма-функции вещественного аргумента
int main(){
        float z;
        printf("enter argument z: ");
        scanf(" %f",&z);
        double g = gamma((double)z);
        printf("gamma(%.2f) = %e\n", z,g);
        system("Pause");
        return 0;
}

var
  P: array of double := (-1.71618513886549492533811e+0,
                          2.47656508055759199108314e+1,
                         -3.79804256470945635097577e+2,
                          6.29331155312818442661052e+2,
                          8.66966202790413211295064e+2,
                         -3.14512729688483675254357e+4,
                         -3.61444134186911729807069e+4,
                          6.64561438202405440627855e+4);
  
  Q: array of double := (-3.08402300119738975254353e+1,
                          3.15350626979604161529144e+2,
                         -1.01515636749021914166146e+3, 
                         -3.10777167157231109440444e+3,
                          2.25381184209801510330112e+4,
                          4.75584627752788110767815e+3,
                         -1.34659959864969306392456e+5,
                         -1.15132259675553483497211e+5);
                       
function GammaApproximation(x: double): double; 
begin
  var (z, a, b) := (x - 1.0, 0.0, 1.0);
  for var i := 0 to 7 do
    (a, b) := ((a + P[i]) * z, b * z + Q[i]);
  Result := a / b + 1.0;
end;
 
function Gamma(z: double): double;
begin
  if (z > 0) and (z < 1) then Result := Gamma(z + 1.0) / z;
  if (z > 2) then Result := (z - 1) * Gamma(z - 1);
  if (z <= 0) then Result := Pi / (Sin(Pi * z) * Gamma(1 - z));
  Result := GammaApproximation(z);
end;
 
begin
  var z := ReadlnReal();
  Writeln($'Gamma({z}) = {Gamma(z)}');
end.