Hello, I am trying to solve the following BVP,

Eq1 := -mu0*M*(diff(theta(eta), eta))*exp(-M*theta(eta))*abs(diff(f(eta), eta))^m*(diff(f(eta), eta))

+mu0*exp(-M*theta(eta))*(m+1)*abs(diff(f(eta), eta))^m*(diff(f(eta), eta, eta))+theta(eta) = 0;

Eq2 := G*(diff(theta(eta), eta, eta))+mu0*exp(-M*theta(eta))*abs(diff(f(eta), eta))^(m+2) = 0;

bc := f(-1) = 0, f(1) = 0, theta(-1) =0.5, theta(1) = -0.5;

m:=-0.3 (or 0.3) :M:=1:mu0:=0.5:G:=3: