Note that a constraint is composed of a constraint expression and forall sets expression.
ConstraintExpr expr = capacity[i, j] >= density * x[i, j];
Constraint con = forall(i, j) | capacity[i, j] >= density * x[i, j];
public static ConstraintExpr operator >=(
Sca lhs,
Term rhs
)