% rows 1 to 7 are sigma1-7, rows 8 to 14 are epsilon 1-7
function [x0] = createInitialValues(sigma0, epsilon_threshold, sigma_threshold, length)
x0 = [];

length = length +  mod(length,2);
for i=1:length/2
    x0(i) = sigma0 + (sigma_threshold - sigma0) * (2*i/length);
end

for i=length/2 + 1:length
    x0(i) = epsilon_threshold * (2*i/length - 1);
end