一维热传导问题求解函数matlab

%%传热模型
clear ;
g = ‘squareg‘;                       %定义以为正方形
b = ‘squareb1‘;                       %定义边界条件为零
c = 1; a = 0; f = 0; d = 1;
%产生初始三角网格
[pet]=initmesh（g）;
%定义初始条件
u0 = zeros（size（p2）1）;
ix = find（sqrt（p（1:）.^2+