procedure AVEVAR(DATA:array of real; N:integer; var AVE, VAR1:real);
var
    J:integer;  S:real;
begin
    AVE:=0; 
    VAR1:=0;
    For J:=1 To N do
        AVE:=AVE + DATA[J];
    AVE:=AVE / N;
    For J:=1 To N do
    begin
        S:=DATA[J] - AVE;
        VAR1:=VAR1 + S * S;
    end;
    VAR1:=VAR1 / (N - 1);
end;