twpopuweb.nlogo

NETLOGO

globals [b b0 b1 b2 b3 b4   

;; System dynamics model globals

  ;; stock values
  Total-popu
  birth-rate
  year
  death-rate
  ;; size of each step, see SYSTEM-DYNAMICS-GO
  dt
  ;; total steps
  system-dynamics-t
]

;; Initializes the system dynamics model.
;; Call this in your model's SETUP procedure.
to system-dynamics-setup 
  set dt 1.0
  set system-dynamics-t 0
  ;; initialize stock values
  set Total-popu 2055
  set birth-rate 0.0157
  set year 80
  set death-rate 0.0052
end

;; Step through the system dynamics model by performing next iteration of Euler's method.
;; Call this in your model's GO procedure.
to system-dynamics-go
  set system-dynamics-t system-dynamics-t + dt

  ;; compute variable and flow values once per step
  let local-popu-birth popu-birth
  let local-birth-rate-change birth-rate-change
  let local-popu-death popu-death
  let local-year-ahead year-ahead
  let local-death-rate-change death-rate-change

  ;; update stock values
  ;; use temporary variables so order of computation doesn't affect result.
  let new-Total-popu max( list 0 ( Total-popu + local-popu-birth - local-popu-death ) )
  let new-birth-rate max( list 0 ( birth-rate - local-birth-rate-change ) )
  let new-year max( list 0 ( year + local-year-ahead ) )
  let new-death-rate max( list 0 ( death-rate + local-death-rate-change ) )
  set Total-popu new-Total-popu
  set birth-rate new-birth-rate
  set year new-year
  set death-rate new-death-rate
end

;; Report value of flow
to-report popu-birth
  report Total-popu * birth-rate  * dt
end

;; Report value of flow
to-report birth-rate-change
  report yearly-birth-rate
 * dt
end

;; Report value of flow
to-report popu-death
  report Total-popu * death-rate * dt
end

;; Report value of flow
to-report year-ahead
  report 1 * dt
end

;; Report value of flow
to-report death-rate-change
  report yearly-death-rate * dt
end

;; Plot the current state of the system dynamics model's stocks
;; Call this procedure in your model's GO procedure.
to system-dynamics-do-plot
  set-current-plot "populations"
  set-current-plot-pen "Total-popu"
  plotxy system-dynamics-t Total-popu
 ; set-current-plot-pen "birth-rate"
  ;plotxy system-dynamics-t birth-rate
  ;set-current-plot-pen "year"
  ;plotxy system-dynamics-t year
  ;set-current-plot-pen "death-rate"
  ;plotxy system-dynamics-t death-rate
end






;;------------------
to setup 
ca
system-dynamics-setup
end

to go
repeat 5 [
            system-dynamics-go
            system-dynamics-do-plot
          ]
end

to-report yearly-birth-rate 
ifelse year >= 80 and year <= 84 [report 0.0001]
[ifelse year >= 85 and year <= 89 [report 0.0007]
  [ifelse year >= 90 and year <= 91 [report 0.00085]
  
      [ifelse year >= 92 and year <= 92 [report 0.0]                                ;[report 0.0]
     [ifelse year >= 93 and year <= 99 [report 0.000257]                              ;[report 0.000257]
      [ifelse year >= 100 and year <= 109 [report 0.00023]                           ;[report 0.00023]
        [ifelse year >= 110 and year <= 119 [report 0.00012]                       ;[report 0.00012]
         [ifelse year >= 120 and year <= 139 [report 0.000075] [report (-1 * yearly-increased-birth-rate-after-140)]         ;[report 0.000075] [report 0.0]
  
;    [ifelse year >= 92 and year <= 92 [report b0]                                ;[report 0.0]
 ;    [ifelse year >= 93 and year <= 99 [report b1]                              ;[report 0.000257]
  ;    [ifelse year >= 100 and year <= 109 [report b2]                           ;[report 0.00023]
   ;     [ifelse year >= 110 and year <= 119 [report b3]                       ;[report 0.00012]
    ;     [ifelse year >= 120 and year <= 139 [report b4] [report 0.0]         ;[report 0.000075] [report 0.0]
             ]
           ]
        ]
      ]
    ]
  ]
]
end


to-report yearly-death-rate 
ifelse year >= 80 and year <= 84 [report 0.0001]
[ifelse year >= 85 and year <= 89 [report 0.0]
  [ifelse year >= 90 and year <= 91 [report 0.00005]
    [ifelse year >= 92 and year <= 92 [report 0.0]
      [ifelse year >= 93 and year <= 99 [report 0.0001143]
        [ifelse year >= 100 and year <= 109 [report 0.00008]
           [ifelse year >= 110 and year <= 119 [report 0.00026]
             [ifelse year >= 120 and year <= 139 [report 0.000405] [report yearly-increased-death-rate-after-140]
             ]
           ]
        ]
      ]
    ]
  ]
]
end
@#$#@#$#@
GRAPHICS-WINDOW
900
10
1002
133
1
1
30.7
1
1
1
1
1
0
1
1
1

CC-WINDOW
5
356
1011
451
Command Center
0

BUTTON
12
20
75
53
NIL
setup
NIL
1
T
OBSERVER
T
NIL

BUTTON
80
20
143
53
NIL
go
NIL
1
T
OBSERVER
T
NIL

PLOT
212
19
709
267
populations
NIL
NIL
0.0
10.0
0.0
10.0
true
true
PENS
"Total-popu" 1.0 0 -2674135 true

MONITOR
13
113
143
162
NIL
Total-popu
0
1

MONITOR
13
59
143
108
NIL
year
0
1

SLIDER
17
309
377
342
yearly-increased-birth-rate-after-140
yearly-increased-birth-rate-after-140
-2.0E-4
2.0E-4
7.0E-5
1.0E-5
1
NIL

SLIDER
375
308
713
341
yearly-increased-death-rate-after-140
yearly-increased-death-rate-after-140
-6.0E-4
6.0E-4
3.0E-4
1.0E-6
1
NIL

MONITOR
14
164
143
213
NIL
birth-rate
10
1

MONITOR
15
215
143
264
NIL
death-rate
10
1

TEXTBOX
91
284
309
310
チ瓣140