#*****************************************************
# programme de synthèse de Fourier
# Dominique Lefebvre  mars 2013
# www.tangenteX.com
#*****************************************************

# chargement des packages
with(plots):

# définition de la fonction creneau
creneau := t -> (-1)^floor(t/Pi):

# plot de la fonction creneau
p1 := plot(creneau,-4*Pi..4*Pi, color=black, numpoints=1000):

# plot de différentes sommes de sinusoïdes
#s0 := plot((4/Pi)*sin(t),t = -4*Pi..4*Pi, color = blue, numpoints=1000):
#s1 := plot((4/Pi)*(sin(t) + sin(3*t)/3),t = -4*Pi..4*Pi, color = green, numpoints=1000):
s2 := plot((4/Pi)*(sin(t) + sin(3*t)/3 + sin(5*t)/5),t = -4*Pi..4*Pi, color = red, numpoints=1000):
s10 := plot((4/Pi)*(sin(t) + sin(3*t)/3 + sin(5*t)/5 + 
       sin(7*t)/7 + sin(9*t)/9 + sin(11*t)/11 + sin(13*t)/13 +
       sin(15*t)/15 + sin(17*t)/17 + sin(19*t)/19),
       t = -4*Pi..4*Pi, color = gray, numpoints=1000):

display([p1,s10]);