2 years, 91 days

Maple 18

## HELP ME TO SOLVE THIS ERROR (USING DSOLV...

Maple 18

restart; with(plots);
[animate, animate3d, animatecurve, arrow, changecoords,

complexplot, complexplot3d, conformal, conformal3d,

contourplot, contourplot3d, coordplot, coordplot3d,

densityplot, display, dualaxisplot, fieldplot, fieldplot3d,

interactive, interactiveparams, intersectplot, listcontplot,

listcontplot3d, listdensityplot, listplot, listplot3d,

loglogplot, logplot, matrixplot, multiple, odeplot, pareto,

plotcompare, pointplot, pointplot3d, polarplot, polygonplot,

polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus,

semilogplot, setcolors, setoptions, setoptions3d, spacecurve,

sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]

fixedparameter1 := [n = .3, W[e] = .3, M = .2, gamma = 1, delta = -1, N[r] = .8, Pr = .72, Nb = .5, Nt = .5, Bi = 2, Pr = .72, Le = 5];
[n = 0.3, W[e] = 0.3, M = 0.2, gamma = 1, delta = -1, N[r] = 0.8,

Pr = 0.72, Nb = 0.5, Nt = 0.5, Bi = 2, Pr = 0.72, Le = 5]

eq1 := (1-n)*(diff(f(eta), eta, eta, eta))+f(eta)*(diff(f(eta), eta, eta))-M*(diff(f(eta), eta))+n*W[e]*(diff(f(eta), eta, eta, eta))*(diff(f(eta), eta, eta)) = 0;
/  d   /  d   /  d         \\\
(1 - n) |----- |----- |----- f(eta)|||
\ deta \ deta \ deta       ///

/  d   /  d         \\     /  d         \
+ f(eta) |----- |----- f(eta)|| - M |----- f(eta)|
\ deta \ deta       //     \ deta       /

/  d   /  d   /  d         \\\ /  d   /  d         \\
+ n W[e] |----- |----- |----- f(eta)||| |----- |----- f(eta)|| =
\ deta \ deta \ deta       /// \ deta \ deta       //

0
deq1; eval(eq1, fixedparameter1);
/  d   /  d   /  d         \\\
0.7 |----- |----- |----- f(eta)|||
\ deta \ deta \ deta       ///

/  d   /  d         \\       /  d         \
+ f(eta) |----- |----- f(eta)|| - 0.2 |----- f(eta)|
\ deta \ deta       //       \ deta       /

/  d   /  d   /  d         \\\ /  d   /  d         \\
+ 0.09 |----- |----- |----- f(eta)||| |----- |----- f(eta)|| =
\ deta \ deta \ deta       /// \ deta \ deta       //

0
eq2 := (1+(4/3)*N[r])*(diff(theta(eta), eta, eta))+Pr*f(eta)*(diff(theta(eta), eta))+Nb*(diff(phi(eta), eta))*(diff(theta(eta), eta))+Nt*(diff(theta(eta), eta))*(diff(theta(eta), eta)) = 0;
/    4     \ /  d   /  d             \\
|1 + - N[r]| |----- |----- theta(eta)||
\    3     / \ deta \ deta           //

/  d             \
+ Pr f(eta) |----- theta(eta)|
\ deta           /

/  d           \ /  d             \
+ Nb |----- phi(eta)| |----- theta(eta)|
\ deta         / \ deta           /

2
/  d             \
+ Nt |----- theta(eta)|  = 0
\ deta           /
deq2; eval(eq2, fixedparameter1);
/  d   /  d             \\
2.066666667 |----- |----- theta(eta)||
\ deta \ deta           //

/  d             \
+ 0.72 f(eta) |----- theta(eta)|
\ deta           /

/  d           \ /  d             \
+ 0.5 |----- phi(eta)| |----- theta(eta)|
\ deta         / \ deta           /

2
/  d             \
+ 0.5 |----- theta(eta)|  = 0
\ deta           /
eq3 := diff(phi(eta), eta, eta)+Pr*Le*f(eta)*(diff(phi(eta), eta))+Nt*(diff(theta(eta), eta, eta))/Nb = 0;
/  d   /  d           \\                /  d           \
|----- |----- phi(eta)|| + Pr Le f(eta) |----- phi(eta)|
\ deta \ deta         //                \ deta         /

/  d   /  d             \\
Nt |----- |----- theta(eta)||
\ deta \ deta           //
+ ----------------------------- = 0
Nb
deq3 := eval(eq3, fixedparameter1);
/  d   /  d           \\               /  d           \
|----- |----- phi(eta)|| + 3.60 f(eta) |----- phi(eta)|
\ deta \ deta         //               \ deta         /

/  d   /  d             \\
+ 1.000000000 |----- |----- theta(eta)|| = 0
\ deta \ deta           //
bcs1 := f(0) = 0, D(f)(0) = 1+gamma*(D@D)(F)(0)+delta*(D@D@D)(f)(0), D(f)(8) = 0;
f(0) = 0,

D(f)(0) = 1 + gamma @@(D, 2)(F)(0) + delta @@(D, 3)(f)(0),

D(f)(8) = 0
bc1 := eval(bcs1, fixedparameter1);
f(0) = 0, D(f)(0) = 1 + @@(D, 2)(F)(0) - @@(D, 3)(f)(0),

D(f)(8) = 0
bcs2 := D(theta)(0) = Bi*(theta(0)-1), theta(8) = 0;
D(theta)(0) = Bi (theta(0) - 1), theta(8) = 0
bc2 := eval(bcs2, fixedparameter1);
D(theta)(0) = 2 theta(0) - 2, theta(8) = 0
bcs3 := Nb*D(phi)(0)+Nt*D(theta)(0) = 0, Nb*D(phi)(0)+Nt*D(theta)(0) = 0, phi(8) = 0;
Nb D(phi)(0) + Nt D(theta)(0) = 0,

Nb D(phi)(0) + Nt D(theta)(0) = 0, phi(8) = 0
bc3 := eval(bcs3, fixedparameter1);
0.5 D(phi)(0) + 0.5 D(theta)(0) = 0,

0.5 D(phi)(0) + 0.5 D(theta)(0) = 0, phi(8) = 0
R := dsolve({bc1, bc2, bc3, deq1, deq2, deq3}, [f(eta), theta(eta), phi(eta)], numeric, output = listprocedure);
Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations

 Page 1 of 1
﻿