function[objval] = readmps(file);
%H1 read a model and solve it

model = mindopt_read(file);
result = mindopt(model);

for v=1:length(result.X)
    fprintf('%s %d\n', model.varnames{v}, result.X(v));
end
fprintf('Obj: %e\n', result.ObjVal);

objval = result.ObjVal;