To obtain second order accuracy in SPH scheme, sufficient conditions for smooth functions are discussed and formula are presented. Uniqueness is analyzed for solution of the second order particle approximation. SPH approximate solutions are numerically validated to arrive at nearly second accuracy in the case of periodic boundary conditions. While particle number increases, L₁-error decreases and L₁-order approaches to the second order particle approximation of SPH scheme.