clc;

%% Propagacion de ondas en una cuerda .....................................

h=0.005;        %delta x 
k=0.01;         %delta t 
c=1/3;          %velocidad en metros por segundo
lambda=c*k/h;   %lambda<1 para que sea convergente
L=1;            %largo de la cuerda en metros

x=0:h:L;

% Condiciones de contorno: cuerda fija en ambos extremos ..................
u_0=zeros(1,length(x)); 
u_1=zeros(1,length(x)); 
u_2=zeros(1,length(x)); 

% Simulacion basica con reflexion a los 2/3 de la cuerda ..................
figure(1);

% Parametros de perturbacion en el extremo izquierdo 
beta=0.002; 
i0 = 40;         %centro del pulso en tiempo 

% Reflexion 
R = 2/5;          %probar 1.5 
R = 1;
% R = 1.5;
for i=1:10000
    plot (x(1:2*length(u_2)/3),u_2(1:2*length(u_2)/3),'b');   %dibujo
    hold on
    plot (x(2*length(u_2)/3:length(u_2)),u_2(2*length(u_2)/3:length(u_2)),'r');   %dibujo 
    hold off
    axis([0 1 -2 2]);
    u_0=u_1;            %actualizo u_0
    u_1=u_2;            %actualizo u_1
    
    %Introduzco una perturbacion
    u_2(1)=exp(-beta.*((i-i0).^2)); 

    %Calculo el nuevo instante de tiempo u_2 
    for j=2:length(x)-1 
        if(j<length(x)*2/3)
          lambda=c*k/h;     %Velocidad c1
        else
          lambda=c*k/h*R;   %Velocidad c2=c1*R
        end
        u_2(j)=(lambda^2)*u_1(j+1)+2*(1-lambda^2)*u_1(j)+(lambda^2)*u_1(j-1)-u_0(j);
    end
    pause(0.001);
end