**> restart:with(plots):blt:=7:**

**> lambda:=2:m:=3:s:=1:**

**> Eq1:=(diff(f(eta),eta,eta,eta))+(f(eta)*(diff(f(eta),eta,eta)))-((diff(f(eta),eta))^2)+lambda*(((f(eta)*(diff(f(eta),eta,eta,eta))))-2*(diff(f(eta),eta))*(diff(f(eta),eta,eta,eta))^2)-(M/(1+m^2))*((diff(f(eta),eta)+ms))=0;**

**> **

**> Eq2:=(diff(h(eta),eta,eta))+(f(eta)*(diff(h(eta),eta)))-((diff(f(eta),eta))*(h(eta)))+lambda*(((f(eta)*(diff(h(eta),eta,eta,eta))))+(h(eta)*(diff(f(eta),eta,eta,eta)))+(diff(f(eta),eta,eta))*(diff(h(eta),eta))-2*(diff(f(eta),eta))*(diff(h(eta),eta,eta))+(M/(1+m^2)))*(m*(diff(f(eta),eta)-h))=0;**

**> **

**> Eq3:=((f(eta))*(diff(theta(eta),eta)))+Pr*((diff(theta(eta),eta,eta)))=0;**

**> **

**> bcs1 := f(0) = 0, (D(f))(0) = 1, (h)(0) = 0, (theta)(0) = 1, (D(f))(blt) = 0, h(blt) = 0, theta(blt) = 0, D(D(f))(blt) = 0, (D(h))(blt)=0;**

**> **

**> L := [1,2,3];**

**> for k to 3 do R := dsolve(eval({Eq1, Eq2, Eq3, bcs1}, M = L[k]), [f(eta),h(eta),theta(eta)], numeric, maxmesh=10000, output = listprocedure);;Y || k := rhs(R[3]); YA || k := rhs(R[6]);YB || k := rhs(R[5]);YC || k := rhs(R[4]);YD || k := rhs(R[7]);end do:**

**> **

**> **

**> **

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