
index = indexes;
for i = 1:length(DensityRange)
    
    Capacity(i)=(W./NoOfBandsRange(index(i))).*log2(1+Ptopt(i,index(i)).*lambda1^alpha1./(InterfOptimal(i,index(i))));

    SINRdB(i) = 10*log10(Ptopt(i,index(i)).*lambda1^alpha1./(InterfOptimal(i,index(i))));
    
    SNRneeded(i) = 2^(2*Ractual*NoOfBandsRange(index(i))/W)-1;
    
    SNRneededdB(i) = 10*log10(SNRneeded(i));
end