clear; close all; clc;

load data_problem_1.mat;

n      = length(c_bar);
n_iter = 1000;


% LP
cvx_quiet true;
cvx_begin

% your solution

cvx_end

disp(x_LP);


% QP
% run s5p1d1


% comparison between LP and QP
cost_LP = zeros(n_iter, 1);

for k = 1 : n_iter
    
    c = c_bar + sqrtm(Sigma) * randn(n,1);
    
    % expected cost with LP design
    cost_LP(k,1) = c' * x_LP;
    
%     run s5p1d2
    
end

% run s5p1e