% initU.m
% initiates the 20 initial state vectors (rows)
%
% i = index of initial state to use

function u = initU (i)
if ~(i <= 20 && i >= 1)
    error('the specified index must be between 1 and 20 inclusive')
else
    init = [
        1, 1, 1, 0, 0, 1, 1, 0, 1;
        1, 0, 0, 1, 0, 1, 0, 0, 1;
        1, 1, 0, 0, 1, 0, 0, 0, 1;
        0, 0, 0, 1, 0, 1, 0, 1, 1;
        1, 0, 0, 1, 0, 0, 0, 0, 0;
        1, 1, 0, 1, 1, 0, 1, 1, 1;
        0, 0, 1, 0, 0, 0, 1, 1, 1;
        1, 1, 0, 1, 0, 0, 0, 1, 0;
        0, 1, 0, 1, 0, 1, 0, 0, 0;
        0, 1, 0, 1, 1, 0, 1, 0, 0;
        0, 0, 1, 0, 1, 0, 0, 0, 1;
        0, 1, 0, 1, 1, 1, 0, 1, 1;
        1, 0, 1, 1, 1, 1, 1, 1, 0;
        1, 1, 0, 0, 1, 0, 1, 0, 0;
        1, 0, 0, 1, 0, 1, 0, 0, 0;
        1, 0, 0, 1, 1, 0, 0, 0, 1;
        0, 0, 1, 0, 0, 0, 1, 0, 1;
        1, 0, 0, 1, 1, 0, 1, 0, 1;
        0, 0, 0, 0, 1, 0, 1, 1, 1;
        1, 0, 0, 1, 1, 0, 0, 1, 0];
    u = init(i, 1:9);
end