Код на python для нелинейного уравнения(численное моделирование)

Распределение температуры в стержне с несколькими слоями, методом конечных разностей численно решить уравнение(есть файл в котором более подробно все описано: постановка задачи, даны прогоночные коэффициенты и тд). Далее надо будет сделать точное решение(тоже на python).

В целом сделать все то, что указано в файле используя python