define f(x) {return  re( 1 + power(x,2)) ;}
 
define F(x) {return re( x + (1/3.0)*power(x,3)) ;}

define Df(x) {return  re( 2*x  );}


define listagem(inicio, fim, delta) 
{
 local teste = 0;   
 local dados = fopen("dados","w");
 local teste = 0.0;
 local x = inicio;
 while( x < fim)
 {	 
	fprintf(dados, "%f  %f %f  \n", x,f(x), Df(x));
   	x = x + delta;
 } 
 fclose(dados);
}

define transfere () 
{
 local transfere = fopen("transfere","w");
 fprintf(transfere, "set terminal postscript portrait enhanced  \n");
 fprintf(transfere, "set output 'exer07_02.eps' \n");
 fprintf(transfere, "set pointsize 0.1 \n");  
 fprintf(transfere, "plot 'dados' with points  \n");
 fclose(transfere);
}

define riemann(inicio, fim, delta) 
{
 local soma = 0;   
 local x = inicio;
 while( x < fim)
 {	 
	soma = soma + f(x);
   	x = x + delta;
 } 
 return( soma*delta);
}





## delta = 0.01; precisao = 0.1
##listagem(-3,3,1);
riemann(-3,2,0.01);