function [output] = boundary_condition_BSO(X,Lb,Ub,dim)
for i = 1:dim
    if X(i) > Ub(i)% lon hon can tren
        X(i) = Lb(i) + rand.*(Ub(i)-Lb(i));
    end
    if X(i) < Lb(i)% nho hon can duoi
        X(i) = Lb(i) + rand.*(Ub(i)-Lb(i));
    end
end
% X = min(X, Ub);
% X = max(X,Lb);
output = X;
end


