Hello every one,

I am trying to solve a system of 5 odes numerically but I get an error

"unable to convert to an explicit first-order system".

Please have a look and guide me through the problem.

Thanks

wprob1.mw

Is it possible to have a contour plot of an ode, 

restart:with(plots);alpha:=1:

de1:=alpha*diff(f(y),y$4)+(diff(f(y),y$2)*diff(f(y),y$3)^2+diff(f(y),y$4)*diff(f(y),y$2)^2);

with  such type of bcs?

bc1:= f(1+sin(x)) = .95, f(-.1-.1*sin(x)) = -.95, (D(f))(1+sin(x)) = -1, (D(f))(-.1-.1*sin(x)) = -1;

contourplot({de1,bc1}, [x,y, f(y)],x=1..10, y=1..10);

Thanks

eq1:=diff(f(eta),eta,eta,eta)+f(eta)*diff(f(eta),eta,eta)+beta*(1-diff(f(eta),eta)^2)=0;

eq2:=diff(theta(eta), eta, eta)+pr*f(eta)*(diff(theta(eta), eta))+pr*Ec*(diff(f(eta), eta, eta))^2

+(2-beta)*pr*Q*theta(eta) = 0;

bc:=theta(9)=0,D(theta)(0)=-Nc*(1-theta(0)),f(0)=1,D(f)(0)=lambda,D(f)(9)=1};

Para:={beta=1,Ec=0.5,Q=-0.5,Nc=5,pr=3}:

How to plot -D(theta)(0) vs lambda?

Thanks

Hi every one,

I am facing a situation here. Whenever, I sub n = 0 or 1, the following ode has a solution. 

But when I adopt other values except the above ones then the maple gives me the error

"division by zero". 

eq11:=(diff(((abs(diff(f(eta),eta,eta)))^(n-1)*diff(f(eta),eta,eta)),eta)

+n/(n+1)*f(eta)*diff(f(eta),eta,eta))=0;

resx:=dsolve({eq11,f(0)=1,D(f)(0)=0,D(f)(12)=0},numeric,method = bvp[midrich],

I am looking for a numerical scheme which can easly handle such type of problem.

I already had the analytical results of this problem, now I want to know how to treat this problem numerically.

Any idea, then please let me know.

Thanks