Краткое описание

Mathematica 11.0 code

np = 200;
m = 1;
pos = Table[RandomReal[{-1, 1}, 2], {np}];
acc = Table[
   Sum[If[j == k, 
     0, -(1/m) (
      pos[[j]] - 
       pos[[k]])/((pos[[j, 1]] - pos[[k, 1]])^2 + (pos[[j, 2]] - 
          pos[[k, 2]])^2)^(3/2)], {k, 1, np}], {j, 1, np}];
dt = 0.01;
posold = pos;
pos = posold + acc/2 dt^2;
clip = 10^0;
range = 20;
plots = {
   Graphics[Point[posold] , PlotRange -> {{-10, 10}, {-10, 10}}],
   Graphics[Point[pos] , PlotRange -> {{-10, 10}, {-10, 10}}]
   };
Do[
  acc = Table[
    Sum[If[j == k, 
      0, -(1/m) (
       pos[[j]] - 
        pos[[k]])/((pos[[j, 1]] - pos[[k, 1]])^2 + (pos[[j, 2]] - 
           pos[[k, 2]])^2)^(3/2)], {k, 1, np}], {j, 1, np}];
  acc = Clip[acc, {-clip, clip}];
  posoldold = posold;
  posold = pos;
  pos = 2 posold - posoldold + acc dt^2;
  AppendTo[plots, 
   Graphics[Point[pos] , 
    PlotRange -> {{-range, range}, {-range, range}}]];
  , 250];

Лицензирование

Я, владелец авторских прав на это произведение, добровольно публикую его на условиях следующей лицензии:

Этот файл доступен на условиях Creative Commons CC0 1.0 Универсальной передачи в общественное достояние (Universal Public Domain Dedication).

Лица, связанные с работой над этим произведением, решили передать данное произведение в общественное достояние, отказавшись от всех прав на произведение по всему миру в рамках закона об авторских правах (а также связанных и смежных прав), в той степени, которую допускает закон. Вы можете копировать, изменять, распространять, исполнять данное произведение в любых целях, в том числе в коммерческих, без получения на это разрешения автора. http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse

Этот файл, изначально опубликованный на внешнем сайте https://twitter.com/j_bertolotti/status/1037655339968147456, был проверен 19 октября 2018 проверяющим Ronhjones, подтвердившим, что файл был там доступен в ту дату на условиях указанной лицензии.