function y = circconv(x1,x2)
% circconv(x1,x2) obtains the circular convolution of vectors
% x1 and x2.  x1 and x2 must be the same length.

X1 = fft(x1);
X2 = fft(x2);
Y = X1.*X2;
y = real(ifft(Y));
	
