using higher-order discretization for second-order derivative