function ret=hh(r)
if r<7
   ret=0.05*exp(r.^0.554);
else
   ret=exp(-9*(log10(r)-log10(9)).^(2.3));
end