clear; clc; close all;

% dimension of the constellation
dimension = 2;

% set S of possible noise values
noise_levels = [-2; 0; 2];
S = combn(noise_levels, dimension)';

% set P of possible constellation points
constellation_points = [-3; -1; 1; 3];
P = combn(constellation_points,dimension)';

for j = 1 : size(P,2) % over all symbols
    
    % initialize the empty sets
    C  = []; % set C_j
    SC = []; % set S \ C_j
    
    
    %%% here: fill the sets C and SC %%%
    
    
    C % print the set C_j
    SC % print the set SC_j
    
%     run s4p1fg % uncomment this for parts f) and g)
    
end
