param d := 5;

param n := 100;

param a {1..d, 1..n}, default Uniform01();

var y >= 0;
var x {1..d} >= 0;

param i0 within {1..n};

maximize ydist: sum{j in 1..d} a[j,i0]*x[j] / y;

subject to one_sided {i in 1..n}: sum {j in 1..d} a[j,i]*x[j] <= y;

option solver loqo;
option loqo_options "pred_corr=0 mufactor=0 verbose=2";

for {i in 1..n} {
    let i0 := i;
    solve;
}