Dirac function (Dt) calculation

Readded the limit function to the code