>> [u1,u1_true,x1,t] = beam_warming(0,1,20,20); >> [u2,u2_true,x2,t] = beam_warming(0,1,40,40); >> [u3,u3_true,x3,t] = beam_warming(0,1,80,80); >> figure(1) >> plot(x1,u1_true(:,21)); >> hold on >> plot(x1,u1(:,21)); >> plot(x2,u2(:,41),'--'); >> plot(x3,u3(:,81),':'); >> axis([0,2*pi,-.2,1.2]); >> title('Beam Warming method','FontSize',16) >> xlabel('space','FontSize',16) >> xlabel('x','FontSize',16) >> gtext('solid: true solution and dt=.05 approximation'); >> gtext('dashed: dt=.025 approximation'); >> gtext('dotted: dt=.0125'); >> print -dps beam_warming.ps >> figure(2) >> [u3,u3_true,x3,t] = lax_wendroff(0,1,80,80); >> [u2,u2_true,x2,t] = lax_wendroff(0,1,40,40); >> [u1,u1_true,x1,t] = lax_wendroff(0,1,20,20); >> figure(2) >> plot(x3,u3(:,81),':'); >> hold on >> plot(x2,u2(:,41),'--'); >> plot(x1,u1(:,21)); >> plot(x1,u1_true(:,21)); >> axis([0,2*pi,-.2,1.4]); >> title('Lax Wendroff method','FontSize',16) >> xlabel('x','FontSize',16) >> gtext('dotted: dt=.0125'); >> gtext('dashed: dt=.025'); >> gtext('solid: true and dt=.05'); >> print -dps lax_wendroff.ps >> figure(3) >> [u1,u1_true,x1,t] = lax_friedrichs(0,1,20,20); >> [u2,u2_true,x2,t] = lax_friedrichs(0,1,40,40); >> [u3,u3_true,x3,t] = lax_friedrichs(0,1,80,80); >> plot(x1,u1_true(:,21)); >> hold on >> plot(x1,u1(:,21)); >> plot(x2,u2(:,41),'--'); >> plot(x3,u3(:,81),':'); >> figure(3) >> axis([0,2*pi,-.1,1.1]); >> title('Lax Friedrichs method','FontSize',16) >> xlabel('x','FontSize',16) >> gtext('solid: true and dt=.05'); >> gtext('dashed: dt=.025'); >> gtext('dotted: dt=.0125'); >> print -dps lax_friedrichs.ps >> diary off