function c = fast_multiply(a,b)
	n = length(a);
	pad_a = [a, zeros(1,n)];
	pad_b = [b, zeros(1,n)];
	c = ifft( fft(pad_a) .* fft(pad_b) );
	c = c(1:2*n-1);
end
