%% programa cuerda/presion 01

clear all;
close all;
clc;

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

%% ciclo principal de animacion

N = 6;          %numero de terminos de la sumatoria

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

for t = 0:0.01:5
    
    pxt(1:Nx) = 0;   %incializacion de presion
    
    for n = 1:N
        pxt(1:Nx) = pxt(1:Nx) + (8*A)/(n^2*pi^2)*sin(n*pi/2)*sin(n*pi*x/L)*cos(n*pi*c*t/L);
    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;