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