Вычисление определенного интеграла методом парабол (Симпсона)

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
 
namespace ConsoleApplication33
{
    class Program
    {
        static void Main(string[] args)
        {
            double a = 1, b = 2, h;
            Console.WriteLine("Вычисление определенного интеграла методом парабол (Симпсона)");
            Console.Write("Введите четное количество ординат разбиения n/2=");
            int n = Convert.ToInt32(Console.ReadLine());
            Console.WriteLine("Получившиеся промежутки разбиения");
            double[] promX = new double[2 * n + 1];
            for (int i = 0; i <= 2 * n; i++)
            {
                h = (b - a) / (2 * n);
                promX[i] = a + i * h;
                if (i == 0) Console.Write("a(" + a);
                if ((i != 0) && (i != 2 * n)) Console.Write(")___(" + promX[i]);
                if (i == 2 * n) Console.WriteLine(")___(" + b + ")b");
            }
            double fa, fb, f2 = 0, f1 = 0;
            fa = a * Math.Pow(2 - a, 1 / 2);
            fb = b * Math.Pow(2 - b, 1 / 2);
            for (int i = 1; i <= 2 * n - 1; i += 2)
            {
                f1 += (promX[i] * Math.Pow(2 - promX[i], 1 / 2));
            }
            for (int i = 2; i <= 2 * n - 2; i += 2)
            {
                f2 += (promX[i] * Math.Pow(2 - promX[i], 1 / 2));
            }
 
            double integral = ((b - a) / (6 * n) * ((fa + fb) + 4 * f1 + 2 * f2));
            Console.WriteLine("Результат: " + integral);
            Console.ReadKey();
        }
    }
}