Numerical Investigation on Dirac Delta Function