We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent a75aa2e commit ca6f92dCopy full SHA for ca6f92d
src/+otp/+nbody/+presets/Pleiades.m
@@ -0,0 +1,18 @@
1
+classdef Pleiades < otp.nbody.NBodyProblem
2
+ methods
3
+ function obj = Pleiades()
4
+ params.spacialDim = 2;
5
+ params.masses = 1:7;
6
+ params.gravitationalConstant = 1;
7
+ params.softeningLength = 0;
8
+
9
+ positions = [3; 3; 3; -3; -1; 2; -3; 0; 2; 0; -2; -4; 2; 4];
10
+ velocities = [0; 0; 0; 0; 0; 0; 0; -1.25; 0; 1; 1.75; 0; -1.5; 0];
11
+ y0 = [positions; velocities];
12
13
+ tspan = [0, 3];
14
15
+ obj = obj@otp.nbody.NBodyProblem(tspan, y0, params);
16
+ end
17
18
+end
0 commit comments