i did a solution of this ODE equation but the solution of paper is different from mine also in other some equation i have same problem i can't get exactly and pretty solution

how  define a function for computing multi-variable adomian polynomial  what is wrong with this? what i did mistake