Springmotions

globals [
  tyngde-acceleration
  y         ; lodret position
  vy        ; lodret hastighed
  ay        ; lodret acceleration
  dt        ; tidsskridt
  k         ; fjederkonstant
]

to nulstil
  clear-all

  ask patches [ set pcolor cyan ] ;
  create-turtles 1 ; opret en kugle og tildel startværdier
  [
    setxy 0 30 ; start position
    set shape "square"
    set color red
    set size 8
    set heading 180 ; nedadrettet bevægelse
  ]
  set dt 0.1
  set y ([ycor] of turtle 0)
  set vy 0
  set k 2
  draw-spring
  reset-ticks
end 

; beregn ny position og fart

to go
  set vy vy - y * 2 / masse-af-lod * dt    ;   Hookes lov accelerationen er a=F/m=k*x/m
  set y y + vy * dt
  ask turtles
  [
    setxy xcor y
  ]
  draw-spring
  tick
end 

; beregn kinetisk energi: Ekin = (1/2)*m*v^2

to-report Ekin
  report 1 / 2 * masse-af-lod * ( vy ^ 2 )
end 

to draw-spring
  clear-drawing
  create-turtles 1 [
    setxy 0 50
    pendown
    set pen-size 3
    set color blue
    setxy xcor y
    die
  ]
end