1 ' 3D gravitation model
\r
2 ' made by Svjatoslav Agejenko
\r
4 ' H-Page: svjatoslav.eu
\r
5 ' E-Mail: svjatoslav@svjatoslav.eu
\r
8 DECLARE SUB adda (x!, y!, z!, s!)
\r
9 DECLARE SUB display ()
\r
10 DIM SHARED atx(1 TO 100)
\r
11 DIM SHARED aty(1 TO 100)
\r
12 DIM SHARED atz(1 TO 100)
\r
13 DIM SHARED atxs(1 TO 100)
\r
14 DIM SHARED atys(1 TO 100)
\r
15 DIM SHARED atzs(1 TO 100)
\r
16 DIM SHARED ats(1 TO 100)
\r
18 DIM SHARED myx, myy, myz
\r
19 DIM SHARED oax(1 TO 100)
\r
20 DIM SHARED oay(1 TO 100)
\r
21 DIM SHARED oas(1 TO 100)
\r
32 adda RND * 6 - 3, RND * 6 - 3, RND * 4, 50
\r
38 IF INKEY$ <> "" THEN SYSTEM
\r
41 SUB adda (x, y, z, s)
\r
62 x1 = x / z * 100 + 160
\r
63 y1 = y / z * 100 + 100
\r
65 CIRCLE (oax(a), oay(a)), oas(a), 0
\r
66 CIRCLE (x1, y1), ats(a) / z, 15
\r
87 IF b = a THEN GOTO 2
\r
88 v = SQR((atx(b) - x) ^ 2 + (aty(b) - y) ^ 2 + (atz(b) - z) ^ 2)
\r
91 pxs = pxs + (atx(b) - x) / v2 / 10000
\r
92 pys = pys + (aty(b) - y) / v2 / 10000
\r
93 pzs = pzs + (atz(b) - z) / v2 / 10000
\r
98 atxs(a) = atxs(a) / 1.01 + pxs
\r
99 atys(a) = atys(a) / 1.01 + pys
\r
100 atzs(a) = atzs(a) / 1.01 + pzs
\r
105 atx(a) = atx(a) + atxs(a)
\r
106 aty(a) = aty(a) + atys(a)
\r
107 atz(a) = atz(a) + atzs(a)
\r