{***  programa para calcular integrais com somas de Riemann
 ***  autores Geronimo/Praciano-Pereira data 26/01/92    **
 ***  programa testado - sem erros encontrados           **}

program Riem02;

{Uses Crt;}

Var contador  : Integer;
    delta,x, somaanterior, soma, inicio, fim, deltax : Real;

Function f(x:Real):Real;
Begin
   f := sin(x);
End;

Begin
  {ClrScr}; soma := 0; contador := 1;
  Writeln('Inicio e fim do intervalo de integracao');
  Write('Inicio: ');Readln(inicio);
  Write('Fim: ');Readln(fim);
  Write('deltax: 0.001  < deltax  < 0.1  ');Readln(deltax);
  x := inicio; {ClrScr};
  Repeat
    somaanterior := soma; soma := 0;
    x := inicio;
    While x < fim Do
      Begin
        soma := soma + f(x)+f(x+deltax);
        x := x + deltax;
        {Writeln('S=  ', soma*deltax:9:4);}
      End;
    contador := contador + 1;
    soma := 0.5*soma*deltax;
    WriteLn('Experiencia no. ',contador);
    WriteLn('Valor da Integral: ',soma:9:4);
    WriteLn('deseja repetir os calculos: forneca o deltax,');
    Write('(deltax = 0 para o programa), deltax:  ');
    ReadLn(deltax);
  Until deltax = 0;
  {ClrScr};
End.