Problem with solving DE

Ask a Question

Create a Post

Question:Problem with solving DE

Posted: diandra 5 Product: Maple

February 11 2016

>	$Shootlib := "C:\\diana/":$

(1)

$ODE := {A*(diff(j(eta), eta))+f(eta)*j(eta)-h(eta)*v(eta) = 0, diff(t(eta), eta)+Pr*Le*t(eta)+Nt*(diff(q(eta), eta))/Nb = 0, diff(q(eta), eta)+Pr*f(eta)*q(eta)+Nb*q(eta)*t(eta)+Nt*t(eta)^2 = 0, diff(y(eta), eta)+Sb*f(eta)*y(eta)-Pe*(w(eta), eta)*(diff(t(eta), eta))-Pe*t(eta)*y(eta) = 0, diff(v(eta), eta)+f(eta)*v(eta)-u(eta)^2+M*i(eta)^2-M*h(eta)*j(eta)-M+1 = 0, diff(f(eta), eta) = u(eta), diff(h(eta), eta) = i(eta), diff(i(eta), eta) = j(eta), diff(p(eta), eta) = q(eta), diff(r(eta), eta) = t(eta), diff(u(eta), eta) = v(eta), diff(w(eta), eta) = y(eta)};$

${.1*(diff(j(eta), eta))+f(eta)*j(eta)-h(eta)*v(eta) = 0, diff(t(eta), eta)+12.4*t(eta)+5*(diff(q(eta), eta)) = 0, diff(q(eta), eta)+6.2*f(eta)*q(eta)+2*q(eta)*t(eta)+10*t(eta)^2 = 0, diff(y(eta), eta)+2*f(eta)*y(eta)-.1*(w(eta), eta)*(diff(t(eta), eta))-.1*t(eta)*y(eta) = 0, diff(v(eta), eta)+f(eta)*v(eta)-u(eta)^2+.1*i(eta)^2-.1*h(eta)*j(eta)+.9 = 0, diff(f(eta), eta) = u(eta), diff(h(eta), eta) = i(eta), diff(i(eta), eta) = j(eta), diff(p(eta), eta) = q(eta), diff(r(eta), eta) = t(eta), diff(u(eta), eta) = v(eta), diff(w(eta), eta) = y(eta)}$

(2)

IC := {Pr*f(0)+m*q(0) = 0, f(0) = s*q(0)/(Pr*Le), h(0) = 0, i(etainf) = 1, j(0) = 0, p(etainf) = 1, q(0) = beta, r(etainf) = 0, t(0) = gamma, t(0) = -Nd*(1-v(0)), u(0) = b*v(0), u(etainf) = 1, v(0) = alpha, w(0) = 1, w(etainf) = 0, y(0) = delta};