Сформировать и вывести массив Х из 50 вещественных чисел в пределах от 20.0 до 40.0.Вычислить среднее арифметическое Х и среднеквадратическое отклонение s — Pascal(Паскаль)

Program P3;

var
  x: array [1 .. 50] of real;
  a, xm, sum, s, ks, ds: real;
  i, kk, dd: integer;

begin
  for i := 1 to 50 do
    x[i] := random(20) + 20;
  xm := 0;
  kk := 0;
  ks := 0;
  for i := 1 to 50 do
    write(x[i], ' ');
  writeln;
  for i := 1 to 50 do
    sum := sum + x[i];
  xm := sum / 50;
  writeln('Srednee arifmeticheskoe X =', xm:0:2);
  s := (sqr(sum) / 50) - sqr(xm); { формула дисперсии }
  s := sqrt(s);
  s := sqrt(s);
  { среднеквадратическое отклонение s равно квадратному корню от дисперсии }
  writeln('Srednekvadraticheskoe otklonenie s =', s:0:2);

  for i := 1 to 50 do
    if x[i] > (xm + s) then
    begin
      inc(kk);
      ks := ks + x[i];
    end;
  writeln('Elementov massiva, prevyshayuschih X+s =', kk);
  writeln('Ix summa =', ks);

  for i := 1 to 50 do
    if x[i] < (xm - s) then
    begin
      inc(dd);
      ds := ds + x[i];
    end;
  writeln('Elementov massiva, prevyshayuschih X+s =', dd);
  writeln('Ix summa =', ds);
  readln;

end.