model obstacle.mod; if $presolve == '0' then { fix{i in {0,M+1}, j in X} v[i,j]; fix{i in Y, j in {0,N+1}} v[i,j]; } solve; display max{i in 1.._nccons} abs(_ccon[i]), min{i in 1.._ncons} _con[i].slack;