% 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
 u=heaviside(1-sqrt(x^2+y^2))*besselj(0,5.5207*sqrt(x^2+y^2))
ezsurf(u,[-1.5,1.5],[-1.5,1.5])
axis([-1.5,1.5,-1.5,1.5,-1 ,1 ])
pause
u_n=u*cos(5.5207*t)
% en t=0
ezsurf(subs(u_n,t,0),[-1.5,1.5],[-1.5,1.5])
axis([-1.5,1.5,-1.5,1.5,-1 ,1 ])
pause
% en t=0.1
ezsurf(subs(u_n,t,0.1),[-1.5,1.5],[-1.5,1.5])
axis([-1.5,1.5,-1.5,1.5 ,-1 ,1 ])
pause
% en t=0.2
ezsurf(subs(u_n,t,0.2),[-1.5,1.5],[-1.5,1.5])
axis([-1.5,1.5,-1.5,1.5,-1 ,1 ])
pause
% en t=0.3
ezsurf(subs(u_n,t,0.3),[-1.5,1.5],[-1.5,1.5])
axis([-1.5,1.5,-1.5,1.5,-1 ,1 ])
pause
moviein(20);
for j=1:81
   dt=(j-1)*0.1
w=subs(u_n,t,dt);
ezsurf(w,[-1.5,1.5],[-1.5,1.5])
axis([-1.5,1.5,-1.5,1.5,-1 ,1 ])
    % pause
   M(:,j)=getframe;
end
 
