Raket opsendelse

globals [mass acceleration velocity height dm/dt N]

to setup
  clear-all
  reset-ticks
  create-turtles 1 [
    set shape "rocket"
    set size 3
    set heading 0
    setxy 0 1
    set color red
  ]
  set velocity 0
  set acceleration 0
  set height 0
  set mass raket_masse + brændstofmasse      ; startsmasse
  set dm/dt masseændring_per_tid   ; masseændring  per tid
  set N 10
  setup-world
end 

to setup-world
ask patches [
    set pcolor 99.8
    if pycor < 1 [ set pcolor green]
  ]
end 

to go
  let deltat 1 / N
  repeat N [
  let thrust ( dm/dt * udstødningshastighed )
  ifelse mass >  raket_masse [
  set mass (mass - dm/dt * deltat)
  ]
  [ set dm/dt 0]
  set acceleration (thrust / mass -  9.82)
  set velocity (velocity + acceleration * deltat)
  set height (height + velocity * deltat)
  ]
  if height < 0 [setup-world stop]
  ask turtles [
    setxy xcor height / 1000
  ]
  paint-sky
  tick
end 

to paint-sky
    ask patches [ set pcolor 100 - ( min list  (height / 10000)  10 )]
end