Из круглой, доверху наполненной жидкостью емкости (высота 10м) выливают ее содержимое через отверстие в днище площадью 0.1 кв.м. Скорость вытекания V=sqrt(2gb),где g-постоянная сила тяжести, b-высота столба жидкости. Рассчитать изменение уровня жидкости в емкости в зависимости от времени до полного вытекания жидкости — Pascal(Паскаль)

const
  G = 9.8;
  S = 0.1;
  D = 2;
  dt = 0.5;
 
var
  i: integer;
  diam: double;
 
  t, h: double;
  dV: double;
begin
  diam := 2 * sqrt(S / pi);
 
  t := 0.0;
  h := 10;
  i := 0;
 
  while h > 0 do begin
    dV := (pi*sqr(diam) / 4) * sqrt(2 * h * G) * dt;
    h := h - 4 * dV / (pi * sqr(D));
    if h < 0 then h := 0;
 
    t := t + dt;
    inc(i);
 
    writeln('t = ', t:10:5, ' h = ', h:10:5);
 
    if i mod 20 = 0 then begin
      writeln('press <enter> to continue...');
      readln;
    end;
  end;
  readln;
end.