Вычислением математического ожидания случайной величины в Unite

program Project1;
 
{$APPTYPE CONSOLE}
 
uses
  SysUtils,
  Windows,
  Math,
  Unit1 in 'Unit1.pas';
 
type  mu= array [1..3] of Extended;
 
var k:char;
    i,j,a,b,c:integer;
    sumA,sumB,sumC: extended;
    l,r1,r2,r3:extended;
    x:masa;
    y:masb;
    z:masc;
    g:mu;
    q:boolean;
 
begin
  repeat
    repeat
      writeln ('Как решать будем?');
      writeln ('1 - Значения с клавиатуры');
      writeln ('2 - Рандомные значения');
      readln (j);
    until (j=1) or (j=2);
 
    repeat
    writeln;
    write ('Введите размерность первого массива (не более 10): ');
    readln (a);
    until a<=10;
 
    repeat
    writeln;
    write ('Введите размерность второй массива (не более 10): ');
    readln (b);
    until b<=10;
 
    repeat
    writeln;
    write ('Введите размерность третей массива (не более 10): ');
    readln (c);
    until c<=10;
    writeln;
 
    case j of
 
      1:begin
 
        begin
        sumA:=0;
          for i := 1 to a do
            begin
            write ('masA[',i,']-');
            readln (x[i]);
            sumA:=sumA+x[i];
            end;
        end;
        writeln;
        begin
        sumB:=0;
          for i := 1 to b do
            begin
            write ('masB[',i,']-');
            readln (y[i]);
            sumB:=sumB+y[i];
            end;
        end;
        writeln;
         begin
        sumC:=0;
        for i := 1 to c do
          begin
          write ('masC[',i,']-');
          readln (z[i]);
          sumC:=sumC+z[i];
          end;
        end;
      end;
 
      2:begin
 
        begin
        sumA:=0;
        randomize;
         for i := 1 to a do
          begin
          x[i]:=random(10)-random(10);
          sumA:=sumA+x[i];
          writeln ('masA[',i,']-',x[i]:5:1);
          end;
        end;
        writeln;
        begin
        sumB:=0;
        randomize;
         for i := 1 to b do
          begin
          y[i]:=random(10)-random(10);
          sumB:=sumB+y[i];
          writeln ('masB[',i,']-',y[i]:5:1);
          end;
        end;
        writeln;
        begin
        sumC:=0;
        randomize;
         for i := 1 to c do
          begin
          z[i]:=random(10)-random(10);
          sumC:=sumC+z[i];
          writeln ('masC[',i,']-',z[i]:5:1);
          end;
        end;
      end;
    end;
 
    writeln;
    Ma(x,a);
    Mb(y,b);
    Mc(z,c);
    writeln;
 
    repeat
    q:=true;
    for  i:= 1 to 3-1 do
    if g[i] < g[i+1] then
      begin
      l:=g[i];
      g[i]:=g[i+1];
      g[i+1]:=l;
      q:=false;
      end;
    until q = true;
    writeln ('Максимальное из чисел M: ', g[1]:0:2);
 
    writeln;
    writeln ('Заново?  Y/N');
    readln(k);
  until  (upcase(k)='N');
 
end.

unit Unit1;
 
interface
 
uses
  SysUtils,
  Windows,
  Math;
type  masa=array [1..10] of  Extended;
      masb=array [1..10] of  Extended;
      masc=array [1..10] of  Extended;
      mu= array [1..3] of Extended;
 
procedure Ma(x:masa;var a:integer);
procedure Mb(y:masb;var b:integer);
procedure Mc(z:masc;var c:integer);
 
var i:integer;
    g:mu;
    q:boolean;
 
 
implementation
 
procedure Ma (x:masa;var a:integer);
var r1: extended;
    g:mu;
  begin
   r1:=0;
  for i := 1 to a do
   r1:=r1+(1/a)*x[i];
   r1:= r1 + x[i]/a;
   g[1]:=r1;
   writeln ('X (a)-[',r1:0:2,']');
  end;
 
procedure Mb(y:masb;var b:integer);
var r2:extended;
  begin
   r2:=0;
  for i := 1 to b do
   r2:=r2+(1/b)*y[i];
   r2:= r2 + y[i]/b;
   g[2]:=r2;
   writeln ('Y (b)-[',r2:0:2,']');
  end;
 
procedure Mc(z:masc;var c:integer);
var r3:extended;
  begin
   r3:=0;
  for i := 1 to c do
   r3:=r3+(1/c)*z[i];
   r3:=r3 + z[i]/c;
   g[3]:=r3;
   writeln ('Z (c)-[',r3:0:2,']');
  end;
  end.