Требуется решить задачу численно, результатом будет две матрицы T и Phi, пока решаю в маткаде используя конечные разности и найдя все эти производные вручную, но может в питоне есть способ попроще? Каппа, альфа и сигма - тензоры соответствующие, зависящие от координаты и температуры, функции найдены. Задаем начальную температуру, которая будет в нулевом столбце одинаковая и потенциал, первую итерацию, а дальше расчет оставшихся.

http://python.su/forum/topic/20185/
а вообще смотри в сторону библиотеки math она довольно много что умеет

насколько я понял это задача в частных производных. Методы
попроще:
http://sfepy.org/doc-devel/index.html
http://www.ctcms.nist.gov/fipy/
http://fenicsproject.org/
Мы пользуемся в основном фениксом.
В wolfram mathematica тоже очень просто pde решать. начиная с 10 там тоже сделали фемовский решатель.