Question: Collecting coefficient of perturbation equations

I have a system of three non-linear coupled ODEs. I am trying to obtain the zeroth order, first order and second order equations and their boundary conditions using the perturbation technique. I then need to solve them. I wrote the following MAPLE lines for the purpose.

restart; N := 2

Eq[1] := diff(f(eta), `$`(eta, 3))+(m+4)*f(eta)*(diff(f(eta), `$`(eta, 2)))-(2*(m+1))*(diff(f(eta), eta))^2-(2*(m+1))*h(eta)-(.5*(m-2))*eta*(diff(h(eta), eta))

Eq[2] := g(eta)-(diff(h(eta), eta))

Eq[3] := (diff(g(eta), `$`(eta, 2)))/Pr+(m+4)*f(eta)*(diff(g(eta), eta))-(5*m+2)*g(eta)*(diff(f(eta), eta))

BCs := [f(0), (D(f))(0), (D(f))(infinity), (D(g))(0)+1, g(infinity), h(infinity)]

eqs := Eq[1], Eq[2], Eq[3]

f(eta):=sum('f[i](eta)*m^(i)','i' =0..N);

h(eta):=sum('h[i](eta)*m^(i)','i' =0..N);

g(eta):=sum('g[i](eta)*m^(i)','i' =0..N);

EQS := expand(eqs)

MAPLE says 

Error, (in sum) too many levels of recursion

How do I proceed further?perturbation_equatio.mwperturbation_equatio.mw