Симуляция физических явлений с VPython

[N-Cube]

Полученная визуализация совсем не впечатляет. А вот, скажем, ноутбук Jupyter с примером визуализации с помощью библиотеки VTK: 3D Density Inversion

[netricks]

А систему тел со связями можно промоделировать?

[mphys]

Статью не читал ( :-) ) но мне кажется пакет просто решает систему дифф. уравнений. Сможете написать систему — будет моделирование.

[netricks]

Систему я и без библиотеки написать могу...

[omxela]

Ну, предположим, можете. Хотя тут без оговорок не обойтись, имхо. Типа, а в какой системе координат мы эту систему запишем? К Солнцу привяжем? Оно, конечно, массивное относительно. Но тем самым мы выкинем некоторое дрожание центра масс системы связанных тел. Ладно, записали систему. Расставили планеты, задали начальные скорости. Предположим, все это у нас есть. Но тут начнется хорошо известный цирк, связанный с тем, что в такой системе некое малое рассогласование начальных условий может вызвать качественное различие итоговых стационарных решений. В таком случае просто решать задачу Коши не получится. Это не значит, что все пропало, шеф, но нужны спецметоды численного интегрирования плюс меры по демпфированию зависимости от ошибок начальных условий. Это я к тому, что вряд ли это есть в предлагаемом ресурсе. Хотя, кто его знает, может, и есть. Тогда я буду сильно потрясен.

[lz961]

В том то дело, что не решает. Автор метод Эйлера ручками пишет. Т.е. берем интерпретатор, который сам по себе для вычислений не пригоден в виду низкой производительности, цепляем к ниму средство визуализации, и используем самую неэфективную численную схему для решения ОДУ.

[kovserg]

Засимулируте нам шаровую молнию. Чем не физическое явление?

[ku4in]

У Вас опечатка вначале в коде, пропущен перенос строки:

import vpython as vpvp.cylinder() 

И еще непонятно, зачем два раза силу считать, если можно просто знак изменить:

star.force = gravitationalForce(star,planet)
planet.force = gravitationalForce(planet,star)

Хотя, этот вопрос больше к оригиналу.