breed [bunnies]
bunnies-own [fur allele1 allele2]
globals [furallele nofurallele breed?]

to setup
  ca
  reset-ticks
  ask patches [; create background
    set pcolor 9.9 - (random 5)
    ]
  diffuse pcolor 1
  create-bunnies amount-fur-alleles [; creates alleles
    set size 1.5
    set color brown
    set shape "circle"
    setxy random-xcor random-ycor
    set fur true
    ]
  create-bunnies amount-no-fur-alleles [
    set size 1.5
    set color pink
    set shape "circle"
    setxy random-xcor random-ycor
    set fur false
    ]
   set furallele count bunnies with [fur = true]; counts alleles for first generation
   set nofurallele count bunnies with [fur = false]
end 

to generation
  ask bunnies [die]
  set breed? true; starts loop
  while [breed? = true] [
    create-bunnies 1 [
      set size 2
      set shape "rabbit"
      setxy random-xcor random-ycor
      if furallele > 0 and nofurallele > 0 [
        ifelse (random (furallele + nofurallele)) < furallele
          [set allele1 true]
          [set allele1 false]
        ]; sets allele1 random depending on possible alleles left
      if furallele > 0 and nofurallele = 0 [set allele1 true]
      if furallele = 0 and nofurallele > 0 [set allele1 false]
      if allele1 = true [set furallele (furallele - 1)]; removes allele1 from possible alleles left
      if allele1 = false [set nofurallele (nofurallele - 1)]
      if furallele > 0 and nofurallele > 0 [
        ifelse (random (furallele + nofurallele)) < furallele
          [set allele2 true]
          [set allele2 false]
        ]; sets allele2 random depending on possible alleles left
      if furallele > 0 and nofurallele = 0 [set allele2 true]
      if furallele = 0 and nofurallele > 0 [set allele2 false]
      if allele2 = true [set furallele (furallele - 1)]; removes allele2 from possible alleles left
      if allele2 = false [set nofurallele (nofurallele - 1)]
      ifelse allele1 = true or allele2 = true [set color brown][set color pink]; sets color depending on alleles
      ]
    if (furallele + nofurallele) = 0 [set breed? false]; if no more bunnies can be created, stops cycle
    ]
  set furallele (furallele + count bunnies with [allele1 = true and allele2 = true] * 2); sets number of possible alleles for next generation
  set furallele (furallele + count bunnies with [allele1 = true and allele2 = false] + count bunnies with [allele1 = false and allele2 = true])
  set nofurallele (nofurallele + count bunnies with [allele1 = true and allele2 = false] + count bunnies with [allele1 = false and allele2 = true])
  tick 
end