% OCW Universidad de Cantabria
% Cálculo Simbólico y Numérico en Ecuaciones Diferenciales
% © Cursos 2001-2014,  M. Eugenia Pérez, UC
%
syms x y  t r theta
r=sqrt(x^2+y^2);
theta=atan(y/x);
formamembra_en_t=heaviside(1-r)*besselj(3,13.0152*r)*cos(3*theta)*cos(13.0152*t)
 % 
 u=subs(formamembra_en_t,t,0)% 
ezsurf(u,[-1.2,1.2,-1.2,1.2])
axis([-1.2,1.2,-1.2,1.2,-0.5,0.5])
pause
figure(2)
moviein(80);
for j=1:81
   dt=(j-1)*0.1
w=subs(formamembra_en_t,t,dt);
ezsurf(w,[-1.2,1.2],[-1.2,1.2])
axis([-1.2,1.2,-1.2,1.2,-0.5,0.5])
    % pause
   M(:,j)=getframe;
end
   
