Question: I have tried to solve the following ode equation, but I have got error. How can i fix it?


 

restart

mu:=0.0025

0.25e-2

(1)

A:=1*10^(-20); Aprim:=4.46*10^(-20);

1/100000000000000000000

 

0.4460000000e-19

(2)

eta:=1.003*10^(-3)

0.1003000000e-2

(3)

V:=2*10^(-9)

1/500000000

(4)

ode:=mu*h(x)*diff(h(x),x,x)+mu*A/(6*pi*h(x)^2)-eta*V

0.25e-2*h(x)*(diff(diff(h(x), x), x))+0.4166666667e-23/(pi*h(x)^2)-0.2006000000e-11

(5)

ode1:=mu*h(x)^3*diff(h(x),x,x)-eta*V*h(x)^2+mu*A/(6*pi)=0

0.25e-2*h(x)^3*(diff(diff(h(x), x), x))-0.2006000000e-11*h(x)^2+0.4166666667e-23/pi = 0

(6)

ic1:=h(0)=3/1000000000, D(h)(-0.20e-4)=0,h(-0.20e-4)=1e-9

h(0) = 3/1000000000, (D(h))(-0.20e-4) = 0, h(-0.20e-4) = 0.1e-8

(7)

dsol1 := dsolve({ic1, ode}, numeric)

Error, (in dsolve/numeric/bvp) system is singular at left endpoint, use midpoint method instead

 

``


 

Download velocity_tehta1.mwvelocity_tehta1.mw

Please Wait...