function u=BeamWarming(a,dt,n,x1,x2,m)
a=1;
dt=0.005;
n=101;
x1=0;
x2=1;
m=100;
format long;
h=(x2-x1)/(n-1)
for j=1n+2*m)
    [b]Sample Text[/b]  [b]Sample Text[/b]
end
u1=u0;
for k=1:m
    for[b]Sample Text[/b]      [b]Sample Text[/b]
        u1(i)=u0(i)-(dt/h)*a*(u0(i)-u0(i-1))-a*(dt/h)*(1-a*(dt/h))*(u0(i)-2*u0(i-1)+u0(i-2))/2
    end
    u0=u1;
end
x=0:0.01n-1)/100;
u=u1((m+1)m+n));
plot(x,u,'r');
axis([0 1 -0.5 1.5]);
hold on
y=1.*(x<=0.5)+0.*(x>0.5)
plot(x,y,'c')
title('Beam-Warming')
format short;
function ux=U(x)
format long;
if x<=0
    ux=1
else
    ux=0
end