Semianalytical method for the Graetz problem

This PDE is first order in z. The z variable can be treated as a time variable. The temperature profiles depend on the Peclet number Pe. For convenience we introduce the variable transformation z = 2Pet, which converts equation (5.29) to 

Equation (5.30) is solved below in Maple using the Maple program developed for example 5.2 by making very few changes as  

5 Method of Lines for Parabolic Partial Differential Equations 

---