In this article, a modern analytical technique called Homotopy Analysis Method (HAM) has been applied to several partial differential equations to obtain their sequential solution. The performance of the method was analysed on partial differential equations and found to be computationally efficient. The homotopy analysis method (HAM) has been found to be effective. The numerical results have been improved by proposing a new technique that combines the numerical method with the Harris Hawks Optimization algorithm. The current numerical results are compared to the exact solutions and it's found to be in very good results.