Qmd1asymp[x_,r_] := ((1-r)/(r E^z-1)*Exp[-z*x])/. FindRoot[r(E^z-1)-z,{z,(3-r)(1-r)Log[10,10/r]}, AccuracyGoal->15,MaxIterations->100]