In this research, we propose fourth-order non-uniform Hermitian differencing with a fifth-order adaptive time integration method for pricing system of free boundary exotic power put options consisting of the option value, delta sensitivity, and gamma. The main objective for implementing the above scheme is to carefully account for the irregularity in the locality of the left corner point after fixing the free boundary. Specifically and mainly, we stretch the performance of our proposed method threefold. First, we exploit the non-uniform fourth-order Hermitian scheme to locally concentrate space grid points arbitrarily close to the left boundary. Secondly, we further leverage the adaptive nature of the embedded time integration method, which allows optimal selection of a time step based on the space grid point distribution and regional variation. Thirdly, we introduce a fourth-order combined Hermitian scheme, which requires fewer grid points for computing the near boundary point of the delta sensitivity and gamma. Another novelty is how we approximate the optimal exercise boundary and its derivative using a fifth-order Robin boundary scheme and fourth-order combined Hermitian scheme. Our proposed method consistently achieves reasonable accuracy with very coarse grids and little runtime across the numerical experiments. We further compare the results with existing methods and the ones we obtained from the uniform space grid.