with plots

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

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

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

BCs := [f(0), (D(f))(0), (D(f))(18), g(0)-1, g(18), h(18)]

pars := {Pr = .1, m = 0}

for i to 3 do eq[i] := subs(pars, Eq[i]) end do