无粘burgers方程差分求解matlab代码

function [ ux ] = burgers（ dt N T）
%求解无粘burgers方程的一阶差分法
%   使用迎风差分格式 将C转到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））; %初始函数
end
for j=2:1:T
    for i=1:1:N+1
        if （i==1） pre=N; else pre=i-1;end