То ли MATLAB неправильно считает, то ли я чего-то не понимаю.  Попробуем решить в MATLAB вот такую простенькую систему одномерных дифференциальных уравнений:      с краевыми условиями 3-го рода    и однородными начальными условиями    Фактически начальное условие для функции v  определяется по начальному условию для функции u .  Первое уравнение параболического типа, второе — эллиптического.  Применим встроенный решатель MATLAB, а именно функцию pdepe. В документации сказано, что эта функция может решать системы уравнений параболического и эллиптического типа, в которых есть по крайней мере одно параболическое уравнение.  Но оказывается, что эта функция считает абсолютно неправильно.  
Блог Глеба Гренкина