We present a computational study of the solution of the Falkner-Skan equation (a thirdorder boundary value problem arising in boundary-layer theory) using high-order and high-order-compact finite differences schemes. There are a number of previously reported solution approaches that adopt a reduced-order system of equations, and numerical methods such as: shooting, Taylor series, Runge-Kutta and other semi-analytic methods. Interestingly, though, methods that solve the original non-reduced third-order equation directly are absent from the literature. Two high-order schemes are presented using both explicit (third-order) and implicit compact-difference (fourth-order) formulations on a semi-infinite domain; to our knowledge this is the first time that high-order finite difference schemes are presented to find numerical solutions to the non-reduced-order Falkner-Skan equation directly. This approach maintains the simplicity of Taylor-series coefficient matching methods, avoiding complicated numerical algorithms, and in turn presents valuable information about the numerical behaviour of the equation. The accuracy and effectiveness of this approach is established by comparison with published data for accelerating, constant and decelerating flows; excellent agreement is observed. In general, the numerical behaviour of formulations that seek an optimum physical domain size (for a given computational grid) is discussed. Based on new insight into such methods, an alternative optimisation procedure is proposed that should increase the range of initial seed points for which convergence can be achieved.