Dirac function (Dt) calculation

I found that matlab has a built in dirac function, to avoid a name collision I renamed this function in case possibly a want occurs to use this alternative one.  I also updated the code to work with the limit function.