function?u=peYF（adtnminxmaxxM）?
%%方程中的常數；時間步長；空間節點數；求解區間的左端；求解區間的右端；迭代次數
format?long;
h=（maxx-minx）/（n-1）;
if?a>0
???for?j=1:（n+M）
??????u0（j）=IniU（minx+（j-M-1）*h）;%向左延拓M個節點的函數值
???end
else
???for?j=1:（n+M）
??????u0（j）?=?IniU（minx+（j-1）*h）;%向左延拓M個節點的函數值
???end
end
u