Adomian Decomposition

Ask a Question

Create a Post

Question:Adomian Decomposition

Posted: Gabriel samaila 35 Product: Maple

September 21 2019

Hi guys,

I'm trying to apply boundary condition on the seris i generated using Adomian decomposition. Could some one please check it for me?

A[1] := R*(diff(u[0], y))+A-Gr*(B*phi[0]+theta[0])/Ree+Ha^2*(alpha*u[0]+beta*w[0])/(alpha^2+beta^2)

R*a2+A-Gr*(a5+a6*y+B*(a7+a8*y))/Ree+Ha^2*(alpha*(a1+a2*y)+beta*(a3+a4*y))/(alpha^2+beta^2)

(1)

u[1] := int(A[1], y = 0 .. y)

(1/2)*(-Gr*(a6+B*a8)/Ree+Ha^2*(alpha*a2+beta*a4)/(alpha^2+beta^2))*y^2+R*a2*y+A*y-Gr*(a5+B*a7)*y/Ree+Ha^2*(alpha*a1+beta*a3)*y/(alpha^2+beta^2)

(2)

u[11] := int(u[1], y = 0 .. y)

(1/3)*(-(1/2)*Gr*(a6+B*a8)/Ree+(1/2)*Ha^2*(alpha*a2+beta*a4)/(alpha^2+beta^2))*y^3+(1/2)*(R*a2+A-Gr*(a5+B*a7)/Ree+Ha^2*(alpha*a1+beta*a3)/(alpha^2+beta^2))*y^2

(3)

u = a1+a2*y+(1/3)*(-(1/2)*Gr*(a6+B*a8)/Ree+(1/2)*Ha^2*(alpha*a2+beta*a4)/(alpha^2+beta^2))*y^3+(1/2)*(R*a2+A-Gr*(a5+B*a7)/Ree+Ha^2*(alpha*a1+beta*a3)/(alpha^2+beta^2))*y^2

(4)

A[2] := R*(diff(w[0], y))-Ha^2*(beta*u[0]-alpha*w[0])/(alpha^2+beta^2)

R*a4-Ha^2*(beta*(a1+a2*y)-alpha*(a3+a4*y))/(alpha^2+beta^2)

(5)

w[1] := int(A[2], y = 0 .. y)

-(1/2)*Ha^2*(beta*a2-alpha*a4)*y^2/(alpha^2+beta^2)+R*a4*y-Ha^2*(beta*a1-alpha*a3)*y/(alpha^2+beta^2)

(6)

w[11] := int(w[1], y = 0 .. y)

-(1/6)*Ha^2*(beta*a2-alpha*a4)*y^3/(alpha^2+beta^2)+(1/2)*(R*a4-Ha^2*(beta*a1-alpha*a3)/(alpha^2+beta^2))*y^2

(7)

w = a3+a4*y-(1/6)*Ha^2*(beta*a2-alpha*a4)*y^3/(alpha^2+beta^2)+(1/2)*(R*a4-Ha^2*(beta*a1-alpha*a3)/(alpha^2+beta^2))*y^2

(8)

A[3] := R*Pr*(diff(theta[0], y))-2*Br*((diff(u[0], y))^2+(diff(w[0], y))^2+M^2*C*(u[0]^2+w[0]^2))-gamma*R*Pr*theta[0]

R*Pr*a6-2*Br*(a2^2+a4^2+M^2*C*((a1+a2*y)^2+(a3+a4*y)^2))-gamma*R*Pr*(a5+a6*y)

(9)

theta[1] := int(A[3], y = 0 .. y)

-(2/3)*Br*M^2*C*(a2^2+a4^2)*y^3+(1/2)*(-2*Br*M^2*C*(2*a1*a2+2*a3*a4)-gamma*R*Pr*a6)*y^2+R*Pr*a6*y-2*Br*(a2^2+a4^2+M^2*C*(a1^2+a3^2))*y-gamma*R*Pr*a5*y

(10)

theta[11] := int(theta[1], y = 0 .. y)

-(1/6)*Br*M^2*C*(a2^2+a4^2)*y^4+(1/3)*(-Br*M^2*C*(2*a1*a2+2*a3*a4)-(1/2)*gamma*R*Pr*a6)*y^3+(1/2)*(R*Pr*a6-2*Br*(a2^2+a4^2+M^2*C*(a1^2+a3^2))-gamma*R*Pr*a5)*y^2

(11)

theta = a5+a6*y-(1/6)*Br*M^2*C*(a2^2+a4^2)*y^4+(1/3)*(-Br*M^2*C*(2*a1*a2+2*a3*a4)-(1/2)*gamma*R*Pr*a6)*y^3+(1/2)*(R*Pr*a6-2*Br*(a2^2+a4^2+M^2*C*(a1^2+a3^2))-gamma*R*Pr*a5)*y^2