%
%  wave plot script
%
%  hw09a.m -- wave equation with zero BCs
%

disp(' ')
disp('wave equation experiment')
disp('   work = 0/1/2 (poor/better/best)')

work = 1

Lx =  6;  Nx =  48*(2^work);  dx = Lx/Nx;
Ly =  8;  Ny =  64*(2^work);  dy = Ly/Ny;
Lt =  8;  Nt = 100*(2^work);  dt = Lt/Nt;

Ntsave = Nt/25;  lam2 = (dt/dx)^2
pause(2)

% coordinates

xc =  0:dx:Lx;  yc =  0:dy:Ly;  [x y] = meshgrid(xc,yc);

% sample colorplot

load u1 -ascii

clf;figure(1)
pcolor(xc,yc,u1);axis equal;axis([0 Lx 0 Ly])
shading flat;caxis([-1 1]);colorbar
hold on;
title('\bf waves');xlabel('\bf x-axis');ylabel('\bf y-axis')
pause(0.0)