#include <iostream>
#include <stdlib.h>
#include <windows.h>
#include <math.h>
using namespace std;

double f(double x)
{
	return (1+log(x)-exp(x))/(x*x+3);  // expression de la fonction (log représente ici ln)
}
int main(int argc, char* argv[])
{
double a,b,dt,g,h,k;
const int K= 100;						// nombre de parabole * 2
double resultat= 0.0;					// résultat approche
b =2;								    // borne superieure
a =1;									// borne inferieure
dt=(b-a)/K;								// notre intervalle entre 2 points
short int stwa,stout;
stwa =GetTickCount();
for(int n = 1; n<= (K/2); n++)
	{
	k= a+ 2*n*dt;
	g= k- 2*dt;
	h= (g+k)/2;
	resultat += (f(g)+4*f(h)+f(k));
	}
resultat *= (dt/3);
stout=GetTickCount();
stout = stout-stwa;
if(!stout)
	cout <<"Temps necessaire en millisecondes pour le calcul inferieur a 0 ms " << endl;
else
	cout <<"Temps necessaire en millisecondes pour le calcul: " << stout << "ms" << endl;
cout <<"Resultat approche: "<< resultat;
cout << endl;
system("PAUSE");
return 0;
}