I want to get the zeroth order and first order perturbation equations of the following nonlinear coupled ODE.
restart; N := 1
alias(eta = e, theta = t)
eq[1] := 5*(diff(F(e), `$`(e, 3)))+(1+3/m)*F(e)*(diff(F(e), `$`(e, 2)))-(2+1/m)*(diff(F(e), e))^2-(4+2/m)*H(e)-(1-2/m)*e*(diff(H(e), e)) = 0
eq[2] := diff(H(e), e) = t(e)
eq[3] := 5*(diff(t(e), `$`(e, 2)))/Pr+(1+3/m)*F(e)*(diff(t(e), e))-5*(diff(F(e), e))*t(e)
eqs := eq[1], eq[2], eq[3]