We consider the numerical solution of a third-order Falkner-Skan-like boundary value problem arising in boundary layer theory. The problem is defined on a semi-infinite interval [0,∞) with a condition given at ∞. We first transform the problem into a second-order boundary value problem defined on a finite interval [γ,1]. To solve the resulting boundary value problem, we developed an iterative finite-difference scheme based on Newton's quasilinearization. At every step, the linearized differential equation is approximated using the finite-difference method. Numerical results will be presented to demonstrate the efficiency of the method and will be compared with other results presented in the literature.