i have problem expressed in mathprog doesn't seem describable using c api. specifically, have constraints between variables. i've generated mathprog file , passed glpk. finds correct solution, don't see how access solution programmatically. returned glp_prob struct has no rows or columns. parse solution printed solver, i'm hoping there's better way. alternatively, if it's possible express constraints between variables using c api, suspect solve problem. mathprog code below. param t := 200; set b, dimen 3; set c, dimen 2; set s, dimen 2; set q := setof{(i,j,c) in b : j == 1} i; set := setof{(i,s) in s} i; set e := setof{(i,j) in q cross i} (i, j); var x{(i,j) in e}, >=0, <=1, binary; var y{i}, >=0, <=1, binary; maximize obj : sum{(i,j,c) in b} x[i,j] * c; s.t. q1c: sum{(i,s) in s} x[1,i] <= 1; s.t. q2c: sum{(i,s) in s} x[2,i] <= 1; s.t. size : sum{(i,c) in c} c * y[i] <= t; s.t. c111avail : x[1,1] <= y[1]; s.t. ...