

t_0 = 0;
t_f = 1;
M = 100;
N = 100;


[u_eu,u_true,x,t] = explicit_upwind_periodic(@g,t_0,t_f,M,N);
[u_lf,u_true,x,t] = lax_friedrichs_periodic(@g,t_0,t_f,M,N);
[u_lw,u_true,x,t] = lax_wendroff_periodic(@g,t_0,t_f,M,N);
[u_bw,u_true,x,t] = beam_warming_periodic(@g,t_0,t_f,M,N);
for j=0:20
	subplot(2,1,1)
    plot(x,u_true(:,1+5*j),'r','LineWidth',2);
    hold on
    plot(x,u_eu(:,1+5*j),'LineWidth',2);
    axis([0,2*pi,-1.3,1.3]);
    title('Explicit Upwind','FontSize',16)
    xlabel('true solution is in red','FontSize',16)
    hold off
	subplot(2,1,2)
    plot(x,u_true(:,1+5*j),'r','LineWidth',2);
    hold on
    plot(x,u_lf(:,1+5*j),'LineWidth',2);
    axis([0,2*pi,-1.3,1.3]);
    title('Lax Friedrichs','FontSize',16)
    xlabel('true solution is in red','FontSize',16)
    hold off

    figure(1)
    pause(1)
end

pause

for j=0:20
	subplot(2,1,1)
    plot(x,u_true(:,1+5*j),'r','LineWidth',2);
    hold on
    plot(x,u_lw(:,1+5*j),'LineWidth',2);
    axis([0,2*pi,-1.3,1.3]);
    title('Lax Wendroff','FontSize',16)
    xlabel('true solution is in red','FontSize',16)
    hold off
	subplot(2,1,2)
    plot(x,u_true(:,1+5*j),'r','LineWidth',2);
    hold on
    plot(x,u_bw(:,1+5*j),'LineWidth',2);
    axis([0,2*pi,-1.3,1.3]);
    title('Beam Warming','FontSize',16)
    xlabel('true solution is in red','FontSize',16)
    hold off

    figure(1)
    pause(1)
end

