%Practica 3-Ecuaciones Diferenciales y Metodos Numericos. Curso 1011
%%Sobre ondas y vibraciones 
%funcion dalembert1011H% a partir de la version 7 de MATLAB
%Dibuja la forma de la cuerda de longitud infinita en el instante t%%para t a determinar en clase de practicas
%con condiciones inciales f y g=0,
%f(x) dada en el ejercicio clase %%%%ecuacion u_{tt}-25 u_{xx}=0
syms  x;
f=-sym('Heaviside(x+pi)'*(x^2-pi^2)-'Heaviside(x-pi)'*(x^2-pi^2))
%
%instante incial t=0
ezplot(f,[-18,18])
AXIS([-18 18 0 15])%ajustando ejes
pause
%t=0.25 
%w es la solucion en t=0.25
w2=subs(f,x,x-5*0.25);
w1=subs(f,x,x+5*0.25);
w=(w1+w2)/2;
ezplot(w,[-18,18])
AXIS([-18 18 0 15])
pause

% superficie solucion
syms x t
solucion=(subs(f,x,x+5*t)+subs(f,x,x-5*t))/2;
ezsurf(solucion,[0,3],[-18,18])

 

