|
 
![geneticdrift := proc (num) local a, b, i, numgen; a := [seq(op([A, B]), i = 1 .. num)]; print(a); numgen := 1; do b := [seq(op(a), i = 1 .. 2)]; a := combinat:-randcomb(b, (1/2)*nops(b)); print(a); if combinat:-numbperm(a) = 1 then break end if; numgen := numgen+1 end do; print("After ", numgen, "generations one genotype has become extinct") end proc:](/view.aspx?sf=143661/454231/c06a16f87ebd2e2dc59179cc13770265.gif)
[A, B, A, B, A, B, A, B] [B, B, A, B, A, A, A, B] [B, B, A, A, B, A, B, B] [B, A, A, B, B, B, B, A] [A, A, A, A, B, B, B, B] [A, A, A, A, A, B, B, B] [A, A, B, B, A, A, B, B] [A, B, A, B, A, B, B, B] [B, B, A, B, B, A, B, B] [A, B, B, A, B, B, B, B] [A, B, A, B, B, A, B, B] [A, B, B, A, B, B, A, B] [B, A, B, B, A, B, B, A] [A, A, B, A, B, A, B, A] [B, A, B, B, B, A, A, A] [A, A, A, A, B, B, A, A] [A, A, A, A, A, A, A, A]
|

|
(1) |


|