用迎风离散格式解对流方程（matlab）

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