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

     end;