Сгенерировать две последовательности

unit Unit1;
interface
 
uses
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, StdCtrls, TeEngine, Series, ExtCtrls, TeeProcs, Chart;
 
type
TForm1 = class(TForm)
    Memo1: TMemo;
    gg: TChart;
    Series1: TBarSeries;
    Button1: TButton;
    Label1: TLabel;
    procedure Button1Click(Sender: TObject);
 
private
{ Private declarations }
public
{ Public declarations }
    end;
    var
    Form1:TForm1;
    mas1,mas2 : array [1..50] of real;
    fff,tr:real;
    i:integer;
implementation
 
{$R *.dfm}
procedure gen;
begin
Form1.Memo1.Lines.Add('------------------------------------------');
for i:=1 to 50 do
begin
mas2[1]:=mas1[1];
mas2[i]:=(mas2[i-1]+mas1[i]);
end;
for i:=1 to 50 do
begin
Form1.Memo1.Lines.Add(IntToStr(i)+' -- '+FloatToStr(mas2[i]));
tr:=tr+mas2[i];
end;
tr:=tr/50;
end;
procedure gistogr;
var
i, j:integer;
max, min,zh, h, interval: real;
x, y:array [0..6] of real;
 
begin
Form1.gg.Series[0].Clear;
max:=mas2[1]; min:=mas2[1]; h:=0;
for i:=1 to 50 do begin
if (mas2[i] > max) then max:= mas2[i];
end;
for i:=1 to 50 do begin
if (mas2[i] < min) then min:= mas2[i];
end;
interval:=max-min;
h:=interval/10;
zh:=h;
x[0]:=0;
for i:=1 to 10 do
x[i]:=x[i-1]+h;
for i:=1 to 10 do
begin
for j:=1 to 50 do
begin
if (mas2[j]>x[i-1]) and (mas2[j]<=x[i]) then
y[i]:=y[i]+1;
end; zh:=h+zh;
end;
 
for i:=1 to 10 do begin
Form1.gg.Series[0].AddXY(x[i], y[i]);
end;
 
end;
procedure TForm1.Button1Click(Sender: TObject);
begin
gen;
gistogr;
label1.Caption:='Ñðåäíåå çíà÷åíèå ' + floattostr(tr);
   end;
end.

unit Unit1;
interface
 
uses
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, StdCtrls, TeEngine, Series, ExtCtrls, TeeProcs, Chart;
 
type
TForm1 = class(TForm)
    Chart1: TChart;
    Series1: TBarSeries;
    Button1: TButton;
    Label1: TLabel;
    Label2: TLabel;
 
    procedure Button1Click(Sender: TObject);
 
private
{ Private declarations }
public
{ Public declarations }
    end;
    var
    Form1:TForm1;
 
implementation
 
{$R *.dfm}
procedure TForm1.Button1Click(Sender: TObject);
var
x,r:array [1..30] of real;
counts: array [1..10] of Integer;
i,j:integer;
lambda,min,max,sz,mx,width1,st:real;
begin
randomize;
 
lambda:=0.8;
For i:= 1 to 30 do
begin
x[i]:=-ln(random)/lambda;
 
end;
 
min:=x[1];
max:=x[1];
for i:= 1 to 29 do
begin
for j:=1 to 29 do
begin
if (x[j]>x[j+1]) then
begin
min:=x[j];
x[j]:=x[j+1];
x[j+1]:=min
end;
 
end;
end;
 
for i := 1 to 30-1 do
begin
r[i]:=x[i+1]-x[i];
end;
 
//ìàò îæèäàíèå
for i := 1 to 29 do
mx:=mx+r[i];
mx:=(1/29)*mx;
//ñð çíà÷åíèå
for i := 1 to 29 do
sz:=sz+x[i];
sz:=sz/29;
Label2.Caption:='';
label2.Caption:='ñðåäíåå çíà÷åíèå = '+floattostr(sz);
 
min:=x[1];
max:=x[1];
for i:= 2 to 29 do
begin
if (r[i]<min) then min:=r[i];
if (r[i]>max) then max:=r[i];
end;
width1:=(max-min)/10;
st:=Min;
for j:= 1 to 10 do
begin
counts[j]:=0;
for i:=1 to 29 do
if (r[i]>=St) and (r[i]<St+Width1) then
Inc(counts[j]);
St:=St+Width1;
end;
if r[29]=Max then
Inc(counts[10]);
Chart1.SeriesList[0].Clear;
for i:= 1 to 10 do
Chart1.SeriesList[0].AddXY(i,counts[i]/29);
 
end;
 
end.