STABW berechnen-Qualität des Zufallszahlengenerators beurteilen

[mabstudent]

Hallo, ich verzweifle an einer Teilaufgabe und wäre für ein paar Tipps sehr dankbar.
Die Gesamtaufgabe besteht darin den Zufallszahlengenerator der in Delphi/Lazarus über random(bereich) aufgerufen werden kann auf seine Qualität zu untersuchen und grafisch zu visualisieren. Am Ende soll die Standardabweichung berechnet werden und dabei hakt es. Zum besseren Verständnis poste ich mal den kompletten Quelltext und ein Bild des laufenden Progs im Anhang.
x Achse = Schleifenindex
y Achse = Zufallszahl

Code:

procedure TForm1.Button1Click(Sender: TObject);
var i,bereich,r,summe:integer;
    y:longint;
    zaehler: array[0..600]of longint;
    mittelwert:double;

begin
  randomize;
  image1.Canvas.brush.color:=clwhite;
  image1.canvas.fillrect(rect(0,0,600,200));
  bereich:=strtoint(edit1.text); //Zufallszahlenbereich

  IF checkbox1.checked=true THEN //von unten nach oben auftragen
  begin

  FOR i:=0 to image1.width-1 DO
    begin
    zaehler[i]:=0;         //alle Zähler i..imagebreite nullen
    end;

  FOR i:=0 TO image1.width-1 DO
    begin
    r:=random(bereich);
    inc(zaehler[r]);      //inkrementiert den Zählerplatz i um die Zufalsszahl r
    IF checkbox3.checked=true THEN image1.canvas.pen.color:=clred;
    IF checkbox4.checked=true THEN image1.canvas.pen.color:=clgreen;
    IF checkbox5.checked=true THEN image1.canvas.pen.color:=clblue;
    image1.canvas.moveto(i,200);   //nach unten links fahren
    image1.canvas.lineto(i,200-r); //Linie ziehen zur Zufallszahl
    end;
  end;
  end;
label3.caption:=floattostr(mittelwert);

Frage: Wie kann ich die STABW, also [STABW=sqrt(summe(jede Zufallszahl-mittelwert aller)/Anzahl von Zufallszahlen-1))] berechnen? -d.h.wie den Mittelwert, wie die Summe aller Zufallszahlen,...
Sollte ich das in die vorhandene Schleife integrieren?
Vorab schonmal vielen Dank und entschuldigt den Quellcode(denke nach 1Semester Maschinenbauinformatik noch nicht wirklich informatisch )