План работы по моделированию движения в поле тяжести с помощью Python :

1. Выбор пакета для моделирования. На GitFlic есть пакет pySimSheduler , который позволяет рассчитывать и моделировать взаимодействие нескольких тел под влиянием силы тяжести. Он написан на библиотеках Pymunk, Scipy.odeint, matplotlib и numpy.
2. Определение входных данных задачи. Нужно задать начальные условия и другие параметры движения, например, начальную скорость и угол наклона тела.
3. Определение числа интервалов сетки, на которой будет искаться приближённое решение.
4. Выделение памяти под массив сеточных значений решения системы ОДУ. В строке с номером m этого массива хранятся сеточные значения решения, соответствующие моменту времени t_m.
5. Задание начальных условий. Записываются значения в строку с номером 0 массива u.
6. Реализация схемы Эйлера. Для этого нужно выполнить цикл, в котором u[m + 1] = u[m] + tau*f(u[m], g, k, mass).
7. Отрисовка решения. Нужно построить график, который отображает траекторию движения тела.