Solution of the Falkner–Skan equation by recursive evaluation of Taylor coefficients

We present a computational method for the solution of the third-order boundary value problem characterized by the well-known Falkner–Skan equation on a semi-infinite domain. Numerical treatments of this problem reported in the literature thus far are based on shooting and finite differences. While maintaining the simplicity of the shooting approach, the method presented in this paper uses a technique known as automatic differentiation, which is neither numerical nor symbolic. Using automatic differentiation, a Taylor series solution is constructed for the initial value problems by calculating the Taylor coefficients recursively. The effectiveness of the method is illustrated by applying it successfully to various instances of the Falkner–Skan equation.