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