x = optimvar('x');
y = optimvar('y');
z = optimvar('z');
prob = optimproblem;
H=[146 187 145; 187 854 104; 145 104 289];
prob.Objective = [x y z]*H*[x;y;z];
prob.Constraints.cons1 = x+y+z==1;
prob.Constraints.cons2 = 0.175*x+0.24*y+0.227*z>=0.2;

sol = solve(prob)