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