##
## wuerfeln.m 
##   Octave-Script zur Demonstration der Berechnung der Verteilung
##   der Augensumme beim n-fachen Würfeln
##
## Erstellt für die Veranstaltung "Statistik" am Fachbereich
## "Elektrotechnik und Informatik" der Hochschule Niederrhein, Krefeld
##
## Autor:   Christoph Dalitz
## Datum:   20.10.2008
## Lizenz:  frei kopierbar unter der GNU General Public License Version 2
##          (siehe http://www.gnu.org/licenses/)
## Voraussetzung:
##          Octave (frei erhältlich von http://www.octave.org/)
##          ab Version 2.9


n = 3; # Anzahl Würfe

p = (1/6) * ones(1,6);
q = [1];
for i = 1:n
  q = conv(q,p);
endfor

# für Octave < 2.9:
#titlestr = sprintf("^;%d-fache Augensumme;", n);
#plot(n:6*n, q, titlestr);
# ab Octave 2.9:
stem(n:6*n, q);
legend(sprintf("%d-fache Augensumme;", n));
axis([n-0.5,6*n+0.5]);