Question: solving a differential equation

How can I solve the following differential equation?

with plots

Eq[1]:=diff(f(eta),eta$3)+g(eta);

Eq[2] := (diff(g(eta), `$`(eta, 2)))/Pr+f(eta)*(diff(g(eta), eta));

BCs := [f(0), (D(f))(0), (diff(f(eta), `$`(eta, 2)))(2), g(0)-1, g(2)]

pars := {Pr = 1}

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

eqs := eq[1], eq[2]

vars := f(eta), g(eta)

bcs := op(subs(pars, BCs))

sol := dsolve({bcs, eqs}, {vars}, type = numeric, output = array([seq(i*1, i = 0 .. 2)]))

 

Maple is saying: Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

Kindly help me.

Thanking you 

Please Wait...