function y = dmpcos(t, a, b) if nargin == 1, a = 1; b = 1; end y = ustep(t) .* exp(-a*t) .* cos(2*pi*b*t);