import numpy as np
import matplotlib.pylab as plt
from scipy import signal

t=np.linspace(0,4,1000, endpoint=False)
tp=signal.square(np.pi*t)

ts=4/3*np.sin(np.pi*t)+(1/3)*np.sin(np.pi*3*t)+(1/5)*np.sin(np.pi*5*t)+(1/7)*np.sin(np.pi*7*t)+(1/9)*np.sin(np.pi*9*t)+(1/11)*np.sin(np.pi*11*t)+(1/13)*np.sin(np.pi*13*t)+(1/15)*np.sin(np.pi*15*t)+(1/17)*np.sin(np.pi*17*t)+(1/19)*np.sin(np.pi*19*t)+(1/21)*np.sin(np.pi*21*t)+(1/23)*np.sin(np.pi*23*t)+(1/25)*np.sin(np.pi*25*t)+(1/27)*np.sin(np.pi*27*t)+(1/29)*np.sin(np.pi*29*t)+(1/29)*np.sin(np.pi*29*t)+(1/31)*np.sin(np.pi*31*t)+(1/33)*np.sin(np.pi*33*t)

plt.plot(t,tp)
plt.plot(t,ts)
plt.show()