'Written by Marsha Bowen as assignment #1
'Printing routines edited by WB
'declare function and integer variables
DECLARE FUNCTION CountEq! (x1!, x2!, x3!, x4!, x5!)
'prepare screen
CLS
PRINT "       This program simulates roll of five dice a 1000 times"
PRINT "       and  estimates probability   of  our-of-a-kind  event"
PRINT "        This  will be repeated 10 times to see variability"
PRINT
PRINT "Theoretical probability of four-of-a-kind is "; USING "#.####"; 150 / 6 ^ 5
PRINT "Please wait patiently, until the first thousand rolls..."
PRINT
PRINT "Attempt #", "Probability (Simulation estimate)"
RANDOMIZE TIMER

'Repeat process ten times
FOR m = 1 TO 10
total = 0
'Simulate 1000 rolls of five dice
FOR n = 1 TO 1000
'Generate random integer xi, 0<xi<7, i=1,...,5
x1 = INT(1 + RND(1) * 6)
x2 = INT(1 + RND(1) * 6)
x3 = INT(1 + RND(1) * 6)
x4 = INT(1 + RND(1) * 6)
x5 = INT(1 + RND(1) * 6)
z = CountEq(x1, x2, x3, x4, x5)
IF z = 4 THEN total = total + 1
NEXT n
'Print result
PRINT m, total / 1000
NEXT m

FUNCTION CountEq (x1, x2, x3, x4, x5)
'Count how many of five numbers are the same
DIM q(5)
q(1) = x1
q(2) = x2
q(3) = x3
q(4) = x4
q(5) = x5
c = 0
FOR j = 1 TO 5
        ck = 0
        FOR k = 1 TO 5
            ck = ck + (q(j) = q(k))
        NEXT k
        ck = -ck
        IF ck > c THEN c = ck
NEXT j
'assign value to function
CountEq = c
END FUNCTION