#Serie de Fourier V=pulso(x, a) (extension impar en (-a,a)) desarrollada en sen(m*pi/a*x)

#Parametros
V0 = 1 ;
a = 5 ;


#Cantidad de puntos x y numero de terminos de la sdf
Nsum = 150 ;
Npx = 10000 ;

#Vectores de coordenadas
xvec = linspace(-3*a, 3*a, Npx) ;
xvecREP = -a + mod((xvec+a), 2*a) ;

#Extension de la funcion V
V = V0*sign(xvecREP) ;

#Vector de coeficientes de la sdf
lind = 0:Nsum ;
kl = (2*lind+1)*pi/a ;
Cl =  1./(2*lind+1) ;

#Inicializacion de la sdf
sdf = zeros(1,length(xvec)) ;
#Evaluacion para cada x
for ix=1:length(xvec)
    sdf(ix) = sum( Cl.*sin(kl*xvec(ix)) ) ;
end

#Factor
sdf=4*V0/pi*sdf ;

#Grafica
figure
plot(xvec, sdf, "color", "red", "linewidth", 1.5)
hold
plot(xvec, V, ".", "color", "blue", "linewidth", 1.5)
grid
set(gca(), "fontsize", 18)