Question: Problem for solving system of ODEs

Dear friends! I am facing problem to solve the below system of ODEs numerically please find the mistake and correct it.

alpha := -1; R := 2; m := 2; Pr := 7; Le := 1.25; Nt := .2; Nb := .2; g := .5; K1 := .1; Q := .5

Eq1 := eta^3*(diff(F(eta), eta, eta, eta, eta))+alpha*(eta^4*(diff(F(eta), eta, eta, eta))+eta^3*(diff(F(eta), eta, eta))-eta^2*F(eta))-2*eta^2*(diff(F(eta), eta, eta, eta))+3*eta*(diff(F(eta), eta, eta))-3*(diff(F(eta), eta))+eta*R*(diff(F(eta), eta))^2-3*eta*R*F(eta)*(diff(F(eta), eta, eta))+3*R*F(eta)*(diff(F(eta), eta))+3*eta^2*R*F(eta)*(diff(F(eta), eta, eta, eta))-eta^2*(diff(F(eta), eta))*(diff(F(eta), eta, eta))-M^2*(eta^3*(diff(F(eta), eta, eta))-eta^2*(diff(F(eta), eta))); Eq2 := eta*(diff(G(eta), eta, eta))+alpha*Pr*eta^2*(diff(G(eta), eta))+R*Pr*F(eta)*(diff(G(eta), eta))+Nb*eta*(diff(G(eta), eta))*(diff(H(eta), eta))+Nt*eta*(diff(G(eta), eta))^2+diff(G(eta), eta)+Q*Pr*eta*G(eta) = 0; Eq3 := eta*(diff(H(eta), eta, eta))+alpha*Le*Pr*eta^2*(diff(H(eta), eta))+R*Le*Pr*F(eta)*(diff(H(eta), eta))+Nt*eta*(diff(G(eta), eta, eta))/Nb+Nt*(diff(G(eta), eta))/Nb+diff(H(eta), eta)-g*Le*Pr*eta*H(eta)-Le*Pr*K1*eta = 0;

IC1 := F(0) = 0, F(1) = 1, (D(F))(0) = 0, (D(F))(1) = 0, (D(G))(0) = 0, G(1) = 1, (D(H))(0) = 0, H(1) = lambda; dsys1 := {Eq1, Eq2, Eq3, IC1}; dsol1 := dsolve(dsys1, numeric, continuation = lambda, range = 0 .. 1);

dsol1x := subs(dsol1, F(eta));

dsol1y := subs(dsol1, G(eta)); dsol1z := subs(dsol1, H(eta));
 

With my best regards and sincerely.

Muhammad Usman

School of Mathematical Sciences 
Peking University, Beijing, China