%
%   Delta Function
%

function Cd = cdelta(t)
    Cd = t*0;
    i = find(t>=0);
    Cd(i(1))=1/(t(2)-t(1));
end