Advanced Search

Question: not answer for dsolve differential equation

Question:not answer for dsolve differential equation

torabi 85 Maple 2016
dsolve numeric bvp
August 17 2016
1

hi

please help me for dsolve differential equations...after much time dont answer!!!

thanks..

dsolve.mw

Digits := 15; SYS := [.16783*h1(theta)-0.96238e-3*(diff(h1(theta), theta, theta))+0.61603e-1*(diff(h2(theta), theta))+0.14870e-4*(diff(h2(theta), theta, theta, theta))-.23703*h3(theta)-0.84431e-3*(diff(h3(theta), theta, theta))+3.4919*10^(-7)*(diff(h1(theta), theta, theta, theta, theta)) = 0, 2.3940*h2(theta)-.35329*(diff(h2(theta), theta, theta))-0.68260e-1*(diff(h1(theta), theta))-0.16526e-4*(diff(h1(theta), theta, theta, theta))+3.0808*(diff(h3(theta), theta))-0.17833e-2*(diff(h3(theta), theta, theta, theta)) = 0, 9.4813*10^(-7)/((1.+1.5802*10^(-8)*h3(theta))*ln(10.+1.5802*10^(-7)*h3(theta))^2)-3.1867/((1.-0.26556e-1*h3(theta))*ln(10.-.26556*h3(theta))^2)-7.6530/((1.-0.31888e-1*h3(theta))*ln(10.-.31888*h3(theta))^2)-4.2551/((1.-0.35459e-1*h3(theta))*ln(10.-.35459*h3(theta))^2)-9.0315/((1.-0.37632e-1*h3(theta))*ln(10.-.37632*h3(theta))^2)-4.6587/((1.-0.38822e-1*h3(theta))*ln(10.-.38822*h3(theta))^2)-9.4520/((1.-0.39384e-1*h3(theta))*ln(10.-.39384*h3(theta))^2)+0.74999e-1/((1.+0.12500e-2*h3(theta))*ln(10.+0.12500e-1*h3(theta))^2)-.69143/((1.-0.28810e-2*h3(theta))*ln(10.-0.28810e-1*h3(theta))^2)-.12945/(1.-0.38836e-1*h3(theta))^4-0.38260e-1/(1.-0.11478e-1*h3(theta))^4-0.24826e-1/(1.-0.37240e-2*h3(theta))^4+0.74712e-3*(diff(h3(theta), theta, theta, theta, theta))+2.6337*10^(-8)/(1.+1.5802*10^(-8)*h3(theta))^4-0.27242e-1*(diff(h3(theta), theta, theta))-3.0707*(diff(h2(theta), theta))+0.17833e-2*(diff(h2(theta), theta, theta, theta))-.23618*h1(theta)-0.84126e-3*(diff(h1(theta), theta, theta))-0.89222e-1/(1.-0.26767e-1*h3(theta))^4-.21340/(1.-0.32010e-1*h3(theta))^4-.13154/(1.-0.19732e-1*h3(theta))^4-0.88519e-1/(1.-0.26556e-1*h3(theta))^4-4.7454/((1.-0.39545e-1*h3(theta))*ln(10.-.39545*h3(theta))^2)-0.36390e-1/(1.-0.10917e-1*h3(theta))^4-1.3100/((1.-0.10917e-1*h3(theta))*ln(10.-.10917*h3(theta))^2)-.89374/((1.-0.37240e-2*h3(theta))*ln(10.-0.37240e-1*h3(theta))^2)-1.3773/((1.-0.11478e-1*h3(theta))*ln(10.-.11478*h3(theta))^2)-.11820/(1.-0.35459e-1*h3(theta))^4-.26259/(1.-0.39388e-1*h3(theta))^4-4.7356/((1.-0.19732e-1*h3(theta))*ln(10.-.19732*h3(theta))^2)+28.586*h3(theta)-.21258/(1.-0.31888e-1*h3(theta))^4-9.4531/((1.-0.39388e-1*h3(theta))*ln(10.-.39388*h3(theta))^2)-4.6603/((1.-0.38836e-1*h3(theta))*ln(10.-.38836*h3(theta))^2)-9.0393/((1.-0.37664e-1*h3(theta))*ln(10.-.37664*h3(theta))^2)-4.2631/((1.-0.35526e-1*h3(theta))*ln(10.-.35526*h3(theta))^2)-7.6824/((1.-0.32010e-1*h3(theta))*ln(10.-.32010*h3(theta))^2)-.25088/(1.-0.37632e-1*h3(theta))^4-.12919/(1.-0.19378e-1*h3(theta))^4-.11842/(1.-0.35526e-1*h3(theta))^4-.12941/(1.-0.38822e-1*h3(theta))^4+0.20833e-2/(1.+0.12500e-2*h3(theta))^4-0.19206e-1/(1.-0.28810e-2*h3(theta))^4-4.6507/((1.-0.19378e-1*h3(theta))*ln(10.-.19378*h3(theta))^2)-.26255/(1.-0.39384e-1*h3(theta))^4-.13182/(1.-0.39545e-1*h3(theta))^4-.25109/(1.-0.37664e-1*h3(theta))^4-3.2120/((1.-0.26767e-1*h3(theta))*ln(10.-.26767*h3(theta))^2) = 0]

ode2 := diff(SYS[2], theta); SYS2 := {ode2, SYS[1], SYS[3]}; bcs2 := {h1(0) = 0, h1(1) = 0, h2(0) = 0, h2(1) = 0, h3(0) = 0, h3(1) = 0, ((D@@1)(h1))(0) = 0, ((D@@1)(h1))(1) = 0, ((D@@1)(h3))(0) = 0, ((D@@1)(h3))(1) = 0}; bcs22 := eval[recurse](convert(SYS[2], D), `union`({theta = 1}, bcs2)); res2 := dsolve(`union`(`union`(SYS2, bcs2), {bcs22}), 'maxmesh' = 2024, numeric, method = bvp[middefer], range = 0 .. 1, abserr = 0.1e-3, output = listprocedure)

NULL



Download dsolve.mw

 
Please Wait...