求助 求教用fem解旋转对称热传导方程的方法

paddycy

是作业啦，可是搞不定，大侠们帮帮忙吧！
要求是写一个matlab的程序，用有限元方法，解一个热传导问题，区域有旋转对成性。
我是用了decsg和initmesh来生成网格的，刚度阵打算用一维存储，
可是套上自己构造的算例，解出来的就是不对……
对于最后解线性方程的部分我还是有一点信心的，
问题应该是出在中间生成刚度阵和右端项的时候了。
哪位能指导一下这个过程么？谢谢！
如果能给在下一个程序参考，更加感激：）
在下附上自己写的那个，望有前辈指正！多谢！
这个当中的含有我自己用的例子，
答案应该是u=2*z+sin(z),
子午面的右边一半是[0,1]*[0,pi],横坐标r，纵坐标z，
对z＝0用第一类边界条件u＝0，
对z＝pi和r＝1的边界用第三（二？）类边界条件：
(u的外法向偏导)+3*u=g(z),
其中，g=6z+3sin(z) （在r＝1处）； g=6*pi+1（在z＝pi处）
再次感谢！！