r=[29 18 -5; 30 16 4; 6 13 10];
fun = @(x) -mean(x*r);  %
x0 = [0,0,0];       % Starting point
lb = [-1,-1,-1];    % Lower bound limit
ub = [1,1,1];       % Upper bound limit
A=[1 1 1];  b=[1];  % Inequality constraints, A*x <= 1
Aeq=[1 1 1];  beq=1;    % Linear equality constrints, Aeq*x = 1
x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub);