Wolfram Cloud Document

In[]:=

$Version

Out[]=

12.1.0 for Microsoft Windows (64-bit) (March 14, 2020)

In[]:=

SNRdB=Table[i,{i,-5,45,1}];M1=4;Var=4;SNR:=10^(SNRdB/10);f0=(1/(SNR*2*Var))*Exp[-x/(SNR*2*Var)];f1=Exp[-((y-Sqrt[2*x])^2/2)]*f0;f2=1-((1-Erfc[y/Sqrt[2]]/2)^(M1-1));

In[]:=

Pb=(1/Sqrt[2*Pi])*NIntegrate[f2*f1*f0,{y,-Infinity,Infinity},{x,0,Infinity},Method"LocalAdaptive"];

In[]:=

ListLogPlot[Transpose[{SNRdB,Pb}],JoinedTrue]

Out[]=