function y1=qh(a1,b1)

y1=a1+rand*(b1-a1);

end