Programa de La regla del Trapecio

iogha
10 de Junio del 2006
Hola, alguien me podria ayudar a hacer un programa en c++ que calcula las integrales mediante la regla del trpecio simple y multiple??
Encontre este en internet, pero la funcion esta predeterminda, se los agradeceria mucho
#include <math.h>
#include <stdio.h>
#include <iostream.h>
#include <conio.h>
float funcion(float);
float trapecio(float, float);
float trapecio_extendida(float, float, int);
void main(void);
float funcion(float x)
{
return pow(x,2)/2+3*x+2;
}

float trapecio(float a, float b)
{
return (b-a)*(funcion(a)+funcion(b))/2;
}
float trapecio_extendida(float a, float b, int n)
{
float h=(b-a)/n;
float sumatoria=0;
int j;
for (j=1;j<=n-1;j++)
sumatoria=sumatoria+funcion(a+j*h);

return (h/2)*(funcion(a)+2*sumatoria+funcion(b));
}
void main()
{

printf("\nResultado trapecio--> %f",trapecio(1,5));
printf("\nResultado trapecio extendida--> %f",trapecio_extendida(1,5,8));
scanf("\n");
}

Noel Solw
10 de Junio del 2006
A continuacion un programa que encuentra la integral usando el metodo del trapecio,
No conozco la diferencia entre trapecio simple y multiple.

// program trapecio.cpp
// integracion usando el metodo de los trapecios
// 10/6/2006 - writtem im bprland c++, ver 4.52

#include <iostream.h>
#include <iomanip.h>

double f(double x) // aqui intoduces la funcion que quieres integrar
{
return x*x;
}

int main()
{
double h = 0.25,a = 3,b = 6,sum;
sum = (f(a)+f(b))/2.0;
for(double x = a + h;x < b;x += h)
sum += f(x);
sum = sum*h;
cout << "integral = " << sum << endl << endl;
cout << "end of ptogram - good bye ! ! !" << endl;
return 0;
} // MAIN

// con los datos del programa obtendras integral = 63.0312
// la integral excata seria ; [x^3/3] entre 3 y 6 = 63.00