I am trying to simulate a paper . Please have a look to the problem (13) in the paper.

http://www.sce.carleton.ca/faculty/yanikomeroglu/Pub/ICC2016-ibyaekhy.pdf.

Is the following formulation correct?

x_max=3000;

y_max=3000;

h_min=100;

h_max=1000;

M1 = 1e4;

M2 = 1e4;

cvx_begin quiet

cvx_precision medium

cvx_solver mosek

variable U(1,NUM_USER) binary; %users

variable xd;

variable yd;

variable h

variable R;

expression d(NUM_USER,1);

expression th(NUM_USER,1);

expression rmin(NUM_USER,1);

maximize sum(U)

subject to

for i = 1:NUM_USER

d(i,1)=R^2;

th(i,1) = Gamma(i)+(M1*(1-U(i)));

d(i,1) <= th(i,1);

```
rmin(i) = norm([ (Uloc(i,1) - xd); (Uloc(i,2) - yd)]) - (M2* (1-U(i)));
R >= rmin(i);
```

end

xd^2 <= x_max;

yd^2 <= y_max;

h >= h_min;

h <= h_max;

R >= 0;

cvx_end