Dear all,

I am trying to solve the following system of equations by using dsolve, but I get the error: error, (in RootOf) expression independent of, _Z, could you please help me to solve it. Thank you.

restart;

Digits := 20;

with(plots);

Nr := .1; Nb := .3; Nt := .1; Rb := 0; Lb := 1; Le := 10; Pe := 1; ss := .2; aa := .1; bb := .2; cc := .3; nn := 1.5;

Eq1 := nn.(diff(f(eta), eta))^(nn-1).(diff(f(eta), `$`(eta, 2)))-(nn+1)/(2.*nn+1).eta.(diff(theta(eta), eta)-Nr.(diff(h(eta), eta))-Rb.(diff(g(eta), eta))) = 0;

Eq2 := diff(theta(eta), `$`(eta, 2))+nn/(2.*nn+1).f(eta).(diff(theta(eta), eta))+Nb.(diff(theta(eta), eta)).(diff(h(eta), eta))+Nt.((diff(theta(eta), eta))^2) = 0;

Eq3 := diff(h(eta), `$`(eta, 2))+nn/(2.*nn+1).Le.f(eta).(diff(h(eta), eta))+Nt/Nb.(diff(theta(eta), `$`(eta, 2))) = 0;

Eq4 := diff(g(eta), `$`(eta, 2))+nn/(2.*nn+1).Lb.f(eta).(diff(g(eta), eta))-Pe.((diff(g(eta), eta)).(diff(h(eta), eta))+(diff(h(eta), `$`(eta, 2))).g(eta)) = 0;

etainf := 10;

bcs := f(0) = ss/Le.(D(h))(0), theta(0) = lambda+aa.(D(theta))(0), h(0) = lambda+bb.(D(h))(0), g(0) = lambda+cc.(D(g))(0), (D(f))(etainf) = 0, theta(etainf) = 0, h(etainf) = 0, g(etainf) = 0;

dsys := {Eq1, Eq2, Eq3, Eq4, bcs};

dsol := dsolve(dsys, numeric, continuation = lambda, output = procedurelist);

Error, (in RootOf) expression independent of, _Z