import numpy as np
import matplotlib.pyplot as plt
from scipy.fft import fft, ifft
from scipy import signal
t=np.arange(0,100,0.1)
x= np.cos(2*np.pi*t)
y= fft(x)
p= np.linspace(0,1000,1,endpoint=True)
u = signal.square( np.pi*5*p)
e = np.fft.fft(u)
plt.plot(p,e)