function X = uniformPointsWithinBounds(N, LB, UB)
Unifs = rand(N, size(LB,2));
X = Unifs .* repmat(UB - LB, N, 1);
X = X + repmat(LB, N, 1);
end
