Crazy ODE, can it be solved?

  • Thread starter Brin
  • Start date
  • #1
73
0
Can anyone tell me if this ODE has an analytic solution? And if it does, how the heck might I go about it?

[tex]
\left(\frac{1}{y^{2}}\frac{dy}{dx}\right)^{2}-\frac{A}{y^{3}}-\frac{B}{y^{2}}=D
[/tex]
 
Last edited:

Answers and Replies

  • #2
Ben Niehoff
Science Advisor
Gold Member
1,879
162
It looks like it is separable. Just isolate dy/dx. It looks like the solution will be some kind of elliptic function.

(Mathematica can find an analytic solution for x as a function of y. It involves elliptic functions and finding roots of a cubic polynomial.)
 
  • #3
798
34
The integral of :
dx = Sqrt[Ay+By²+Dy^4] dy
involves elliptic integrals.
In the general case, the function x(y) which can be obtained would be rather complicated. Then, inverting it in order to express y(x) would be a big chore. Better use numerical integration.
 
  • #4
1,796
53
. . . ohhhhhh . . . you guys give up too easy. Assume we are given that the solution to:

[tex]\left(\frac{dz}{dx}\right)^2=(1-z^2)(1-k^2z^2)=\Delta^2(z)[/tex]

is:

[tex]y=sn(z,k)[/tex]

where sn is the Jacobi elliptic function. Then we seek a transformation z=z(y) that transforms:

[tex]
\begin{align*}
\left(\frac{dy}{dx}\right)^2&=a+by+cy^2+dy^3+ey^4\\
&=h^2(y-\alpha)(y-\beta)(y-\gamma)(y-\delta) \\
&=h^2 \Delta_2^2(y)
\end{align*}
[/tex]

into this standard form. To this end we let:

[tex]z^2=\frac{(\beta-\gamma)}{(\alpha-\delta)}\frac{(y-\alpha)}{(y-\beta)}=s\frac{(y-\alpha)}{(y-\beta)},\quad k^2=\frac{(\beta-\delta)}{(\alpha-\gamma)}\frac{(\alpha-\delta)}{(\beta-\delta)},\quad M^2=\frac{(\beta-\delta)(\alpha-\delta)}{4}[/tex]

for which we obtain:

[tex]\frac{1}{\Delta(z)}\frac{dz}{dx}=\frac{M}{\Delta_2(y)}\frac{dy}{dx}=Mh[/tex]

so that:

[tex]\frac{dz}{dx}=Mh\Delta(z)=Mh\sqrt{(1-z^2)(1-k^2z^2)}[/tex]

and therefore:

[tex]z=sn(hMv,k),\quad v=x-x_0[/tex]

or:

[tex]y=\frac{z^2\beta-s\alpha}{z^2-s}[/tex]

I believe though the actual implementation of this would be difficult as I have never worked a real problem using this method but I think would be a nice project for someone taking non-linear DEs next semester. :)
 
  • #5
798
34
Hi !
Very nice job jackmell, but . . .
. . . ohhhhhh . . . you guys give up too easy.
. . . ohhhhhh . . . even easier . . . just a few seconds to have the explicit formula :
I let you try WolframAlpha and see the result :rofl: :surprised
 

Attachments

  • WolframAlpha.JPG
    WolframAlpha.JPG
    10.4 KB · Views: 360

Related Threads on Crazy ODE, can it be solved?

  • Last Post
Replies
4
Views
2K
  • Last Post
Replies
2
Views
2K
  • Last Post
Replies
2
Views
2K
Replies
2
Views
2K
  • Last Post
Replies
4
Views
1K
Replies
4
Views
2K
  • Last Post
Replies
3
Views
2K
Replies
1
Views
648
Replies
2
Views
2K
  • Last Post
Replies
1
Views
1K
Top