Advanced Search

Question: error for dsolve differential equation....

Question:error for dsolve differential equation....

farahnaz 15 Maple 2016
ode differential-equation dsolve nonlinear numeric
October 12 2016
0

hi..i have a problem for solving this nonlinear differential equationerror.mw

restart; Digite := 200; L := 100*10^(-9)

1/10000000

 (1)

EQ11 := -3.000000000*10^(-8)+3.815358072*sin(3.141592654*10^7*x)+9.534375000*10^(-30)*(diff(w(x), x, x, x, x))-2.383593750*10^(-60)*(diff(w(x), x, x, x, x, x, x))-5.085000000*10^(-13)*(diff(w(x), x))*(diff(u(x), x, x))-7.627500000*10^(-13)*(diff(w(x), x))^2*(diff(w(x), x, x))-5.085000000*10^(-13)*(diff(w(x), x, x))*(diff(u(x), x))+0.2410290000e-5*(diff(w(x), x, x)):

EQ2 := 5.650000000*10^(-20)*(diff(u(x), x, x, x, x))-226000000000*(diff(u(x), x, x))-226000000000*(diff(w(x), x))*(diff(w(x), x, x)):

dsys3 := {EQ11, EQ2, u(0) = 0, u(L) = 0, w(0) = 0, w(L) = 0, (D(u))(0) = 0, (D(u))(L) = 0, ((D@@2)(w))(0) = 0, ((D@@2)(w))(L) = 0, ((D@@4)(w))(0) = 0, ((D@@4)(w))(L) = 0}; res := dsolve(dsys3, numeric, initmesh = 1024, abserr = 0.1e-4); res(0.1e-9)

res(0.1e-9)

 (2)

############################################################CHANGE OF VARIABLE:::           x=y*L

bcs := {u(0) = 0, u(L) = 0, w(0) = 0, w(L) = 0, (D(u))(0) = 0, (D(u))(L) = 0, ((D@@2)(w))(0) = 0, ((D@@2)(w))(L) = 0, ((D@@4)(w))(0) = 0, ((D@@4)(w))(L) = 0}; sys := {EQ11, EQ2}; sys2 := PDEtools:-dchange({x = L*y, u(x) = g2(y), w(x) = g1(y)}, sys, [g1, g2, y]); solve(sys2, {diff(g2(y), y, y, y, y), diff(g1(y), y, y, y, y, y, y)}); bcs3 := {g1(0) = 0, g1(1) = 0, g2(0) = 0, g2(1) = 0, (D(g2))(0) = 0, (D(g2))(1) = 0, ((D@@2)(g1))(0) = 0, ((D@@2)(g1))(1) = 0, ((D@@4)(g1))(0) = 0, ((D@@4)(g1))(1) = 0}; res3 := dsolve(`union`(sys2, bcs3), numeric, maxmesh = 2024, abserr = 0.1e-4); plots:-odeplot(res3, [seq([y, (cat(g, i))(y)], i = 1 .. 2)], 0 .. 1)

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

  

``

 

Download error.mw

please help me

thanks...

 
Please Wait...