Программа на с++ или Python для решения задачи о работе ударной трубы

В файлах прикреплено пособие
Страницы 66-69: постановка задачи, точное решение, задание расчетной области.
Страница 80: на ней указаны параметры тестовой задачи, по которым надо строить решения
Страницы 113-123: выбрать одну из неявных разностных схем для нахождения приближенного решения. 
Нужно выбрать одну из схем 18, 19, 21-24
Нужно написать программу на С++ или Python(лучше Python), которая будет находить точное решение задачи, приближенное решение по одной из схем(на ваше усмотрение); 
Построить графики решений как на странице 81(точное решение и приближенное на одном графике для каждого параметра)
Найти максимальное отклонение приближенного решения от точного по каждому параметру
Оставить комментарии в коде что за что отвечает и что означает
Сделать записи на листочке о том, как была преобразована схема для использования в программе
Вечер 1-го июня крайний срок выполнения.