Advanced Search

Question: Error, (in dsolve/numeric/bvp) division by zero

Question:Error, (in dsolve/numeric/bvp) division by zero

J4James 355 Maple 16
numeric
May 30 2013
1

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:

dsolve({bc, Eq1, Eq2},type=numeric,method=bvp[midrich]): 

Error, (in dsolve/numeric/bvp) division by zero

I need to solve the above system for m=0.3 and -0.3 but I do not know how to make it work.

Please help me out.

Thanks
Please Wait...