%%傳熱模型
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+