#include <iostream>
#include <cmath>

using namespace std;

//declaration de la fonction a integrer          
double fonc(double);

int main() {
  //saisi des parametres  
  double min, max;    
  cout << "Entrez un interval d'integration : " << endl;
  cout << "min = ";   cin >> min;
  cout << "max = ";   cin >> max;
  int div;
  cout << "Entrez un nombre de suos-divisions :" << endl;
  cout << "N = ";   cin >> div;
  double deltax = (max-min)/div;
  
  double intTrapeze = 0.;
  double intSimpson = 0.;
  for (int i=0; i<=div; i++) {
  	//methode du trapeze
    intTrapeze += (fonc(min+i*deltax) + 
                   fonc(min+(i+1)*deltax));
    //methode de Simpson
    intSimpson += (fonc(min+deltax*i) + 
                   4.*fonc(min+(i+0.5)*deltax) +
                   fonc(min+(i+1)*deltax));
  }                                        
  intTrapeze *= deltax/2.;
  intSimpson *= deltax/6.;
  cout << "Integrale avec la methode du trapeze: " << intTrapeze << endl;
  cout << "Integrale avec la methode de Simpson: " << intSimpson << endl;

  return 0;
}

//definition de la fonction a integrer 
double fonc(double x) {
  //dans cet example nous avons choisi un polynome de 2eme degree;
  double a, b, c;
  a = 10.;
  b = -3.;
  c = 7.;

  double y = a*(x*x) + b*x + c;
  
  return y;
}