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. 
Hamid Semiyari
math.GM