The piecewise quadratic polynomial collocation is used to approximate the nonlocal model, which generally leads to a nonsymmetric indefinite system (Chen et al. (2021) [5]). In this case, the discrete maximum principle is not satisfied, which might be trickier for the stability analysis of the high-order numerical schemes (D'Elia et al. (2020) [10]; Leng et al. (2021) [26]). Here, we present a modified (shifted-symmetric) piecewise quadratic polynomial collocation for solving the linear nonlocal diffusion model, which leads to a symmetric positive definite system and satisfies the discrete max...