Всё это и многое другое хорошо реализовано (в смысле достаточно доступно, просто и ясно) в MathCAD.
Всё в смысле- решение методом Рунге-Кутта и Эйлера.
Для подробного изучения, а также реализации можно почитать литературу (например С.В.Поршнев Компьютерное моделирование физических процессов с использованием пакета MathCAD,
MathCad12 для студентов).
А лучше посетить
Математический образовательный сайт и на нем зайти на Форум по MathCAD.