function?[pnextrcount]=koskie_gajic1
sir_tar=9;%目標信干比
k=5000;%bi/2ci
w=5e+6;%擴頻帶寬
eps=1e-15;%迭代誤差
v=1e+4;%傳輸速率
q=5e-15;%背景噪聲
n=1:1:50;
d=301:14:1000;%隨機50個用戶離基站的距離

for?i=1:50
????dd=d（1:i）;
????h=0.0075*dd.^（-3.6）;
????power0=5e-15*ones（1i）;
????ponce=power0;
????
j=sum（ponce.*h）-ponce.*h+q;
r=w/v*ponce.*h./j;
pnext=sir_tar*（ponce./r）-k*（ponce./r）.^2;

while?sum（abs（pnext-ponce）>=eps）>0
????ponce=pnext;
????j=sum（ponce.*h）-ponce.*h+q;
????r=w/v*ponce.*h./j;
????pnext=sir_tar*（ponce./r）-k*（ponce./r）.^2;
????[a]=find（pnext<0）