set i /1*150/;
alias (i,j);

scalar thresh /0.5/;
parameter D(i,j);
option seed = 101;

D(i,j)$(ord(i) ge ord(j)) = uniform(0,1);
D(i,j)$(ord(i) lt ord(j)) = D(j,i);