Написать программу для решения системы ДУ на Python

Необходимо написать программу на Python, которая решает систему диф.ур. методом Рунге-Кутты 4-ого порядка. Имеется лит-ра (15 страниц), которая описывает данную систему. И пример программы на Pascal ,на котором можно базироваться при выполнении работы.  Необходимо реализовать ее с вводом параметров (теплоемкости, массы воды, коэф. теплопередачи и тд) на входе и получением результатов и графиков на выходе. Система представлена на картинке.