xvalues = [list all possible x values]
yvalues = [list all possible y values]
zvalues = [list all possible z values];
[X, Y, Z] = ndgrid(xvalues, yvalues, zvalues);
mask = X + Y + Z <=1;
x = X(mask);
y = Y(mask);
z = Z(mask);
f = x.^3 - log(y+z) + sin(z);
[bestf, bestidx] = min(f);
bestx = x(bestidx);
besty = y(bestidx);
bestz = z(bestidx);
Best Answer