Hi,

I have the following set of equations and boundary conditions. The problem is I cannot be able to make accurate guesses with the shooting method. I request your assistance. 

> k1 := diff(X[1](t), t) = X[2](t);
> k2 := diff(X[2](t), t) = M*(X[1](t)-1)-(2*(eta+b))*X[2](t);
> k3 := diff(X[3](t), t) = X[4](t);
> k4 := diff(X[4](t), t) = (2*Sc*Du*(eta+b)*X[6](t)-Du*lambda*X[5](t)-2*Pr*(eta+b)*X[4](t)-Pr*Ec*X[2](t)^2-Pr*Ec*M*(X[1](t)-1)^2)/(1-Du*Sr);