samer@4: function x=closed2hopen(x)
samer@4: % closed2hopen - convert closed interval to half-open interval
samer@4: %
samer@4: % closed2hopen :: [[N,2]] -> [[N,2]]
samer@4: %
samer@4: % The second column is assumed to be the upper limit of a closed
samer@4: % interval. These upper limits are incremented to the next floating
samer@4: % point number so that the result can be interpreted as a half-open
samer@4: % interval equivalent to the original closed interval.
samer@4: 
samer@4: mx=x(:,2);
samer@4: x(:,2)=mx+eps(mx);