The application of finite differences methods to solve boundary value problems typically involves a discretization of such a problem across an orthogonal array of discrete grid points. This leads to an array of difference equations which is solved numerically within the constraints of the boundary conditions to yield solutions at the grid point locations. However, the accuracy of the solutions is limited with conventional finite differences methods when the boundary conditions are not represented exactly within the orthogonal array of discrete grid points, as when the boundary conditions are curved or slanted surfaces. The invention described herein provides finite differences methods for solving boundary value problems more accurately than with conventional finite differences methods, particularly when curved or slanted boundary surfaces correspond to terminations of a known analytical function.