Maybe some fun, maybe not

[romichess]

I was inspired today to write this silly simulation when I watched Fellippe's pandemic simulation video. It is not about the pandemic. It is about something different. It is not really a true simulation. It is a silly simulation. And hopefully it may also be fun. Some silly fun!

Code: QB64: [Select]

1. TYPE OBJECT
2.     sex AS INTEGER
3.     age AS INTEGER
4.     x AS SINGLE
5.     y AS SINGLE
6.     vx AS SINGLE
7.     vy AS SINGLE
8.     dt AS INTEGER
9. END TYPE
10.  
11. DIM SHARED object(1000) AS OBJECT
12.  
13. DIM SHARED sw AS INTEGER
14. DIM SHARED sh AS INTEGER
15.  
16. sw = _DESKTOPWIDTH
17. sh = _DESKTOPHEIGHT
18.  
19. SCREEN _NEWIMAGE(sw, sh, 256)
20.  
21. _FULLSCREEN
22.  
23. Initialize
24.  
25. repeat = 1
26.  
27. WHILE repeat
28.     CLS
29.  
30.     FOR i = 0 TO 999
31.         IF object(i).age > 6000 THEN
32.             object(i).sex = 0
33.             _CONTINUE
34.         END IF
35.         IF object(i).sex = 1 AND object(i).age > 960 AND object(i).age < 5400 THEN
36.             FOR j = 0 TO 999
37.                 IF object(i).age > 6000 THEN
38.                     object(i).sex = 0
39.                     _CONTINUE
40.                 END IF
41.                 IF object(j).sex = 2 AND object(j).age > 960 AND object(j).age < 2760 THEN
42.                     d = SQR((ABS(object(i).x - object(j).x) ^ 2) + (ABS(object(i).y - object(j).y) ^ 2))
43.                     IF d < 13 AND object(j).dt = 0 THEN
44.                         object(j).dt = 46
45.                         FOR k = 0 TO 999
46.                             IF object(k).sex = 0 THEN
47.                                 object(k).sex = RND * 2 + 1
48.                                 object(k).age = 0
49.                                 object(k).x = RND * (sw - 24) + 12
50.                                 object(k).y = RND * (sh - 24) + 12
51.                                 EXIT FOR
52.                             END IF
53.                         NEXT
54.                     END IF
55.                 END IF
56.             NEXT
57.         END IF
58.     NEXT
59.  
60.     FOR i = 0 TO 999
61.         IF object(i).sex = 0 THEN _CONTINUE
62.         object(i).age = object(i).age + 1
63.         IF object(i).dt > 0 THEN object(i).dt = object(i).dt - 1
64.         c = 11
65.         IF object(i).sex = 2 THEN c = 13
66.         object(i).vx = object(i).vx + RND * .1 - .05
67.         object(i).vy = object(i).vy + RND * .1 - .05
68.         object(i).x = object(i).x + object(i).vx
69.         object(i).y = object(i).y + object(i).vy
70.         IF object(i).x < 13 THEN object(i).vx = 1
71.         IF object(i).x > sw - 13 THEN object(i).vx = -1
72.         IF object(i).y < 13 THEN object(i).vy = 1
73.         IF object(i).y > sh - 13 THEN object(i).vy = -1
74.         CIRCLE (object(i).x, object(i).y), 12, c
75.         'PAINT (object(i).x, object(i).y), c
76.     NEXT
77.  
78.     IF INKEY$ = CHR$(27) THEN repeat = 0
79.  
80.     _LIMIT 60
81.  
82.     _DISPLAY
83. WEND
84.  
85. SYSTEM
86.  
87. SUB Initialize
88.     RANDOMIZE TIMER
89.     FOR i = 0 TO 999
90.         object(i).sex = 0
91.     NEXT
92.     FOR i = 0 TO 40
93.         object(i).sex = 1
94.         object(i).age = RND * 100
95.         object(i).x = RND * (sw - 24) + 12
96.         object(i).y = RND * (sh - 24) + 12
97.  
98.     NEXT
99.     FOR j = i + 1 TO i * 2
100.         object(j).sex = 2
101.         object(j).age = RND * 100
102.         object(j).x = RND * (sw - 24) + 12
103.         object(j).y = RND * (sh - 24) + 12
104.     NEXT
105. END SUB
106.  

[romichess]

Sixty frames = one simulated year.

Changed the mechanics a bit. Added more color information.
Green - too young
Grey - too old
Pink - female
Blue - male
Solid Pink - take a guess

Things that affect the outcome.
Population density - the larger the starting population the more likely the simulation will sustain
Size of the circles - the larger they are the more plentiful will be the collisions
CONST min - The lower the more fruitful the society will be
CONST maxm - The higher "
CONST maxf - The higher "

The goal is to achieve simulation sustainability without too explosive growth and with not too many failures.

DISCLAIMER ABOUT THE CURRENT SETTINGS
It is just a numbers game. The settings are restricted to what is commonly possible. They are not 100% realistic. It is not meant to be politically correct. So just consider these settings being for a primitive society without written laws or moral standards.

The morally correct may want to increase the circle size and increase the minimum age! Etc.

Code: QB64: [Select]

1. TYPE OBJECT
2.     sex AS INTEGER
3.     age AS INTEGER
4.     x AS SINGLE
5.     y AS SINGLE
6.     vx AS SINGLE
7.     vy AS SINGLE
8.     dt AS INTEGER
9. END TYPE
10.  
11. DIM SHARED object(1000) AS OBJECT
12.  
13. CONST min = 750
14. CONST maxm = 4800
15. CONST maxf = 2760
16.  
17. DIM SHARED sw AS INTEGER
18. DIM SHARED sh AS INTEGER
19.  
20. sw = _DESKTOPWIDTH
21. sh = _DESKTOPHEIGHT
22.  
23. SCREEN _NEWIMAGE(sw, sh, 256)
24.  
25. _FULLSCREEN
26.  
27. Initialize
28.  
29. repeat = 1
30.  
31. WHILE repeat
32.     CLS
33.  
34.     FOR i = 0 TO 999
35.         IF object(i).age > 6000 THEN
36.             object(i).sex = 0
37.             _CONTINUE
38.         END IF
39.         IF object(i).sex = 1 AND object(i).age > min AND object(i).age < maxm THEN
40.             FOR j = 0 TO 999
41.                 IF object(j).age > 6000 THEN
42.                     object(j).sex = 0
43.                     _CONTINUE
44.                 END IF
45.                 IF object(j).sex = 2 AND object(j).age > min AND object(j).age < maxf AND object(j).dt = 0 THEN
46.                     d = SQR((ABS(object(i).x - object(j).x) ^ 2) + (ABS(object(i).y - object(j).y) ^ 2))
47.                     IF d < 13 THEN object(j).dt = 46
48.                 END IF
49.                 IF object(j).dt = 1 THEN
50.                     object(j).dt = 0
51.                     twins = RND * 250
52.                     count = 0
53.                     IF twins = 1 THEN count = 1
54.                     FOR l = 0 TO count
55.                         FOR k = 0 TO 999
56.                             IF object(k).sex = 0 THEN
57.                                 object(k).sex = RND * 2 + 1
58.                                 object(k).age = 0
59.                                 object(k).x = RND * (sw - 24) + 12
60.                                 object(k).y = RND * (sh - 24) + 12
61.                                 object(k).dt = 0
62.                                 EXIT FOR
63.                             END IF
64.                         NEXT
65.                     NEXT
66.                 END IF
67.             NEXT
68.         END IF
69.     NEXT
70.  
71.     FOR i = 0 TO 999
72.         IF object(i).sex = 0 THEN _CONTINUE
73.         object(i).age = object(i).age + 1
74.         IF object(i).sex = 2 AND object(i).dt > 1 THEN object(i).dt = object(i).dt - 1
75.         c = 2
76.         IF object(i).sex = 1 THEN
77.             IF object(i).age > min THEN c = 11
78.             IF object(i).age > maxf THEN c = 7
79.         END IF
80.         IF object(i).sex = 2 THEN
81.             IF object(i).age > min THEN c = 13
82.             IF object(i).age > maxf THEN c = 7
83.         END IF
84.         object(i).vx = object(i).vx + RND * .1 - .05
85.         object(i).vy = object(i).vy + RND * .1 - .05
86.         object(i).x = object(i).x + object(i).vx
87.         object(i).y = object(i).y + object(i).vy
88.         IF object(i).x < 13 THEN object(i).vx = 1
89.         IF object(i).x > sw - 13 THEN object(i).vx = -1
90.         IF object(i).y < 13 THEN object(i).vy = 1
91.         IF object(i).y > sh - 13 THEN object(i).vy = -1
92.         CIRCLE (object(i).x, object(i).y), 12, c
93.         IF object(i).sex = 2 AND object(i).dt > 0 THEN PAINT (object(i).x, object(i).y), c
94.     NEXT
95.  
96.     IF INKEY$ = CHR$(27) THEN repeat = 0
97.  
98.     _LIMIT 60
99.  
100.     _DISPLAY
101. WEND
102.  
103. SYSTEM
104.  
105. SUB Initialize
106.     RANDOMIZE TIMER
107.     FOR i = 0 TO 999
108.         object(i).sex = 0
109.     NEXT
110.     FOR i = 0 TO 100
111.         object(i).sex = 1
112.         object(i).age = RND * 6000
113.         object(i).x = RND * (sw - 24) + 12
114.         object(i).y = RND * (sh - 24) + 12
115.  
116.     NEXT
117.     FOR j = i + 1 TO i * 2
118.         object(j).sex = 2
119.         object(j).age = RND * 6000
120.         object(j).x = RND * (sw - 24) + 12
121.         object(j).y = RND * (sh - 24) + 12
122.     NEXT
123. END SUB
124.  

[AndyA]

That's an interesting take on population growth and decline.

Maybe the solid pink circle should change to solid pink filled ellipses for a short time then be accompanied by smaller green circles. Then like all the others turn gray.