function Pb=myfun(BPH,n,rho_in_db)

%r=[5:0.5:9];
%p=[0.2148,0.1297,0.0723,0.0370,0.0168,0.0069,0.0024,0.0007,0.00018];
%semilogy(r,p);
%hold on
%r1=[5:1:8,8.5,9];
%p1=[0.0597,0.0268,0.0095,0.0017,0.0006,0.00015];
%semilogy(r1,p1,'^');
fanout=2^BPH;
u=10^(rho_in_db/10);
u1=BPH*u;
n1=2*n;
%g=inline('(1-exp(-x).*(1+x+(1/2).*x.^2)).*(1/2).*((x./(n.*u)).^((n-2)./4)).*exp((-1/2).*(x+n.*u)).*besseli(n./2-1,sqrt(n.*u.*x))','x','n','u');
%Q=quad(g,0,100,1e-3,[],n1,u1);
%E=1-Q;
%Pb=(fanout^n/2)*E;