Объём воды в неполной цистерне

var
  l, d, r, h, s, v : real;
begin
  repeat
    write('введите длину цистерны (в метрах): ');
    readln(l);
  until (l > 0) and (l <= 100);
  repeat 
    write('введите диаметр цистерны (в метрах): ');
    readln(d)
  until (d > 0) and (d <= l / 2);  
  repeat 
    write('введите уровень воды (в метрах): ');
    readln(h)
  until (h >= 0) and (h <= d);
  r := d / 2;
  {площадь сегмента круга в разрезе цистерны занятая водой}
  s := (h - r) * sqrt(2 * h * r - sqr(h)) + sqr(r) *
       (pi / 2 + arctan((h - r) / (sqrt(2 * h * r - sqr(h)))));  
  {объём}
  v := s * l;
  writeln('объем воды: ', v:6:3, ' кубических метров');
  readln
end.