%% Guia 3 - Ejercicio 1 - Fisica 2 (1er Cuatrimestre 2018)

L=12;%Largo de la cuerda
c=1; %Velocidad de fase
w1=c*pi/L; %frec primer modo
w2=2*w1; %frec segundo modo
A0=1;%Amplitud inicial de ambos modos

psi=@(x,t) A0*(sin(pi*x/L)*cos(w1*t)+sin(2*pi*x/L)*cos(w2*t)); %Solucion del problema

x=linspace(0,L,400); %Vector que recorre el largo de la cuerda
tau=2*pi/w1; %Periodo del movimiento
t=linspace(0,2*tau,300);%Vector de tiempos

H=zeros(length(t),length(x)); %Defino matrix auxiliar H
for i=1:length(t)
    for j=1:length(x)
        H(i,j)=psi(x(j),t(i));
    end
end
%H(i,:) Vector fila que en cada columna tiene psi(x) para el tiempo t(i)
y1=min(H(:)); %Extremos que uso para fijar los limites del grafico en y
y2=max(H(:));


%%
%Grafico

figure(2);
for i=1:length(t)
    figure(2)
    
   
   % plot(x,H(i,:),x,sin(pi*x/L)*cos(w1*t(i)),'--',x,sin(2*pi*x/L)*cos(w2*t(i)),'--')
   % plotea la solucion y los modos 1 y 2 por separado
   
   
   % plot(x,sin(pi*x/L)*cos(w1*t(i)),'--')
   % plot(x,sin(2*pi*x/L)*cos(w2*t(i)),'--')
    plot(x,H(i,:),'LineWidth',4)
    
    title('Ejercicio 1 - Guia 3','FontSize',22)
    xlim([0 L]);
    ylim([y1 y2]);
   
    xlabel('x','FontSize',18)
    ylabel('\psi (x)','FontSize',18)
    
   hold on
    plot(x,sin(pi*x/L)+sin(2*pi*x/L),'--') %Deja la condicion inicial
    plot(x,0*x,'--','LineWidth',1.2) %eje x
    hold off
    pause(0.04)
end