def delta_eps(t,eps):
d = np.zeros(len(t))<EOL>for k,tt in enumerate(t):<EOL><INDENT>if abs(tt) <= eps/<NUM_LIT>:<EOL><INDENT>d[k] = <NUM_LIT:1>/float(eps)<EOL><DEDENT><DEDENT>return d<EOL>
Rectangular pulse approximation to impulse function.

Parameters
----------
t : ndarray of time axis
eps : pulse width

Returns
-------
d : ndarray containing the impulse approximation

Examples
--------
>>> import matplotlib.pyplot as plt
>>> from numpy import arange
>>> from sk_dsp_comm.sigsys import delta_eps
>>> t = np.arange(-2,2,.001)
>>> d = delta_eps(t,.1)
>>> plt.plot(t,d)
>>> plt.show()