clear;
f=@(x) sin(2*pi*x);
l=@(t) 0*t;
r=@(t) 0*t;
M = [10 50]; N = [10 50]; D = 0.05;

for i = 1:length(M)
    figure(i)
    w = burgers(0,1,0,1,M(i),N(i),f,l,r,D);
    fstr = sprintf('problem_8_4_1_%d',i);
    tstr = sprintf('Problem 8.4.1, h = 1/%d, k = 1/%d',M(i),N(i));
    title(tstr)
    print(fstr,'-dpng')
    print(fstr,'-dpdf')
end