Advanced Search

Question: How do I solve Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

Question:How do I solve Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

jennierubyjane 25 Maple
ode syntax homework dsolve numeric differential-equations bvp
December 19 2021
2

> 

restart
> 

``
> 

with(plots)

 

> 

fixedparameter := [Nb = 0, Nt = 0, Bi = 1000, Le = 10]

[Nb = 0, Nt = 0, Bi = 1000, Le = 10]

 (1)
> 

DE1 := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))-[diff(f(eta), eta)]^2 = 0

diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))-[diff(f(eta), eta)]^2 = 0

 (2)
> 

``
> 

DE2 := diff(theta(eta), eta, eta)+Pr*f(eta)*(diff(theta(eta), eta))+Pr*Nb*(diff(theta(eta), eta))*(diff(phi(eta), eta))+Pr*Nt*[diff(theta(eta), eta)]^2 = 0

diff(diff(theta(eta), eta), eta)+Pr*f(eta)*(diff(theta(eta), eta))+Pr*Nb*(diff(theta(eta), eta))*(diff(phi(eta), eta))+Pr*Nt*[diff(theta(eta), eta)]^2 = 0

 (3)
> 

DE3 := diff(phi(eta), eta, eta)+Le*(diff(phi(eta), eta))+Nt*(diff(theta(eta), eta, eta))/Nb = 0

diff(diff(phi(eta), eta), eta)+Le*(diff(phi(eta), eta))+Nt*(diff(diff(theta(eta), eta), eta))/Nb = 0

 (4)
> 

BC1 := f(0) = 0, f(eta) = 1-exp(-eta), (D(f))(0) = 1, (D(f))(10) = 0

f(0) = 0, f(eta) = 1-exp(-eta), (D(f))(0) = 1, (D(f))(10) = 0

 (5)
> 

BC2 := theta(10) = 0, (D(theta))(0) = -Bi*[1-theta(0)]

theta(10) = 0, (D(theta))(0) = -Bi*[1-theta(0)]

 (6)
> 

BC3 := phi(0) = 1, phi(10) = 0

phi(0) = 1, phi(10) = 0

 (7)
> 

``
> 

L := [0.7e-1, .2, .7, 2, 7, 20, 70]; for k to 7 do R := dsolve(eval({BC1, BC2, BC3, DE1, DE2, DE3}, Pr = L[k]), [f(eta), theta(eta), phi(eta)], numeric, output = listprocedure); Y || k := rhs(R[5]); YP || k := -rhs(R[6]) end do

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

  
> 

 
> 

NULL

Download ODEprobz.mw
Please Wait...