怎样编程求解二维偏微分方程

blueberry

遇到一个二维偏微分方程：
dy/dt = D(d^2y/dr^2+1/r*dy/dr+d^2y/dh^2)

本打算用matlab的pdetool解决，但是因为有一个边界条件是这样的：-D(dy/dr(r=R)+dy/dh(r=R))=Kdy/dt (k=const.)

因为这个边界条件在pdetool中无法输入，所以用准备用有限元差分方法编程，但是我实在缺乏编程技巧，没有能解决这个问题。

希望这里的大虾能不吝赐教一下，能有相似问题的程序供我参考一下更好，感激不尽。