[求助]二维椭圆型偏微分方程求解的问题

zerocpu

-▽·(c▽u)+au=f 形式的方程，用函数assempde(b,p,e,t,c,a,f)求解
其中c a f是偏微分方程的参数
b是边界条件
p e t是网格划分
具体的不是很清楚，看了一些实例，还是不太明白

我想问的是
1.看到的几个实例中c a f 都是常数，
不过我需要求解的方程是-▽·(c▽u)+au=f(x,y)的形式，并且f是一个离散函数，这个要怎么办？
(假设f(x,y)和已知矩阵M的元素一一对应)

2.边界条件的设置，实例中看到有这些
g='circleg';     %  定义单位圆区域
b=’circleb2’;  % x^2 on the boundary
另一种
g='squareg';     %  定义单位方形区域
b='squareb3';    % 左右零边界条件，顶底零导数边界条件
我希望是矩形区域，边界四周函数值都是0，应该怎么定义b？
还有g定义区域，单位方形的时候是不是x，y都是-1~1的范围？

3.网格划分
[p,e,t] = initmesh(g);
这个完全不明白，p,e,t分别代表什么参数？
如果我要按照上面的M矩阵的规模（比如M是256*256的矩阵）来划分网格，该怎么定义？

4.最后我怎么调出u(x,y)的值啊？
貌似实例中都将pedsurf（p，t，u）画个图，不过我要具体数值

初次接触PDE问题，同学介绍来这个论坛，还望各位高手能帮小弟解答一下。万分感谢！

zerocpu

额。。。不知道我有没有把问题讲清楚。。。

ChaChing

个人水平专业有限, PDE Toolbox没使用过!
最近工作很忙, 实在没空看help文件学习! sorry
耐心待高人路过

cjcfrank

可以查一下武汉大学陆君安老师的书。