#Dada la solución analítica vamos a graficar para distintos tiempos. 
#La solición será una combinación de senos finita
L=1; #largo de la soga
K=pi/L; 
#Defino el vector x en donde calculo la solución, tengo que elegir un número discreto de valores, pueden ver cómo cambia al cambiar el último valor
X=linspace(0,L,1000); #es un vector desde 0 a L con 1000 valores equiespaciados

#Defino una función: dado un vector x, un dado k y un wt calcula la función de onda
function [Y]=onda(x,k,wt) #devuelve los valores de la función, la llamé Y
  Y=sin(k*x)*cos(wt)+(1/3)*sin(3*k*x)*cos(3*wt)+ (1/5)*sin(5*k*x)*cos(5*wt);
endfunction

#Vamos a correr la función para distintos tiempos-veo "fotos" de la cuerda.
#La forma es la siguiente: creo un vector de wt (V,B y C) y realizo un barrido en ese vector:
#calculo la función onda en el vector X, con k=K definido al ppcio y para cada wt de V.

V=[0,1/5,1/3,1/2,1]; #Valores de w1t -a un factor pi- para los que me pide graficar 
for w=V
    w1t=pi*w; #defino wt
    [y]=onda(X,K,w1t); #; es para que no lo muestre en la terminal
    hold on #Le pido que espere así grafica todas las curvas juntas
    plot(X,y) #grafico y en función de X
    xlabel('X*L')
    ylabel("Cuerda")
    legend("0","1/5","1/3","1/2","1")
endfor
hold off

#Simetría -impar- alrededor de ?/2
B=[0,1/3,1/2,2/3,1]
for w=B
    w1t=pi*w;
    [y]=onda(X,K,w1t);
    hold on
    plot(X,y)
    xlabel('X*L')
    ylabel("Cuerda")
    legend("0","1/3","1/2","2/3","1")
endfor
hold off
#Por ejemplo, w1t=?/3 y w1t= 2?/3 equidistan de ?/2 en ?/6.

#Simetría -par- alrededor de w1t=?. 
C=V=[0,1/3,1,5/3,2]
for w=C
    w1t=pi*w;
    [y]=onda(X,K,w1t);
    hold on
    plot(X,y)
    xlabel('X*L')
    ylabel("Cuerda")
    legend("0","1/3","1","5/3","2")
endfor
hold off
#w1t=0 y w1t=2? equidistan en ? a w 1 t=?. w1t=?/3 y w1t=5?/3 equidistan de w1t=? en 2?/3.
  
  