/* programma atto alla risoluzione di un integrale
con approssimazione def. in ingresso */

#include <iostream.h>
#include <math.h>
float f(float x);

void main()
{
 //Intervallo di integrazione
 float minimo,massimo,passo,err,eps;
 //Numero di rettangolini
 int num;
 // Risultato finale;
 float res,prev_res;
 int i,rip;
 do
  {
   cout<<"\n\nCALCOLO DI UN INTEGRALE DEFINITO DELLA FUNZIONE f(x)=x \n"<<endl;
   cout<<"Inserire il minimo: ";cin>>minimo;
   cout<<"Inserire il massimo: ";cin>>massimo;
   cout<<"Inserire approssimazione accettabile: ";cin>>eps;

	num=1;  //num e' il numero di intervalli
	res=-1;
	do
	{
		passo=(massimo-minimo)/num;
		prev_res=res;
		res=0;
      //cout<<"\nCalcolo progressivo dell'integrale :"<<endl;
		for(i=0;i<num;i++)
		{
			res=res+f(minimo+i*passo)*passo;
       //cout<<"passo attuale: "<<(minimo+i*passo)<<" e risultato: "<<res<<endl;
		}
		//cout<<"Il risultato temporaneo e': "<<res<<endl;
		num=num*2;
      //calcolo dell'errore
		if(res>prev_res)
			err=res-prev_res;
		else
			err=prev_res-res;
	}
   while(err>eps);//confronto dell'errore con l'approssimazione data in ingresso
   cout<<"\nIl risultato finale e': "<<res<<endl;
   cout<<"\nDigita 0 per uscire o un altro numero per continuare : ";
   cin>>rip;
 }
while(rip!=0);
}


float f(float x)
{
	float tmp;  //tmp é la funzione
	tmp=x; //invece di x si potrebbe mettere,ad esempio x*x+x
	return tmp;
}