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