Non linear ODE: y'' = a y + b y^3


by galuoises
Tags: linear
galuoises
galuoises is offline
#1
Jan12-13, 07:54 AM
P: 8
I would like to solve the non linear ODE
[itex] \frac{d}{dx}f(x)=a f(x)+ b f^3 (x)[/itex]
with the boundary
[itex]f(0)=0\quad f(+\infty)=f_0 [/itex]

How to find analitical solution?
Phys.Org News Partner Science news on Phys.org
Going nuts? Turkey looks to pistachios to heat new eco-city
Space-tested fluid flow concept advances infectious disease diagnoses
SpaceX launches supplies to space station (Update)
JJacquelin
JJacquelin is offline
#2
Jan12-13, 11:42 AM
P: 745
If it is f ' = a f+b f^3 then integrate dx = df/(a f+b f^3)
If it is y'' = a y+b y^3 then the first step is the integration of :
y'' y' = (a y+b y^3) y'
y' = a y + b (y^4)/2 +c
dx = dy/sqrt(a y + b (y^4)/2 +c)
The integration will lead to x(y) on the form of elliptic functions. Very complicated to invert in order to obtain y(x)
galuoises
galuoises is offline
#3
Jan13-13, 06:02 AM
P: 8
Pardon me, I write the uncorrect differential equation: the problem is at the second order

[itex] \frac{d^2}{dx^2}f(x)=a f(x) + b f^3 (x) [/itex]
with the boundary
[itex] f(0)=0,\ f(+\infty)=f_0[/itex]

jackmell
jackmell is offline
#4
Jan13-13, 08:12 AM
P: 1,666

Non linear ODE: y'' = a y + b y^3


In my opinion, [itex]x(y)=g(y,c_1,c_2)[/itex] is an analytical expression for the solution but I think first, just scrap the a and b and look at:

[tex]y''=y+y^3[/tex]

then do what Jacq said and get the expression in terms of:

[tex]x(y)=g(y,c_1,c_2)[/tex]

then try and solve simultaneously the expressions:

[tex]0=g(0,c_1,c_2)[/tex]

[tex]g(f_0,c_1,c_2)\to+\infty[/tex]

for the constants [itex]c_1,c_2[/itex] and if need be, do so numerically for them just for starters.
JJacquelin
JJacquelin is offline
#5
Jan13-13, 11:52 AM
P: 745
Hello !

It is a difficult problem on the analytical viewpoint.
A difficulty is that elliptic functions are involved. But the major difficulty is due to the condition f(x=+infinity)=finite constant.
I am not quite sure that the solution given in attachment is a correct answer to the problem, so it should be carefully verified.
Attached Thumbnails
ODE.JPG  
JJacquelin
JJacquelin is offline
#6
Jan14-13, 03:37 AM
P: 745
In case of some complementary conditions specified in attachment, y(x) can be fully explicited thanks to a Jacobi elliptic function.
Attached Thumbnails
Result.JPG  
jackmell
jackmell is offline
#7
Jan14-13, 07:46 AM
P: 1,666
Nicely done Jacquelin and beautiful too. Thanks.
galuoises
galuoises is offline
#8
Feb9-13, 06:28 AM
P: 8
Thank you so much!


Register to reply

Related Discussions
Linear Algebra question regarding linear operators and matrix rep. relative to basis Calculus & Beyond Homework 6
To Interpret Solving Systems of Linear Equations Geometrically in Terms of Linear Alg Linear & Abstract Algebra 0
Linear Algebra: Linear Independence and writing Matrices as linear combinations Calculus & Beyond Homework 7
Linear Simple Harmonic Oscillator: period a direct linear proportion to mass? General Physics 1
Linear Algebra: Linear Transformation and Linear Independence Calculus & Beyond Homework 8