% OCW Universidad de Cantabria
% Cálculo Simbólico y Numérico en Ecuaciones Diferenciales
% © Cursos 2001-2014,  M. Eugenia Pérez, UC
%
function y= ondasestacionarias(n)
%
syms x t
u_n=sin(n*pi*x)*cos(3*n*pi*t)
ezsurf(u_n,[0,2],[0,1])
%
pause
ezplot(subs(u_n,t,0),[0,1])
axis([0,1,-2,2])
pause
% 
ezplot(subs(u_n,t,0.1),[0,1])
axis([0,1,-2,2])
pause
% 
ezplot(subs(u_n,t,0.2),[0,1])
axis([0,1,-2,2])
pause
%
ezplot(subs(u_n,t,0.3),[0,1])
axis([0,1,-2,2])
pause
% 
ezplot(subs(u_n,t,0.4),[0,1])
axis([0,1,-2,2])
%
pause
ezplot(subs(u_n,t,1/2),[0,1])
axis([0,1,-2,2])
pause
% 
ezplot(subs(u_n,t,1/2+0.1),[0,1])
axis([0,1,-2,2])
pause
% 
ezplot(subs(u_n,t,1/2+0.2),[0,1])
axis([0,1,-2,2])
pause
% 
ezplot(subs(u_n,t,1/2+0.3),[0,1])
axis([0,1,-2,2])
pause
% 
ezplot(subs(u_n,t,1/2+0.4),[0,1])
axis([0,1,-2,2])
% 
pause
ezplot(subs(u_n,t,1),[0,1])
axis([0,1,-2,2])
pause
% 
ezplot(subs(u_n,t,3/2),[0,1])
axis([0,1,-2,2])
pause
% 
ezplot(subs(u_n,t,2),[0,1])
axis([0,1,-2,2])
pause
% 
ezplot(subs(u_n,t,2+1/2),[0,1])
axis([0,1,-2,2])
pause
% 
ezplot(subs(u_n,t,3),[0,1])
axis([0,1,-2,2])
pause

