Модель колебательного звена с применением символьного и численного решений дифференциального уравнения на SymPy и NumPy

Почему вы не оформите формулы в LaTeX? SymPy, кстати, тоже вроде умеет экспортировать для латеха.

Ради интереса не пробывали сравнить с результатми PyDy.

> Процесс решения дифференциальных уравнений становиться наглядным и хорошо контролируемым на каждом этапе вычислений

А вот статье и не скажешь, ничего не наглядно. Даже понять какую задачу человек решает невозможно.

Полученная графическая иллюстрация работы колебательного звена полностью соответствует теории.

Почему? И какой теории? А в каком смысле соответствует? Графики визуально похожи? Или есть какая-то метрика? А что мы вообще решаем?

Хотя бы общую постановку задачи надо было описать — а то статья превращается в «вот код, он работает и, поверьте мне, работает правильно».