# onda quadra 'pari'
n = 11
A.oq.pari <- rep(0,n)
phi.oq.pari <- rep(0,n)    # sono nulli! (messa per generalità)
for (k in seq(1,n,by=2)) {
  A.oq.pari[k] <- 2/(pi*k)*sin(k*pi/2)
}

T.oq <- 1 # secondo (arbitrario...)

t <- seq(-2*T.oq, +2*T.oq, by=T.oq/1000)
nt <- length(t)

fk <- matrix(rep(0, n*nt), c(nt,n))
ft <- rep(0,nt) 

for(k in 1:n) {
  fk[,k] <- A.oq.pari[k]*cos(k*(2*pi/T.oq)*t)
  ft     <-  ft + fk[,k]
}

plot(t,ft, ty='l')