x0{1} = zeros(2,2); x0{2} = zeros(1,2);
option.method = 'ell'; option.Nout = 500;
option.pstar = 0.95; option.delta = 0.001;
option.R = 1000*eye(5);
Sol = feaslmi(@ex3_robust, x0, option);
X = double(Sol{1})
F = double(Sol{2})
K = F*inv(X)