b = 1;
mm = 10;
lambdaB = 1/(π*500^2);
lambdaD = 10/(π*500^2);
alpha = 4;
Noise = 3.9811*10^(-15);
T = 0.001;
Pc2 = 1180591620717411303424/4722366482869645;
Pd2 = 60095708769426201171875/
4611686018427387904 - (252401976831590044921875*Exp[(-16/5)/4611686018427387904]);
f = (1 + 1/((u^(-4)*mm*(1 + b) + 1 I)*(u^(-4)*mm*(1 - b) + 1 I))^(1/2))*u;
Pcn1 = (4722366482869645^(1/2)*π^(1/2))/34359738368;
R = (17179869184*4722366482869645^(1/2)*π^(1/2))/4722366482869645;
ICC[t_] := Exp (-2*π*lambdaB*(mm*t)^(2/alpha)*Pc2*
NIntegrate[f, {u, Pcn1*R/((t*mm)^(1/alpha)), ∞}])
IDC[t_] := Exp (-2*π*lambdaD*(mm*t)^(2/alpha)*Pd2*NIntegrate[f, {u, 0, ∞}])
TT[t_] := Im[Exp (1 I*t*Noise)/((-1 I*t/T*mm*(1 + b) + 1)*(-1 I*t/T*mm*(1 - b) + 1))^(1/2)*ICC[t]*IDC[t]]/t;
coverageCEUE = 1/2 - 1/π*NIntegrate[TT[t], {t, 0, ∞}] // N