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