Adomain decomposition method is a powerful technique to solve nonlinear differential equations and nonlinear Integro-differential equations. This paper presents a novel way of computing Adomian polynomials for nonlinear terms occur in nonlinear differential equations through partial Bell polynomial function(subroutine) in SageMath software. The obtained Adomian polynomials is used to solve nonlinear integral equations and nonlinear Integro-differential equations by Laplace Adomian decomposition method through symbolic computation in SageMath software. © 2020 IAEME Publication. All rights reserved.