I tried to solve a non linear coupled boundary value problem in MAPLE using DSolve command. The code is :
alias(eta = e, theta = t)
Eq[1] := 5*(diff(F(e), `$`(e, 3)))+(m+3)*F(e)*(diff(F(e), `$`(e, 2)))-(2*m+1)*(diff(F(e), e))^2-(4*m+2)*H(e)-(m-2)*e*(diff(H(e), e)) = 0
Eq[2] := diff(H(e), e) = t(e)
Eq[3] := 5*(diff(t(e), `$`(e, 2)))/Pr-(m+3)*F(e)*(diff(t(e), e))-5*m*(diff(F(e), e))*t(e) = 0
BCs := [F(0), (D(F))(0), (D(F))(infinity), t(0)-1, t(infinity), H(infinity)]