“…Let us sub-divide the domain [0, 2] into n uniform meshes of size h, such that D N s = {0 = s 0 , s 1 , ..., s n/2 = 1, s n/2+1 , ..., s n = 2, s i = s 0 + ih, i = 0(1)n, h = 2/n}. According to the procedures in [36], we consider a constant coefficient sub-equation from Equation (14) as as the approximation of W j+1 (s) at the grid point s i , we have…”