% Domain 
N = 8;
n = (0:N-1)';

% Real signals
x1 = randn(N,1);
x2 = randn(N,1);

% DFT of signals
X1 = fft(x1);
X2 = fft(x2);

% Product of DFTs
X = X1.*X2;

% Circular convolution
x = ifft(X);

figure;
subplot(3,1,1);  stem(n,x1);
subplot(3,1,2);  stem(n,x2);
subplot(3,1,3);  stem(n,x);