%% programa solución de ecuacion de onda 01

clearvars;
close all;
clc;

%este programa grafica y anima la solucion de la ecuacion de onda
%mediante su solucion usando serie de fourier

%% datos inciales
L = 1.00;       %longitud de la cuerda
A = 1/4;        %amplitud
c = 344;        %veocidad de propagacion de ondas
x0 = L/2;       %punto fijo del espacio
fs = 44100;
t = 0:1/fs:1;
Nt = length(t);

%% ciclo principal de calculo

N = 100;  %numero de terminos de la sumatoria

p = 0;   %incializacion de presion

for nt = 1:Nt
    
    p(nt) = 0;   %incializacion de presion
    
    for n = 1:N
        
        
        p(nt) = p(nt) + (8*A)/(n^2*pi^2)*sin(n*pi/2)*sin(n*pi*x0/L)*cos(n*pi*c*t(nt)/L);
    end
    
end
%% resultados

figure(1)
plot(t,p)
title('onda cuerda');
xlabel('t (s)');
ylabel('p(x_0,t) (Pa)')
legend(['p(x_0,t)']);
axis([min(t), max(t), -1.1*A, 1.1*A])
grid on;
box on;


sound(p,fs);