function [A, C] = createConstraints(sigma0, epsilon_threshold, sigma_threshold)
A = [];
C = [];


for i = 1:14
    line = zeros(1, 14);
    if i ~= 7 && i ~= 14
        line(i) = 1.00001;
        line(i+1) = -1.0;
    end
    if i == 7 || i == 14
        line(i) = 1; 
    end
    
   A =  [A;line];
end

% sigma1 > sigma0
line = zeros(1,14);
line(1) = -0.9999;
A = [A; line];



%% Do Constraints vector now
for i = 1:6
    C(i) = 0;
end
C(7) = sigma_threshold; %max value allowed

for i = 8:13
    C(i) = 0;
end
C(14) = epsilon_threshold;

C(15) = -sigma0;
