Symbolic Iterative Solution of Boundary Value Problems for Partial Differential Equations

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.