RANDOMIZE
DEF FN_RND_INDEX = INT(2*RND+1) ! either returns 1 or 2,
! which is equivalent to heads and tails
DIM animal1$(2)
DIM animal2$(2)
DIM animal3$(2)
DIM animal4$(2)
DIM animal5$(2)
DIM animal6$(2)
DIM animal7$(2)
DIM animal8$(2)
DIM animal9$(2)
DIM animal10$(2)

DIM results(4)

LET animal1$(1) = "a" 
LET animal1$(2) = "b" 
LET animal2$(1) = "c" 
LET animal2$(2) = "d"

FOR i = 1 TO 100

LET animal3$(1) = animal2$(FN_RND_INDEX)
LET animal3$(2) = animal1$(FN_RND_INDEX)
   
LET animal4$(1) = animal2$(FN_RND_INDEX)
LET animal4$(2) = animal1$(FN_RND_INDEX)
   
LET animal5$(1) = animal4$(FN_RND_INDEX)
LET animal5$(2) = animal3$(FN_RND_INDEX)
   
LET animal6$(1) = animal4$(FN_RND_INDEX)
LET animal6$(2) = animal3$(FN_RND_INDEX)
   
LET animal7$(1) = animal6$(FN_RND_INDEX)
LET animal7$(2) = animal5$(FN_RND_INDEX)
   
LET animal8$(1) = animal6$(FN_RND_INDEX)
LET animal8$(2) = animal5$(FN_RND_INDEX)
   
LET animal9$(1) = animal8$(FN_RND_INDEX)
LET animal9$(2) = animal7$(FN_RND_INDEX)
   
LET animal10$(1) = animal8$(FN_RND_INDEX)
LET animal10$(2) = animal7$(FN_RND_INDEX)
   
LET a_count = 0
LET b_count = 0
LET c_count = 0
LET d_count = 0
IF animal9$(1) = "a" THEN LET a_count = a_count+1
IF animal9$(1) = "b" THEN LET b_count = b_count+1
IF animal9$(1) = "c" THEN LET c_count = c_count+1 
IF animal9$(1) = "d" THEN LET d_count = d_count+1
IF animal9$(2) = "a" THEN LET a_count = a_count+1
IF animal9$(2) = "b" THEN LET b_count = b_count+1
IF animal9$(2) = "c" THEN LET c_count = c_count+1 
IF animal9$(2) = "d" THEN LET d_count = d_count+1
IF animal10$(1) = "a" THEN LET a_count = a_count+1
IF animal10$(1) = "b" THEN LET b_count = b_count+1
IF animal10$(1) = "c" THEN LET c_count = c_count+1 
IF animal10$(1) = "d" THEN LET d_count = d_count+1
IF animal10$(2) = "a" THEN LET a_count = a_count+1
IF animal10$(2) = "b" THEN LET b_count = b_count+1
IF animal10$(2) = "c" THEN LET c_count = c_count+1 
IF animal10$(2) = "d" THEN LET d_count = d_count+1
   
LET diversity_count = 0
IF a_count > 0 THEN LET diversity_count = diversity_count + 1
IF b_count > 0 THEN LET diversity_count = diversity_count + 1
IF c_count > 0 THEN LET diversity_count = diversity_count + 1
IF d_count > 0 THEN LET diversity_count = diversity_count + 1
   
LET results(diversity_count) = results(diversity_count) + 1

NEXT i

FOR i = 1 TO 4

PRINT "Number of"; i;" alleles: "; results(i)

NEXT i

END