#Serie de V=V0*cos(2pi*n*x/a) (extension impar en (-a,a)) desarrollada en sen(m*pi/a*x)

#Parametros
V0 = 1 ;
a = 5 ;
n = 3 ;

#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 = sign(xvecREP).*(V0*cos(2*n*pi*xvecREP/a)) ;

#Vector de coeficientes de la sdf
lind = 0:Nsum ;
kl = (2*lind+1)*pi/a ;
Cl =  (2*lind+1)./( (2*(lind+n)+1).*(2*(lind-n)+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)