dr/dt* ln(ardr/dt)=b/r^7 how to solve this equation

Question

vishal vyas 2018-2-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/383783-dr-dt-ln-a-r-dr-dt-b-r-7-how-to-solve-this-equation

评论： Torsten 2018-2-21

kindly help me to solve this equation a = 0.5, b=2, r(0)=1.2

2 个评论
显示无隐藏无

John D'Errico 2018-2-20

I'll suggest you probably won't get much of an answer here, because this is not a question about MATLAB. You are asking how to solve that nonlinear differential equation. So it is a question purely about mathematics, on a problem with no clear solution and I will guess no direct analytical solution. You might catch someone here with an idea, but far more likely to get a result is by asking on a site where the question is on-topic.

Walter Roberson 2018-2-20

在 MATLAB Online 中打开

There is no easy solution for that. The rule is:

r(t) = RootOf(int(P^7*LambertW(1/P^6), P = Z .. 6/5)+2*t)

which is to say that at each point, t, r(t) is the lower bound of the integral P^7*LambertW(1/P^6) such that integrating over P from lower bound to 6/5, plus 2*t, gives 0. (P is an arbitrary variable name here.)

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Roger Stafford 2018-2-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/383783-dr-dt-ln-a-r-dr-dt-b-r-7-how-to-solve-this-equation#answer_306168

在 MATLAB Online 中打开

Here is how I would approach your problem. First we write

a*r*dr/dt*log(a*r*dr/dt) = a*b/r^6

Now define w:

w = log(a*r*dr/dt)

and therefore

a*r*dr/dt = exp(w)

Thus

exp(w)*w = a*b/r^6

Hence

w = lambertw(a*b/r^6)
a*r*dr/dt = exp(lambertw(a*b/r^6))
dr/dt = 1/(a*r)*exp(lambertw(a*b/r^6))

Now finally you have a differential equation in the form that Matlab's ode functions can evaluate numerically, provided you have the lambertw function available.