A finite volume scheme preserving the invariant region property for the coupled system of FitzHugh-Nagumo equations on distorted meshes. (1st July 2022)