(Subscript[y, j])={(1,if   reaction   flux   (Subscript[v, j])   is   active        0,if   reaction   flux   (Subscript[v, j])   is   not   active),∀j∈M