function y = delta(x,h)

if abs(x) >= 2*h
   y = 0;
else
   y = (1+cos((pi*x)/(2*h)) )/(4*h);
end