/************************************************/
/*   Programa Integracion_Regla_1_3_Simpson;    */
/*										                  */
/*	     	{ Dada una funcion }                   */
/************************************************/

#include 
#include 
#include 

int Lee_Datos(void);
double Ecuacion(double Xe);

  int      Nseg;
  float    a,b;
  double   Xi;

/**************   Programa Principal     *********/

int main (void)
{
  int      i;
  float    Base;
  double   Area;
  double   Sumpar = 0;
  double   Sumimpar = 0;

  clrscr();
  printf("\n     METODO DE INTEGRACION POR 1/3 de SIMPSON ");
  printf("\n==================================================");
  printf("\n         Dada una Funcion cualesquiera");
  printf("\n");
  Lee_Datos();

  Base = (b-a)/Nseg;
  Xi   = a;

/******** Encabezados ************/

  printf("\nDatos Tabulados.......");
  printf("\n-------------------------");
  printf("\n|   i |    Xi  | Funcion");
  printf("\n-------------------------\n");
  printf("|   0 | %.2f   |  %.4lf",a,Ecuacion(a));

/****** Inicia Proceso (Calculo de Sumatorias) **********/

  for ( i=1; i %.8lf",Area);
  getche();

}

int Lee_Datos(void)

{
	printf("\n  Numero de Segmentos (PAR) =");
	scanf("%d",&Nseg);
	printf("\n  Valor de a =>");
	scanf("%f",&a);
	printf("\n  Valor de b =>");
	scanf("%f",&b);
}

double Ecuacion(double Xe)

{
	   return Xe*Xe;
}


    Source: geocities.com/siliconvalley/pines/7894/metodos

               ( geocities.com/siliconvalley/pines/7894)                   ( geocities.com/siliconvalley/pines)                   ( geocities.com/siliconvalley)