function r = randReal(lo,hi)
% r is a random real number in the interval (lo,hi)

r= rand(1)*(hi-lo) + lo;