To solve (partial) differential equations it is necessary to have good numerical approximations. In SPH, most approximations suffer from the presence of boundaries. In this work a new approximation for the second-order derivative is derived and numerically compared with two other approximation methods for a simple test case. The new method is slightly more expensive, but leads to a significantly improved accuracy
