I have a question to solve as shown below. I have tried to solve it by the codes written in Mathematica but something is wrong with it. I am very appreciated if get some help.
Thanks,
SNRdB := Table [i, {i, -5, 45, 1}]
M1 := 4
Var := 4
SNR := N[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))
f4 = NIntegrate[f1*f0, {x, 0, Infinity}]
f5 = NIntegrate[f2*f4, {y, -Infinity, Infinity}]
Pb = (1/Sqrt[2*Pi])*f5
