% solucion de la ecuacion de ondas bi-dimensional
% cuando la condicion inicial es una autofuncion
% velocidad inicial cero

clear all

L = 1; H = 1;
n = 2;m = 2;

anm = 1; 

x = [0:0.05:L];
y = [0:0.05:H];

[xx,yy] = meshgrid(x,y);

cosx = cos(n*pi*xx);
cosy = cos(m*pi*yy);
coscos = cosx .* cosy;

for t=0:.02:4
  u = anm * cos(sqrt((n*pi)^2+(m*pi)^2)*t) * coscos;
 % contour(x,y,u, [-0.01 0.01]);
 % pause
  surfc(x,y,u);
  axis([0 L 0 H -1 1])
  if (t==0)
    pause
  else
    pause(0.1)
  end
end