“…In Table 3, we give the L ∞ and L 2 errors at T = 1 for 𝜏 = 0.0001,Ω = [−20, 20] × [ −20, 20] and increasing values of nodes where one can see that as the number of nodes increase more accurate results can be obtained. In Table 4, a comparison of the present method with the methods of [26], namely Crank-Nicolson finite difference finite element (ECFD-FE), semi-finite difference finite element (SIFD-FE), explicit finite difference finite element (EXFD-FE) and time-splitting Crank-Nicolson finite element (TSCN-FE), is given. In this case our grid consists of 36 × 36 nodes, whereas in [26] the spatial step size is h = 1∕32 which means they used a grid which consists of about 640 × 640 nodes.…”