delm=.1;
delmn=delm/4.34;
alfa=3;%valoarein dB;
alfa1=alfa/4.34;
n=100000;
j=.005;
sigma=10^(-9.5);

pas=.1:.1:1; 
m=length(pas);
r=ones(1,m);
for i=1:m,
   r(i)=(1/(2*alfa1))*log((n^(.5)*j*2*pas(i)*delmn)/(sigma*1.96*(1+exp(4*alfa1*pas(i)))^(.5)))+pas(i);
end
r1=alfa.*r;
plot(pas,r1)