function x=TwoLoop(a,t);
%TwoLoop.m--A double loop
%

n=length(a);
N=2*n-2;
p=length(t);dt=t(2)-t(1);
x=zeros(p,1);
f=2*pi/(N*dt);
for j=1:p
	for k=1:n
		x(j)=x(j)+a(k)*f*(k-1)*t(j);
	end
end