% 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 = 1; m = 1;

anm = 1; 

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

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

senox = sin(n*pi*xx);
senoy = sin(m*pi*yy);
sensen = senox .* senoy;

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