def METHOD_NAME():
	N = 2048
	t = np.linspace(0,1,N)
	Tj = [0.1,0.13,0.15,0.23,0.25,0.40,0.44,0.65,0.76,0.78,0.81]
	Hj = [4,-5,3,-4,5,-4.2,2.1,4.3,-3.1,2.1,-4.2]
	def K(t):
		return (1 + np.sign(t))/2.0
	y = np.zeros(t.shape)
	for tj,hj in zip(Tj,Hj):
		y = y + hj*K(t - tj)
	return t,4*y