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