Dirac function (Di) calculation

Removed the limit function due to matlab compatibility and will look into another method to add it or the effect this will cause