The homotopy analysis method for solving higher dimensional initial boundary value problems of variable coefficients