%function LB = vlb(x,psil)
function LB = vlb(x,psil)
nu_mul = zeros(7,1);
ul = -10*ones(15,1);
xl = x;
LB = [nu_mul;ul;xl;psil];