function [num] = randDouble(lo, hi)
% num = random real number in the interval (lo,hi), exclusive

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