%function UB = vub(x,psiu)
function UB = vub(x,psiu)
nu_muu = ones(7,1);
uu = 10*ones(15,1);
xu = x;
UB = [nu_muu;uu;xu;psiu];