“…If one puts, without loss of generality, ε = δ = 1, then one obtains the more general system 1 + 2g(ρ) ∂ t µ + µ g ′ (ρ) ∂ t ρ − ∆µ = 0, (1.17) ∂ t ρ − ∆ρ + W ′ (ρ) = µ g ′ (ρ), (1.18) which was investigated in the contributions [11,15,17,19], still for no-flux boundary conditions, also from the side of the numerical approximation. The related phase relaxation system (in which the diffusive term −∆ρ disappears from (1.18)), has been dealt with in [12,13,21]. We also mention the recent article [25], where a nonlocal version of (1.17)-(1.18) -based on the replacement of the diffusive term of (1.18) with a nonlocal operator acting on ρ -has been largely investigated.…”