  In this article we introduce a simple straightforward and powerful method
involving symbolic manipulation, Picard iteration, and auxiliary variables for
approximating solutions of partial differential boundary value problems. The
method is easy to implement, computationally efficient, and it is highly
accurate. The output of the method is a function that approximates the exact
solution.
