“…Problem (1.1)-(1.3) has been studied in several analytical contributions [1,5,6,7,8,9,10,12,13,15,16,17]. The existence of weak solutions for ε > 0 has been verified in [8], where the implicit Euler method is used to semidiscretize the problem in time: Let I k = {t j } J j=1 be an equidistant net of mesh width k = t j − t j−1 to discretize [0, T ], with j ≥ 1.…”