Расчет Величины Гамма по закону распределения СВ Фишера

public void Funcf_F()
        {
            //
            //Этот метод предназначен для рассчета величины Гамма с параметрами
            //К1 и К2
            //
            arG_a = new double[3];
            arG = new double[3];
            //Параметры Гаммы, та величина что в скобке в формуле
            arG_a[0] = (K1 + K2) / 2; //G1_a
            arG_a[1] = K1 / 2; //G2_a
            arG_a[2] = K2 / 2; //G3_a
            //Унивирсальный метод рассчета Гаммы
            for (int i = 0; i < 3; i++)
            {
                arG[i] = 1;
                int G_shag = 1; // Изначальный шаг фактариала
                double G_fakt = 1; //Величина фактариала
                if ((arG_a[i] - Math.Floor(arG_a[i])) == 0)
                    //услови, что если крато двум
                {
                    for (int G_i = 1; G_i < (int)Math.Floor(arG_a[i]); G_i++)
                    {
                        arG[i] *= G_i;
                    }
                }
                else
                {
                    //если нет, то:
                    for (int G_i = 1; G_i < (int)(2 * Math.Floor(arG_a[i])); G_i++)
                    {
                        G_shag += 2;
                        G_fakt *= G_shag;
                    }
                    arG[i] = (Math.Sqrt(Math.PI) * G_fakt) / Math.Pow(2, Math.Floor(arG_a[i]));
                }
            }
        }

arG_a[0] = (double)(K1 + K2) / 2

for (int G_i = 1; G_i < (int)(2 * Math.Floor(arG_a[i])); G_i++) {}

for (; G_shag < (int)(2 * Math.Floor(arG_a[i])); G_shag+=2)  G_fakt *= G_shag;

for (; G_shag < (int)(2 * Math.Floor(arG_a[i])); G_shag+=2)  G_fakt *= G_shag;

for (int G_i = 0; G_i < (int)Math.Floor(arG_a[i]); G_fakt *= (1+2 * G_i++));

                    for (; G_shag < (int)(2 * Math.Floor(arG_a[i]) - 1); G_shag += 2)
                    {
                        G_fakt *= G_shag;
                    }

            //Расчет теоретических велечин
            MatOT_F = K2 / (K2 - 2); //Мат ожидание по Фишеру
            DispT_F = 2 * (Math.Pow(K2, 2) * (K1 + K2 - 2)) / (K1 * Math.Pow((K2 - 2), 2) * (K2 - 4)); //Дисперсия по фишеру
            SkoT_F = Math.Sqrt(DispT_F);

for (int G_i = 0; G_i < (int)Math.Floor(arG_a[i]); G_fakt *= (1+2 * G_i++))
;

for (; G_shag < (int)(2 * Math.Floor(arG_a[i])); G_shag+=2)
    G_fakt *= G_shag;

int K1 = 3;
int K2 = 4;
double res = (K1+K2) / 2; // res == 3
res = (K1+K2) / 2.0; // res == 3.5
res = (double)(K1+K2) / 2; // res == 3.5

class MyCoolClassName 
{
    public void PrintHelloWorld() {Console.WriteLine("Hello World!");}
    // Свойства тоже должны обозначаться также:
    public bool MyBoolProperty { get; set; }
}

int counterInForLoop;

int get_value();

for (int a = 0; a < fF._Ntb; a++)
                {
                    this.datResult.Rows[a].Cells[0].Value = fF._arTable[a, 0]; //N
                    this.datResult.Rows[a].Cells[1].Value = Math.Round(fF._arTable[a, 1], 4); //интервалы
                    this.datResult.Rows[a].Cells[2].Value = fF._arTable[a, 2]; //попадания
                    this.datResult.Rows[a].Cells[3].Value = Math.Round(fF._arTable[a, 3], 2); //теор попадания
                    this.datResult.Rows[a].Cells[4].Value = Math.Round(fF._arTable[a, 4], 4); //хи2
                    this.datResult.Rows[a].Cells[5].Value = Math.Round(fF._arTable[a, 5], 4); //частота
                    this.datResult.Rows[a].Cells[6].Value = Math.Round(fF._arTable[a, 6], 4); //вероятность
}