关于各向异性渗透的问题

[mikewg]

此问题我搜索了很多资料，通常都是渗透过程中，渗透率是各向同性的，这个问题通过pdetool是很容易求解的：
可以使用传热模型求解，这里热导率是各向同性的（如钢铁）
对于稳态方程：利用Elliptic函数，可以模拟整个区域的温度场。

但是现在我想求解各向异性渗透率：也就是说在传热模型中，设置c（热导率）为各向异性的，对于二维来说：
我可以假设在x方向cx＝1，在y方向cy＝2.那么我怎么模拟在整个区域的温度场呢？

在pdetool中也有关于diffusion的求解，但是D值不知道可不可以设置成各向异性的，
我尝试了下，对于二维可以设置成[1,0;0,2]但是报错

希望大家能够给出见解，十分感谢！