'   Every generation of red amoeba has   a  5% chance of spawning 5 more red amoeba,
'                                       a 10% chance of spawning 3 blue amoeba,
'                                       a  1% chance of doing both, and
'                                       a 84% chance of not spawning.
'   Every generation of blue amoeba has  a  3% chance of spawning 5 more blue amoeba,
'                                       a  2% chance of spawning 10 red amoeba,
'                                       a  1% chance of doing both, and
'                                       a 94% chance of not spawning.

'   Your petri dish starts with one red amoeba.

'   Can you show what is the expected number of red and blue amoeba observed before all amoeba die off?
'       Can you do it graphically too?

'   NB I've had an off-scale run of over the 960 generations displayed here.....

    nomainwin

    WindowWidth  =1000
    WindowHeight = 550
    UpperLeftX   =int( ( DisplayWidth  -WindowWidth)  /2)
    UpperLeftY   =int( ( DisplayHeight -WindowHeight) /2)

    graphicbox #w.g1,                   10,   52, 980, 400

    textbox    #w.tb1,                  10,   10, 980,  30
    button     #w.exit, "Exit", [quit], LR,   30,  10

    maxgen  =0
    maxpop  =1
    maxpopt =1

    open "Red & Blue Amoeba pool" for window as #w

    #w,    "trapclose [quit]"
    #w.g1, "down; fill white; size 1; flush"

    #w,    "font arial 12"

for i =1 to 100000
    #w.g1, "cls ; color blue;  line 0 200 960 200"
    gen  =0
    red  =1
    blue =0
    while ((red >0) or (blue >0)) and (gen <960)
        r =int( 100 *rnd( 1)) +1
        if  (r <= 5)               then red  =red  + 5
        if  (r >  5) and (r <=15)  then blue =blue + 3
        if  (r > 15) and (r <=16)  then red  =red  + 5: blue =blue +3
        if  (r > 16) and red >0    then red  =red   -1

        r =int( 100 *rnd( 1)) +1
        if (r <= 3)                then blue =blue + 5
        if (r >3) and (r <= 5)     then red  =red  +10
        if (r >5) and (r <=6)      then blue =blue + 5: red  =red  +10
        if (r >6) and blue >0      then blue =blue  -1

        gen =gen +1

        #w.g1, "color  red;  line "; gen; " 199 "; gen; " "; 199 - red *2
        #w.g1, "color blue;  line "; gen; " 201 "; gen; " "; 201 +blue *2

        maxpop =max( maxpop, red +blue)
        maxgen =max( maxgen, gen)

        #w.tb1, "Test #"; right$( "      " +str$( i), 6); " reached generation =";_
         right$( "      " +str$( gen), 6); " before all were dead. Best so far was a run of ";_
         right$( "      " +str$( maxgen), 6); " & max pop'n was ";_
         right$( "      " +str$( maxpop), 6)

        scan
    wend

    if maxpop >maxpopt then
        maxpopt =maxpop
        #w.g1, "getbmp drawing 1 1 719 399"
        s$ ="AmoebaPic" +str$( maxpop) +".bmp"
        bmpsave "drawing", s$
        timer 5000, {MoveOn2]
        wait    '   allow time to save bmp
     [MoveOn2]
        timer 0
        unloadbmp "drawing"
    end if

next i
    wait

[quit]
    close #w
    end