hx[x_,a_,b_]:=100 Sqrt[s[a,b][[1,1,2]][x]];

"The integrand 1./(Sqrt[1+650.099/(1+x)]<<1>>) has evaluated to
non-numerical values for all sampling points in the region with
boundaries {{15546.6,30000}}"