Интегрирование методом Симпсона. Обобщение

static double f(double x) 
        {
            return x*x+2*x; 
        }
        
        static double Simpson(double n)      //Интегрирование по методу Симпсона
        {
            double x, a, b, h, s;
            a = 0;
            b = 0.5;
            h = (b - a) / n;
            s = 0;
            x = a + h;
            while (x < b)
            {
                s = s + 4 * f(x);
                x = x + h;
                s = s + 2 * f(x);
                x = x + h;
            }
            s = h / 3 * (s + f(a) - f(b));
 
            Console.WriteLine(" Интеграл(Метод Симпсона) = " + s);
            double Simpson = s;
            return Simpson;
        }
 
        static void Main(string[] args)
        {
            
        }
    }
}

       
             static double f(double x)
        {
            return x * x + 2 * x;
        }
        static double Simpson(Func<double,double> func,double n)//Интегрирование по методу Симпсона
        {
            double x, a, b, h, s;
            a = 0;
            b = 0.5;
            h = (b - a) / n;
            s = 0;
            x = a + h;
            while (x < b)
            {
                s = s + 4 * func(x);
                x = x + h;
                s = s + 2 * func(x);
                x = x + h;
            }
            return h / 3 * (s + func(a) - func(b));
        }
        static void Main(string[] args)
        {
            Console.WriteLine(" Интеграл(Метод Симпсона) = {0}", Simpson(f, 10));
            Console.WriteLine(" Интеграл(Метод Симпсона) = {0}", Simpson(x => 3 * x + 4, 10));
            Console.ReadKey(true);
        }

using System;
 
namespace ConsoleApplication22
{
    class Program
    {
        static void Main()
        {
            Func<double, double> f = x => Math.Exp(-x*x);
            var res = Simpson(f, 0, 1);
            Console.WriteLine(res);
            var resinf = InfSimpson(f, 0, 0.0001);
            Console.WriteLine(resinf);
            Console.ReadKey();
        }
 
        static double Simpson(Func<double, double> f, double a, double b)
        {
            return (f(a) + 4*f((a + b)/2) + f(b))*(b - a)/6;
        }
 
        static double InfSimpson(Func<double, double> f, double from, double epsilon)
        {
            double result = 0, summand = double.PositiveInfinity;
            for (double a = from, b = a + 1; summand > epsilon && !double.IsInfinity(result); a = b, b *= 2)
            {
                summand = Simpson(f, a, b);
                result += summand;
            }
            return Math.Abs(result);
        }
    }
}