Recently several fast computational methods have been developed for solving finite difference approximations to the standard partial differential equations of mathematical physics. In this paper, the application of a new fast algorithm for solving tridiagonal systems with constant coefficients occurring in the solution of boundary value problems is presented. It is shown that the development of such methods of solution for commonly occurring problems in mathematical and computational physics can yield fruitful gains in efficiency.