Chain reaction with moderator and the two isotopes of Uran

globals
[
  decay-constant number-of-simulations count-simulation average
]

Breed [ Neutrons   Neutron ]
Breed [ Moderators   moderator ]
Breed [ URAN235s  URAN235 ]
Breed [ URAN236s  URAN236 ]
Breed [ URAN238s  URAN238 ]
Breed [ DaugtherXs  DaugtherX ]
Breed [ DaugtherYs  DaugtherY ]   ; Programmet anvender kun en datterkerne, men i virkeligheden, burde der være to

Neutrons-own [ speed ]

;; procedure for setup button

to setup
  clear-all
  set number-of-simulations 10            ;; Her ændres på antallet af simuleringer som man midler over.
  set average 0
  set count-simulation 1
  set-default-shape Moderators "molecule water"
  set-default-shape Neutrons "dot"
  set-default-shape URAN235s "orbit 6"
  set-default-shape URAN236s "orbit 6"
  set-default-shape URAN238s "orbit 6"
  set-default-shape DaugtherXs "dot"
  setup-particles                                        ;; restarts tick-counter
end 

to setup-particles
  create-URAN235s number-URAN235 [setup-URAN235s ]
  create-URAN238s number-URAN238 [setup-URAN238s ]
  create-Moderators number-moderator [setup-moderators ]
  fire-Neutron
  reset-ticks
end 

to setup-Neutrons
  set color red
  set size 1
  set speed 2
  set heading random 360
end 

to setup-moderators
   set color orange
   set size 2
   setxy random-xcor random-ycor
end 

to setup-URAN235s
   set color blue
   set size 1
   setxy random-xcor random-ycor
end 

to setup-URAN238s
   set color yellow
   set size 1
   setxy random-xcor random-ycor
end 

to setup-URAN236
  set size 1
  set color green
end 

;; procedure for go/stop buttom

to go
  move-Neutrons
  decay
  react
  if (not any? Neutrons) and (not any? URAN236s)
  [ if-else count-simulation < number-of-simulations
    [ print ( number-URAN235 - count URAN235s ) / number-URAN235
      set count-simulation (count-simulation + 1)
      set average (average + count URAN235s)
      ask turtles [die]
      setup-particles
      clear-all-plots
    ]
    [ print ( number-URAN235 - count URAN235s ) / number-URAN235
      type "Middelværdi: "
      print (number-of-simulations * number-URAN235 - (average + count URAN235s)) / (number-of-simulations * number-URAN235)
      print "nu stopper programmet"
      stop
    ]
  ]
  tick
end 

to fire-Neutron           ; Her dannes et antal neutroner - evt. undersg om antallet er stabilt.
  create-Neutrons 10
  [
    setup-Neutrons
    setxy 0 0
    set heading random 360
  ]
end 

to move-Neutrons
  ask Neutrons
  [
  if-else not can-move? ( speed * 1 )
    [ die ]
    [forward (speed * 1 )]
  ]
end 

;; procedure for fusion of particles to something kin of compund nucleus

to react
  ask Neutrons with [speed = 2]
  [
   if any? URAN235s-here
    [
     if (random-float 1 < 0.1 )
     [
      hatch-URAN236s 1 [setup-URAN236]
        ask one-of URAN235s-here [ die ]
        ask one-of Neutrons-here [ die ]
     ]
    ]
  if any? URAN238s-here
    [
      die
    ]
  if any? Moderators-here
    [
      set speed  1
    ]
  ]
  ask Neutrons with [speed = 1]
  [
   if any? URAN235s-here
   [
     if (random-float 1 < 0.9 )
     [
       hatch-URAN236s 1 [setup-URAN236]
       ask one-of URAN235s-here [ die ]
       ask one-of Neutrons-here [ die ]
     ]
    ]
  ]
end 

to decay
  ask URAN236s[
    if (random-float 1 < 0.1 )[
      hatch-neutrons 2 [setup-Neutrons]
      set color white
      hatch-DaugtherXs 1
      die
    ]
  ]
end