% OCW Universidad de Cantabria
% Cálculo Simbólico y Numérico en Ecuaciones Diferenciales
% © Cursos 2001-2014,  M. Eugenia Pérez, UC
%
function y = viajaonda(c)
 %lee c "velocidad"
syms x t
%  
u=1/(((x+c*t)^2+1)) % 
% 
figure(1)
ezsurf(u,[0,2],[-8,8])
pause
%  
pause % 
figure(2)
%  
w1=subs(u,t,0)%  
ezplot(w1,[-8,8])
 AXIS([-8 8 0 2])
pause
figure(3)
w2=subs(u,t,0.1)%  
ezplot(w2,[-8,8])
 AXIS([-8 8 0 2])
pause
w3=subs(u,t,0.3)% 
ezplot(w3,[-8,8])
 AXIS([-8 8 0 2])
pause
w4=subs(u,t,1)% 
ezplot(w4,[-8,8])
 AXIS([-8 8 0 2])
pause
w8=subs(u,t,1.8)%  
ezplot(w8,[-8,8])
 AXIS([-8 8 0 2])
pause
 %
figure(4)
ezplot(w1,[-8,8])
 AXIS([-8 8 0 2])
hold on
ezplot(w2,[-8,8])
 AXIS([-8 8 0 2])
 ezplot(w3,[-8,8])
  AXIS([-8 8 0 2])
  ezplot(w4,[-8,8])
   AXIS([-8 8 0 2])
   ezplot(w8,[-8,8])
    AXIS([-8 8 0 2])
hold off

