dftcs_p.m

clear;  %%% Plot data for dftcs.f %%%
load tplot.dat; load xplot.dat; load ttplot.dat;
subplot(121)
  mesh(ttplot,[-70 30]);
  xlabel('x');  ylabel('time');
  title('Diffusion of a delta spike');
subplot(122)
  cs = contour(ttplot,0:.5:5,tplot,xplot);
  xlabel('time');  ylabel('x'); clabel(cs,0:5);
  title('Contour plot');
subplot(111)