LosG=80;                  (*Losgröße*)
M=16;                       (*Reklamationsgrenze*)
p=0.08;                (*Wahrscheinlichkeit*)
n=13;                    (*Stichprobengröße*)
c=1;                     (*Annahmezahl*)
AnzSim=10;

Es wird ein Los mit 80 Zufallszahlen erzeugt.

los=RandomReal[{0,1},LosG]

Wenn immer die erzeugte Zufallszahl < p, wird ihr der Wert 0 zugewiesen, sonst 1

los=If[#<p,0,1]&/@RandomReal[{0,1},LosG]

Die Wahrscheinlichkeit, eine Eins zu findenist hier

Total[los]/LosG//N

Jetzt erzeuge ich n ( n = 13 ) Zufallszahlen zwischen 1 und 80, die die Nummern der aus los zu ziehenden Teile angeben:

stichprobe=RandomInteger[{1,80},n]

Diese Teile entnehme ich aus los

stichprobe1=los[[#]]&/@stichprobe

und zähle die Einsen

c1=Total[stichprobe1]

Das kann man jetzt  natürlich in einer Schleife etliche Male machen und sehen, wie sich c1 verhält.

fc1:=Module[{},
los=If[#<p,0,1]&/@RandomReal[{0,1},LosG];
stichprobe1=los[[#]]&/@RandomInteger[{1,80},n];Total[stichprobe1]
]

fc1tab=Table[fc1,{20}]

ListLinePlot[fc1tab,PlotRange->{8,14}]