function?[?ux?]?=?burgers（?dt?N?T）
%求解無(wú)粘burgers方程的一階差分法
%???使用迎風(fēng)差分格式?將C轉(zhuǎn)到matlab里面
pi=acos（-1）;
dx=2*pi/N;
eps=1e-16;
for?i=1:1:N+1
????x（i）=（i-1）*dx;
end
for?i=1:1:N+1
????u（1i）=cos（x（i））;?%初始函數(shù)
end
for?j=2:1:T
????for?i=1:1:N+1
????????if?（i==1）?pre=N;?else?pre=i-1;end
?????