% Script to run burgers.m
t_0 = 0; 
t_f = 1; 
M = 1000; 
N = 256;
[u,x,t,kk,amp] = burgers(t_0,t_f,M,N);
for i=1:41; 
   plot(x,u(:,25*(i-1)+1)); 
   hold on; 
   plot(x+2*pi,u(:,25*(i-1)+1)); 
   title('cos initial data, shocks at t = 1'); 
   xlabel(t(25*(i-1)+1)); 
   hold off; 
   figure(1); 
   axis([0,4*pi,-1.5,1.5]); 
   pause(1); 
end
