Is there any an alternative way to write this equation:
X = 2 * r1 * r2^(1/d) - r3^(1/d)
r1 r2 and r3 are uniform random numbers in interval (0,1)
d : it can be 10, 30 or 50
X is a vector [1,d],
For example:
X(1) = = 2 * r1 * r2^(1/d) - r3^(1/d)**
generate new r1, r2 r3
X(2) = = 2 * r1 * r2^(1/d) - r3^(1/d)**
generate new r1, r2 r3
X(3) = = 2 * r1 * r2^(1/d) - r3^(1/d)**
.
.
.
generate new r1, r2 r3
X(d) = = 2 * r1 * r2^(1/d) - r3^(1/d)**