dirac delta function matlab code