melika mokhtari

10 Reputation

One Badge

0 years, 194 days

MaplePrimes Activity


These are questions asked by melika mokhtari

Hi,

i am trying to solve a differential equation numercially but since it is second degree,I can not achieve the proper answer.What should I do?

thank you
diff1.mw
 

restart

h := 1-.8*x+.5*(x^2-x)

1-1.3*x+.5*x^2

(1)

z := 3^(3/2)*(((2+k*De^2*(diff(p(x), x))^2*h^2)^3*h^2+108*k*De^2)/k)^(1/2)

3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2)

(2)

c := (1/6)*((-54*De+z)*k^2*h^2)^(1/3)/(k*De*h)-(1/2)*h*(2+k*De^2*(diff(p(x), x))^2*h^2)/(De*((-54*De+z)*k^2*h^2)^(1/3))-(1/2)*h*(diff(p(x), x))

(1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x))

(3)

ode1 := ((1/2)*c*(diff(p(x), x))^2*h^4+(diff(p(x), x))*c^2*h^3+(1/10)*(diff(p(x), x))*h^5+c^3*h^2)*k*De^2+(1/2)*c*h^2+(1/6)*(diff(p(x), x))*h^3+h = 0

((1/2)*((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^4+(diff(p(x), x))*((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))^2*(1-1.3*x+.5*x^2)^3+(1/10)*(diff(p(x), x))*(1-1.3*x+.5*x^2)^5+((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))^3*(1-1.3*x+.5*x^2)^2)*k*De^2+(1/2)*((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))*(1-1.3*x+.5*x^2)^2+(1/6)*(diff(p(x), x))*(1-1.3*x+.5*x^2)^3+1-1.3*x+.5*x^2 = 0

(4)

ivp := {ode1, p(0) = 0}

{((1/2)*((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^4+(diff(p(x), x))*((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))^2*(1-1.3*x+.5*x^2)^3+(1/10)*(diff(p(x), x))*(1-1.3*x+.5*x^2)^5+((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))^3*(1-1.3*x+.5*x^2)^2)*k*De^2+(1/2)*((1/6)*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3)/(k*De*(1-1.3*x+.5*x^2))-(1/2)*(1-1.3*x+.5*x^2)*(2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)/(De*((-54*De+3*3^(1/2)*(((2+k*De^2*(diff(p(x), x))^2*(1-1.3*x+.5*x^2)^2)^3*(1-1.3*x+.5*x^2)^2+108*k*De^2)/k)^(1/2))*k^2*(1-1.3*x+.5*x^2)^2)^(1/3))-(1/2)*(1-1.3*x+.5*x^2)*(diff(p(x), x)))*(1-1.3*x+.5*x^2)^2+(1/6)*(diff(p(x), x))*(1-1.3*x+.5*x^2)^3+1-1.3*x+.5*x^2 = 0, p(0) = 0}

(5)

dsolve(ivp, p(x), numeric, parameters = [k, De])

Error, (in dsolve/numeric/make_proc) Could not convert to an explicit first order system due to 'RootOf'

 

``


 

Download diff1.mw

 

 


 

Download diff1.mw

Page 1 of 1