Taschenrechnerprogramme zum Würfelexperiment FAIR DIE mit dem Sharp EL-9600


Auf der Grundlage von M (z.B. M=300) Würfelexperimenten, wobei in jedem Experiment N (z.B. N=100) Augenzahlen simuliert werden sollen, errechnet man gemäß dem Chi-Quadrat-Anpassungstest M Chi-Quadrat-verteilte Testgrößen. Die Auswertung dieser Testgrößen im Histogramm wird mit der Chi-Quadrat-Dichte-Funktion als statistische Prüfverteilung des Chi-Quadrat-Anpassungstests verglichen.

Programmidee:
C.C.Edwards: Does a TI-8x Cast a Fair Die? in: Eightysomething!, Vol. 6, No.3, 1997, p. 9-10.
siehe Internet: http://www.comcal.net/www.ti.com/calc/docs/act83stat.htm oder als PDF-File.




Literaturhinweis:
Aulenbacher, Paditz, Wabel-Frenk: Lehr- und Übungsbuch Mathematik Bd.3 (Teil Stochastik: Beispiel 6.1 und Aufgabe 19.1),
Fachbuchverlag Leipzig 1996 ( ISBN 3-446-18715-4)

PC-Version des ehemaligen Studenten M.Jany (Imm.-jahrgang 1996): Würfelexperiment für Win95 oder für Win3.11



EL-9600: Program FAIRDIE

ClrT
ClrG
ClrDraw
ClrList L1,L2,L3,L4,L5,L6
SetList L1,L2,L3,L4,L5,L6
Print "FAIR DIE..."
Print "1 FOR YES"
Print "0 FOR NO"
Input A
Print "NUM EXPERIMENTS:"
Input M
Print "NUM ROLLS:"
Input N
N / 6  --->  E
M  --->  dim(L1)
6  --->  dim(L2)
   0  --->  K
Label LOOPA
   K+1  --->  K
   0 * L2  --->  L2
       0  --->  I
Label LOOPB
      I+1  --->  I
      If A  = 0 Goto CCC
         ipart ( 6 * random + 1 )  --->  D
         Goto CC
    Label CCC
      ipart (11 * random )  --->  R
      R - 6 * ipart ( R / 6 ) + 1  --->  D
    Label CC
    L2(D) + 1  --->  L2(D)
If I< N Goto LOOPB
sum( (L2-E)^2 / E )  --->  C
Print "EXP. NO:"
Print K
Print "CHI^2 ="
Print C
round(C,3)  --->  L1(K)
If K< M Goto LOOPA
Print "NUM EXPERIMENTS:"
Print M
Print "NUM ROLLS:"
Print N
Print "CODE FAIRDIE:"
Print A
ClrList L2
sortA(L1)


EL-9600: Program PRIMFREQ

Print "A MOMENT PLEASE, BUSY ..."
M  --->   dim(L2)
M  --->   dim(L3)
L1(1)  --->   L2(1)
1  --->   L3(1)
   1  --->   K
     1  --->   I
Label PPP
   K+1  --->   K
   If L1(K) = L1(K-1) Goto PP1
     I+1  --->   I
     L1(K)  --->   L2(I)
     1  --->   L3(I)
   Goto PP2
Label PP1
     L3(I)+1  --->   L3(I)
Label PP2
If K< M Goto PPP
I  --->   dim(L2)
I  --->   dim(L3)




EL-9600: Program SECUFREQ

Print "A MOMENT PLEASE, BUSY ..."
ClrList L4,L5
SetList L4,L5
ipart ( L2(dim(L2))+2 )  --->   dim(L4)
ipart ( L2(dim(L2))+2 )  --->   dim(L5)
seq( -.5+X, 0, ipart (L2(dim(L2))+1) )  --->   L4
0  --->   L4(1)
0 * L4  --->   L5
   0  --->   K
Label QQQ
   K+1  --->   K
   ipart (L1(K)) + 2  --->   L
   L5(L) + 1  --->   L5(L)
If K< M Goto QQQ
L5(2)  --->   L5(1)
Bilder zum Würfelexperiment ( idealer Würfel oder nichtidealer Würfel).

Im Y= - Menü theoretische Dichte Chi25 programmieren: Y1 = M * pdfX2 ( X, 5 )

Im WINDOW - Menü (Rect) Grafik-Fenster Xmin= 0, Xmax=15, Xscl=1, Ymin= -8, Ymax=55, Yscl=10 einstellen.

Im STAT PLOT - Menü PLOT1 (Histogramm) und PLOT2 (Häufigkeitspolygon) einstellen.

Ludwig Paditz, 10. April 2002