Решение определенных интегралов методом трапеций

Console.WriteLine(3 * Math.Sqrt(1 + 3) / 2.0 * 3.0);

using System;
 
namespace ConsoleApplication1
{
    class Program
    {
        static double Calculate(Func<double, double> func, double a, double b, double step)
        {
            double res = 0;
            while (a < b)
            {
                res += func(a) + func(a = a + step);
            }
            return res * step / 2;
        }
 
        static void Main(string[] args)
        {
            Console.WriteLine(Calculate((x) => x * Math.Sqrt(1 + x), 0, 3, 0.00001));
            Console.ReadKey();
        }
    }
}

        static void Main(string[] args)
        {
            double a = 0;           //начало
            double b = 3;           //конец
            double step = 0.0001;   //шаг интегрирования
            double result = 0;      //результат
 
            while (a < b) //пока не достигли конца
            {
                double f1 = a * Math.Sqrt(1 + a);  //значение функции в начале отрезка
                a += step;                         //сдвигаемся
                double f2 = a * Math.Sqrt(1 + a);  //значение функции в конце отрезка
                double area = (f1 + f2) / 2 * step;//находим площадь трапеции
                result += area;                    //прибавляем полученную площадь к результату
            }
 
            Console.WriteLine(result);
            Console.ReadKey();
        }

func(a = a + step);

func(a += step);

        static double Calculate(Func<double, double> func, double a, double b, double step)
        {
            double res = 0;
            for (double x = a; x < b; x += step)
            {
                res += func(x) + func(x + step);
            }
            return res * step / 2;
        }