Вычисление определенного интеграла методом криволинейных трапеций

 static double Funciton(double x)
        {
            return x * Math.Sqrt(1 + x);
        }
 
        static double Summa(double n,double x,double step)
        {
            double summa= 0;
            for (int i = 1; i < n; i++)
            {
                summa += Funciton(x);
                x += step;
            }
 
            return 2*summa;
        }
 
        static void Main(string[] args)
        {
            double countStep = 1000;
            double leftBorder = 0;
            double rightBorder = 3;
            double step = (rightBorder - leftBorder) / countStep;
 
            double x=leftBorder+step;
 
            double result=(step/2.0)*(Funciton(leftBorder)+Summa(countStep,x+step,step)+Funciton(rightBorder));
 
            Console.WriteLine(result);
 
 
            Console.ReadKey();
        }

      static double f(double x)
        {
            return x * Math.Sqrt(1 + x);
 
        }
        static double Integral(double x )
        {
            double a = 1;
            double b = 1;
            return (2*(3*a*x-2*b)*Math.Pow((a*x+b),3.0/2))/15;
 
        }
        public static double Trapezoidal(double a, double b, int n)
        {
            double sum = 0.0;
            double h = (b - a) / n;
            for (int i = 0; i < n; i++)
            {
                sum += 0.5 * h * (f(a + i * h) + f(a + (i + 1) * h));
            }
            return sum;
        }
        static void Main(string[] args)
        {
            int n = 1000;
            double result, analitic;
            double a = .0;
            double b = 3.0;
 
            result = Trapezoidal(a, b, n);
            analitic = Integral(3) - Integral(0);
            Console.WriteLine("Numerical Result =  " + result.ToString());
            Console.WriteLine("Analitical Result = " + analitic.ToString());
            Console.ReadKey();
        }