%% programa cuerda/presion 02

clear all;
close all;
clc;

%% datos inciales
L = 1;          %longitud de la cuerda
A = 1/4;        %amplitud
c = 1;          %velocidad de propagacion de ondas
x = 0:L/100:L;  %vector de distancia
Nx = length(x); %longitud de vector x

%% ciclo principal de animacion

N = 100;          %numero de terminos de la sumatoria

pxt(1:Nx) = 0;   %incializacion de presion

for t = 0:0.01:20
    
    pxt(1:Nx) = 0;   %incializacion de presion
    
    for n = 1:N
        kn(n) = ( (2*n -1)*pi )/(2*L);
        wn(n) = c*kn(n);
        d = 0.5;
        b = d/2;
        wnd(n) = sqrt( wn(n)^2 - b^2 );
        pxt(1:Nx) = pxt(1:Nx) - (4*A)/(wnd(n)*(2*n-1)*pi)*...
                                 exp(-b*t)*cos(kn(n)*x)*sin(wnd(n)*t);
                                 
    end;
    
    figure(1)
    plot(x,pxt)
    title('onda cuerda');
    xlabel('x (m)');
    ylabel('pxt (m)')
    axis([0 L -1.1*A 1.1*A])
    grid on;
    box on;
    
    pause(0.01);
    
end;