function y = qd(t)
s = (1-exp(-t)).*sin(0.5*pi*t);
y = [s; s];
endfunction