Numerical Integration at different point

syms g;

I0_dBm=-12:22;

for jj=1:length(I0_dBm)

    var_ni(jj)=(6.09*10^-10*(10^(I0_dBm(jj)/10))*g)+(8.28*10^-14);
    fun(jj)=(1/4)*(2.02*10^7)*1.3844*(g^(2.1417))*exp(-(10* (10^(I0_dBm(jj)/10))*(10^-3)*g/(2*sqrt(var_ni(jj))))^2)*erfc((log(g/(0.0068))+0.1821)/sqrt(0.1));
    ber(jj)=sym(int(fun(jj),g, 0, Inf));

end

figure

semilogy(I0_dBm,ber,'g-d','LineWidth',2);

I tried this integration but the error is-

• Error using semilogyConversion to double from sym is not possible.*

Answer 1

[Torsten]

I0_dBm=-12:22;

for jj=1:length(I0_dBm)
    var_ni=@(g)(6.09*10^-10*(10^(I0_dBm(jj)/10))*g)+(8.28*10^-14);
    fun=@(g)(1/4)*(2.02*10^7)*1.3844*(g.^(2.1417)).*exp(-(10* (10^(I0_dBm(jj)/10))*(10^-3)*g./(2*sqrt(var_ni(g)))).^2).*erfc((log(g/(0.0068))+0.1821)/sqrt(0.1));
    ber(jj)=integral(fun, 0, Inf);
end

figure

semilogy(I0_dBm,ber,'g-d','LineWidth',2);

[Pritam Sahoo]

Thank you Sir for your valuable time. I have tried your answer, but the present error is-

_ Nonscalar arrays of function handles are not allowed; use cell arrays instead._ Error in (line 146) var_ni(jj)=@(g)(6.09.*10.^-10.*(10.^(I0_dBm(jj)./10)).*g)+(8.28.*10.^-14);

[Torsten]

Where did you see var_ni(jj) in my code ?

[Pritam Sahoo]

Thank you very much Sir. I am very sorry. Can I have your name, affiliation and email ID. I think, if you have not any prolem i'll add your name as an author in my paper.

[Torsten]

That's very good of you, but not necessary.

Best wishes

Torsten.

[Pritam Sahoo]

Thank You Sir.