using System;

public class Program

{

    public static void Main()

    {

          Console.WriteLine("Método de Simpson");

        Console.WriteLine();

        //Se declaran las variables del algoritmo

        double deltaX, altura_1, altura_2, altura_3, altura_total, x;

        double partes=2, area =0;

        //Se declaron los limites al igual que en los métodos anteriores

        double lim_inferior = -10, lim_superior = 10;            

        // Se calcula el cambio en x es decir cuanto va incrementar basado en el numero de partes con 2 es mas que suficiente

        // por eso este método es bien machin como Sandor "The Hound" Clegane 

        deltaX = (lim_superior - lim_inferior) / partes;            

        x = lim_inferior;            

​

        while (x < lim_superior)

        {

        // Se trabaja igual que en los temas anteriores solo que aqui calcularemos 3 alturas para nuestro objetivo. 

            altura_1 = (Math.Sin(x) / x) + 1;

            x = x + deltaX/2;

            altura_2 = (Math.Sin(x) / x) + 1;

            x = x + deltaX/2;

            altura_3 = (Math.Sin(x) / x) + 1;           

           //Se obtiene la altura total en base ala proporcion de las 3 alturas

            altura_total = (altura_1 + 4*altura_2 + altura_3) / 6;            

            //Se efectua el cálculo de la sumatorias de las áreas (integral numérica)

            area = area + (deltaX) * altura_total;

        }

​

        Console.WriteLine("El valor del area (integral) es " + area);            

        Console.ReadLine();

    }

}