#Serie de V=C*phi en (-a, a) desarrollada en sen(m*pi/a*phi)

#Parametros
V0 = 1 ;
a = pi ;
C = V0/(2*a)

#Cantidad de puntos phi y numero de terminos de la sdf
Nsum = 150 ;
Npphi = 10000 ;

#Vectores de coordenadas
phivec = linspace(-3*a, 3*a, Npphi) ;
phivecREP = -a + mod((phivec+a), 2*a) ;

#Extension de la funcion V
V = C*phivecREP ;

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

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

#Factor
sdf = (V0/pi)*sdf ;

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