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

how i can dsolve this equations and  convert the equations to an explicit first-order system ?

Also parameter 'p' is unknown and shoud be found.

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

``

restart

with(PDEtools, casesplit, declare);L:=0.75e-4:

[casesplit, declare]

(1)

eq1:=0.2079268293e-3*(diff(psi(x), x, x))-101.2059621*(diff(u(x), x, x))-101.2059621*(diff(w(x), x))*(diff(w(x), x, x));
##
eq2:=0.2079248162e-3+29.57317072*psi(x)-29.57317072*(diff(w(x), x))-1.996189024*10^(-9)*(diff(psi(x), x, x))-1.996189024*10^(-9)*(diff(w(x), x, x, x));
##
eq3:=151.8089432*(diff(w(x), x))^2*(diff(w(x), x, x))+101.2059621*(diff(u(x), x, x))*(diff(w(x), x))+101.2059621*(diff(u(x), x))*(diff(w(x), x, x))-0.2079268293e-3*(diff(psi(x), x, x))*(diff(w(x), x))-0.2079268293e-3*(diff(w(x), x, x))*(diff(psi(x), x))+29.57317070*(diff(w(x), x, x))-29.57317072*(diff(psi(x), x))-1.996189024*10^(-9)*(diff(psi(x), x, x, x))-1.996189024*10^(-9)*(diff(w(x), x, x, x, x))-p*(diff(w(x), x, x));

0.2079268293e-3*(diff(diff(psi(x), x), x))-101.2059621*(diff(diff(u(x), x), x))-101.2059621*(diff(w(x), x))*(diff(diff(w(x), x), x))

 

0.2079248162e-3+29.57317072*psi(x)-29.57317072*(diff(w(x), x))-0.1996189024e-8*(diff(diff(psi(x), x), x))-0.1996189024e-8*(diff(diff(diff(w(x), x), x), x))

 

151.8089432*(diff(w(x), x))^2*(diff(diff(w(x), x), x))+101.2059621*(diff(diff(u(x), x), x))*(diff(w(x), x))+101.2059621*(diff(u(x), x))*(diff(diff(w(x), x), x))-0.2079268293e-3*(diff(diff(psi(x), x), x))*(diff(w(x), x))-0.2079268293e-3*(diff(diff(w(x), x), x))*(diff(psi(x), x))+29.57317070*(diff(diff(w(x), x), x))-29.57317072*(diff(psi(x), x))-0.1996189024e-8*(diff(diff(diff(psi(x), x), x), x))-0.1996189024e-8*(diff(diff(diff(diff(w(x), x), x), x), x))-p*(diff(diff(w(x), x), x))

(2)

dsys3 := {eq1, eq2, eq3, psi(0) = 0, psi(L) = 0, u(0) = 0, u(L) = 0, w(0) = 0, w(L) = 0, ((D@@1)(psi))(0) = 0, ((D@@1)(w))(0) = 0, ((D@@1)(w))(L) = 0}; dsol5 := dsolve(dsys3, 'maxmesh' = 1200, numeric, abserr = .1, output = array([L]))

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

 

L

0.75e-4

(3)

``


 

Download error2.mw

 

 

Please Wait...