inner = 8; outer = 12; path_loss = 2; Ki_min = 0; Ki_max = 2; K(1) = 1; power(1) = 10; For i>1: K(i) = K(i-1) + DOI * randn(); K(i) = max(min(K(i),Ki_max),Ki_min); %Restrict to upper&lower power(i) = inner + path_loss * K(i);