Sol = scenlmi(@ex4_opt, 0.95, 0.001);
X = double(Sol{1})
F = double(Sol{2})
gamma = double(Sol{3})
K = F*inv(X)