function y = xtt(tt,a,T,delay);
x = zeros(length(tt),1);
for (i = 1:2:length(a))
  x = x + a(i)*(u_step(tt,(i-1)*T+delay) - u_step(tt,(i-1)*T + 2+delay)); 
end;
y = x;
end;




