Значение интеграла вычислить приближённо по итерационной формуле левых прямоугольников

var
  a,b,x,y,h,r,c,d,e,s,s1 : real;
  i,n : longint;
  bool : boolean;
 
  Function F(x:real):real;
  begin
    if x>2*r then begin F:=c; exit; end;
    if x>0 then begin F:=-sqrt(sqr(R)-sqr(x-r)); exit; end;
    if x<-2*r then begin F:=-d; exit; end;
    if x<=0 then F:=sqrt(sqr(R)-sqr(x+r));
  end;
 
Begin
  repeat
    write('Input r,c,d : '); readln(r,c,d);
    bool:=(r>0)and(c>0)and(d>0);
    if not bool then writeln('wse poloviteln. ese raz.');
  until bool;
  repeat
    write('Input a,b : '); readln(a,b);
    bool:=b>a;
    if not bool then writeln('b dolvno bytx bolse. ese raz.');
  until bool;
  repeat
    write('Input e : '); readln(e);
    bool:=(e>0.00009)and(e<0.1);
    if not bool then writeln('e=[0.0001 .. 0.9] ese raz.');
  until bool;
 
  n:=2; s:=0;
  repeat
    s1:=s;
    h:=(b-a)/n; s:=0; x:=a;
    for i:=1 to n do begin
      s:=s+f(x); x:=x+h;
    end;
    s:=s*h;
  {  writeln('n=',n:8,'  I= ',s:0:6); }
    n:=n*2;
  until abs(s1-s)<e;
  writeln('diapason ',a:0:2,' po ',b:0:2,'   I= ',s:0:6);
  readln;
End.