From: Lukáš Jiřiště <lukas.jiriste@datapartner.cz>
Date: Fri, 5 Dec 2025 13:03:04 +0000 (+0100)
Subject: Save a nice simulation setting
X-Git-Url: https://git.ljiriste.work/?a=commitdiff_plain;h=9792c7526ba322bbb1b5e2163052aad7fa758993;p=nbody.git

Save a nice simulation setting
---

diff --git a/animate_result.py b/animate_result.py
index fc49565..79d411e 100755
--- a/animate_result.py
+++ b/animate_result.py
@@ -41,13 +41,13 @@ class MatAnimation:
     def show(self):
         self.animation = animation.FuncAnimation(
             self.fig, self.update, len(self.data["time"]),
-            interval = 100,
+            interval = 10,
             blit=False);
         plt.show();
 
     def update(self, frame):
         for line in self.lines:
-            from_frame = 0;
+            from_frame = frame - 100;
             to_frame = frame;
             line.line.set_data([
                 self.data[line.name_x][from_frame:to_frame],
@@ -77,7 +77,7 @@ def matAniFromFile(filename):
 
 
 ani = matAniFromFile(sys.argv[1]);
-for i in range(1,11):
+for i in range(1, 11):
     ani.plot(
         "p[" + str(i) + "].pos[1]",
         "p[" + str(i) + "].pos[2]",
diff --git a/nbody.mo b/nbody.mo
index d7f81b1..96d1700 100644
--- a/nbody.mo
+++ b/nbody.mo
@@ -32,7 +32,7 @@ package nbody
 	end unitvec;
 
 	model System
-		parameter Integer	N = 4;
+		parameter Integer	N = 10;
 		Planet				p[N];
 		RelPosition			pos[N, N];
 
@@ -41,9 +41,9 @@ package nbody
 			p[i].pos[1] = sin(i);
 			p[i].pos[2] = cos(i);
 			p[i].pos[3] = sin(3 * i);
-			p[i].speed[1] = cos(2 * i);
-			p[i].speed[2] = sin(3 * i);
-			p[i].speed[3] = cos(i);
+			p[i].speed[1] = cos(2 * i) / 10;
+			p[i].speed[2] = sin(3 * i) / 10;
+			p[i].speed[3] = cos(i) / 10;
 		end for;
 
 	equation