Метод конечных элементов

[_sliver_]

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

Я пишу свою программу для решения двумерного нестационарного уравнения теплопроводности. Со стационарным решением я разобрался, а для нестационарного вычисляю матрицу демпфирования, и далее действую как написано у Сегерлинда.
Но в результате итераций получается совсем не то решение какое должно быть.

Заранее спасибо за помощь!