\START92\
\COMMENT=PROGRAM DATE 04. FEBR. 2001
\NAME=blumes89
\FILE=BLUMES89.9XP
(t)
Prgm
ClrIO
FnOff 
setGraph("6","1")
ClrHome
ClrDraw
PlotsOff 
0\->\xmin
10\->\xmax
10\->\xscl
0\->\ymin
8\->\ymax
8\->\yscl
Disp "START"
Disp "ANZAHL DER BLUMEN ="
Disp t
Pause 
Disp "A MOMENT, PLEASE"
seq(x,x,1,81)\->\list1
seq(0,x,1,81)\->\list2
t\->\list1[81]
0\->\m
Lbl lm
m+1\->\m
10*tistat.rand83()\->\a
8*tistat.rand83()\->\b
PtOn a,b
10*iPart(b)+iPart(a)+1\->\n
list2[n]+1\->\list2[n]
If m<t
Goto lm
Pause 
0\->\n
Lbl lp
n+1\->\n
Line 0,n,10,n
If n<8
Goto lp
0\->\n
Lbl lq
n+1\->\n
Line n,0,n,8
If n<10
Goto lq
Pause 
max(list2)+1\->\p
seq(x-1,x,1,p)\->\list3
seq(0,x,1,p)\->\list4
0\->\n
Lbl ls
n+1\->\n
list2[n]+1\->\q
list4[q]+1\->\list4[q]
If n<80
Goto ls
max(list4)+1\->\ymax
5\->\yscl
max(list3)+1\->\xmax
1\->\xscl
\(-)\.5\->\xmin
NewPlot 1,4,list3,,list4,,,,1
PlotsOn 1
DispG
Pause 
list1[81]/80\->\q
sum((list3-q)^2*list4)/79\->\s
\root\(s)\->\s
Define y1(x)=80*tistat.normpdf(x,q,s) 
FnOff 1
max(list3)\->\n
q/n\->\p
list3\->\list5
80*tistat.binompdf(n,p)\->\list5
NewPlot 2,2,list3,list5,,,,4,1
PlotsOn 2
DispG
Pause 
list3\->\list6
80*tistat.poisspdf(q,list3)\->\list6
NewPlot 3,2,list3,list6,,,,1,1
PlotsOn 3
DispG
Pause 
PlotsOff 2,3
Style 1,"Thick"
FnOn 1
DispG
Pause 
PlotsOn 2,3
DispG
Pause 
EndPrgm

\STOP92\