A numerical algorithm is developed for the approximation of the solution to certain boundary value problems involving the third-order ordinary differential equation associated with draining and coating flows.The authors show that the approximate solutions obtained by the numerical algorithm developed by using nonpolynomial quintic spline functions are better than those produced by other spline and domain decomposition methods.The algorithm is tested on two problems associated with draining and coating flows to demonstrate the practical usefulness of the approach.