w = .0002; k = Pi/w;
                [Eta] = 1*10^-3; bk = Cosh[k*h0]^2/Sinh[2*k*h0]; tg = (
                2*Pi*[Eta]*Bk)/([Gamma]*k^2);

                 cdg[t_?NumericQ] :=
                20*Quiet[NIntegrate[
              Exp[-(x^2)/8]*(1/(1 + B/x^2))*(1 - Exp[-t/tg*x^4*(1 +B/x^2)])/
              x, {x, 0, [Infinity]}, AccuracyGoal -> 5]]