Q_m = struct('name','Q_m','val', [4]);
Q_p = struct('name','Q_p','val',[0]);
R_m = struct('name','R_m','val',[4]);
R_p = struct('name','R_p','val',[0]);
t = struct('name','t','val','s1');

wgdx('gams_matlab.gdx', Q_m, Q_p, R_m, R_p, t);

cost = gams('myopic_gams_matlab.gms LP=cplex', t, Q_m, Q_p, R_m, R_p);
reward = cost.val;Q_m = struct('name','Q_m','val', [4]);
Q_p = struct('name','Q_p','val',[0]);
R_m = struct('name','R_m','val',[4]);
R_p = struct('name','R_p','val',[0]);
t = struct('name','t','val','s1');

wgdx('gams_matlab.gdx', Q_m, Q_p, R_m, R_p, t);

cost = gams('myopic_gams_matlab.gms LP=cplex', t, Q_m, Q_p, R_m, R_p);
reward = cost.val;