
#1
Sep209, 05:57 PM

P: 54

Hello science companions ;)
I was wondering for a while how to solve this ODE (y''=f(y)) and I finally found one way to integrate it. I don't know actually if it is correct but as I see them, all steps are mathematically consistent. Here I show how I came to the general solution (for all steps consider y=y(x): 1st multiply each side of the ODE by 2y' 2·y'·y''=2·f(y)·y' 2nd using (y'^{2})'=2·y'·y'' write (y'^{2})'=2·f(y)·y' 3rd first bothsidewithrespecttox integration [tex]\int[/tex](y'^{2})'dx=2[tex]\int[/tex]f(y)y'dx+A (1st integration costant) 4th cancelling derivative with integral and using dy=y'dx y'^{2}=2[tex]\int[/tex]f(y)dy+A 5th the integral is the primitive of f(y) so we call it F(y), and aplying squareroots to both sides we finally have: y'=[tex]\pm[/tex][tex]\sqrt{2F(y)+A}[/tex] These are two ecuations that only differ in sign and represent two "families of families" of solutions. In any of them the equation is separable as we have the first derivative equal to an only function of y, so reorganizing we may write: [tex]\frac{1}{\sqrt{2F(y)+A}}·y'=\pm1[/tex] Then integration as usual which gives, using again dy=y'dx: [tex]\int\frac{1}{\sqrt{2F(y)+A}}dy=\pm x+B[/tex] (2nd integration constant) What represents the general set of solutions if I wasn't mistaken. Please any comments, corrections, advise... will be more than welcome. Hope I made myself clear enough ;) 



#2
Sep209, 07:02 PM

Sci Advisor
HW Helper
PF Gold
P: 12,016

This looks correct.




#3
Sep209, 07:13 PM

P: 54

Yeh, that's what I thought, I just wasn´t sure. Thanks a lot for reading.




#4
Sep209, 07:21 PM

P: 813

Opinion about y"=f(y)
I think there is a mistake from step three to step four on the left hand side of the equation. I don't think things cancel how you think they do.




#5
Sep209, 07:24 PM

P: 54

Why not? I mean, is [tex]\int\frac{d(y'^2)}{dx}dx=y'^2 [/tex] not true? Is that where you say i've mistaken?




#6
Sep209, 07:28 PM

P: 813





#7
Sep209, 07:32 PM

P: 54

ok, ok nothing :) John , good night




#8
Sep209, 07:35 PM

P: 813

So, it looks right. Why not try try applying it to some differential equations for which we know the solution. I'd be interested to see a few examples worked out.




#9
Sep309, 03:44 AM

P: 54

Yeh! (why didn't I thought about that before!, thanks)For example we know that the solution of ODE
[tex]\frac{d^2y}{dx^2}=y , y(0)=1 , \frac{dy}{dx}(0)=0[/tex] is y=coshx (integrating it with the usual method for linear2ndorederhomogeneus ODE). The method I exposed will give the solution by evaluating the integral (with f(y)=y and, therefore, F(y)=y^{2}/2): [tex]\int\frac{1}{\sqrt{A+2·\frac{y^2}{2}}}dy=\int\frac{1}{\sqrt{A+y^2}}dy[/tex] which can be written as the standard integral for acoshy (by doing A=a^{2}). And hence, all considering inital conditions for y(0) and y'(0), we'll have that y=coshx C.Q.D. So I thinks this really works. In fact, by the way, once I came to it I couldn't help to think of the pendulum equation, not the one for little variations of [tex]\theta[/tex] but the "generic" one, that is to say a pendulum of lenght l with gravity g that satisfies: [tex]\frac{d^2\theta}{dt^2}=\frac{g}{l}sin\theta[/tex] Unfortunately you finally have to deal with the following integral [tex]\int\frac{1}{\sqrt{Csinu}}du[/tex] which, I think, have no closed solution. Salutations, zaphys ;) 



#10
Sep309, 06:30 AM

Math
Emeritus
Sci Advisor
Thanks
PF Gold
P: 38,900





#11
Sep309, 07:20 AM

P: 54

oh! so pleasant to hear that :D, thanks Hallsofivy. I was just preparing ODEs before I get to university so that I can deal there with sightly complicated problems, and wondered if was possible to integrate y''=f(y), and it is, wonderfull :)




#12
Sep309, 08:23 AM

P: 608

Another way to think of this ODE y'' = f(y) ...
Since x does not appear explicitly in the equation, write z = y' and write it as an ODE for z as a function of y. Then: [tex]y'' = \frac{dz}{dx} = \frac{dz/dy}{dx/dy} = \frac{dz}{dy}\;\frac{dy}{dx} = \frac{dz}{dy}\;z[/tex] So we have to solve the firstorder ODE [tex]\frac{dz}{dy}\;z = f(y)[/tex] which is separable. When you solve that as z = G(y), then convert back to the original variables y' = G(y) and the final result is an integration. 



#13
Sep309, 08:01 PM

P: 54

Yes absolutely. And pretty more elegant than my method ;) I can see that the final result has to deal with the same integral as my final result.
I'm so glad to have done this post :D 


Register to reply 
Related Discussions  
Difference between "Identical", "Equal", "Equivalent"  Calculus & Beyond Homework  9  
Anyone familiar with "centrifugal potential" and "brachistochrone" in polar coords?  Advanced Physics Homework  7 