function[out]=ut(x1,y1,d)
out=conv(urect(2*x1/d)/d*ucomb(x1/d,1),urect(2*y1/d)/d*ucomb(y1/d,1))
end