matlab中圆柱体的三维扩散方程

Question

我想知道是否有办法解决PDE在0.3米半径和1.6米高的圆柱体中的扩散问题。我正在尝试使用此代码创建一个圆柱坐标。

%cylinder
% n = linspace(-pi,pi,30);
% m = linspace(0,1.6,20);
% [theta,z] = meshgrid(n,m);
% r = 1 .* cos(theta);
% [X,Y,Z] = pol2cart(theta,r,z); 
% surf(X,Y,Z)
% axis equal

但是这个我无法定义半径。

这是我必须使用的等式（不稳定状态-3d）。 我得到的代码只有2d平原，但我不能使用3D圆柱坐标。 我想与偏微分方程的离散化一起工作。 我如何用Matlab解决这个问题？ 谢谢