#
# function that output a random real number in [a,b]
#
real_ran:=proc(  a :: numeric,                 
                 b :: numeric
              )
   local x;
   x:=rand();
   return evalf( a+ (b-a)*x/999999999999 );
end;