Third order linear ODE Wolfram Math. problem

  • Thread starter Thread starter dperkovic
  • Start date Start date
  • Tags Tags
    Linear Ode
AI Thread Summary
The discussion centers on solving the homogeneous ordinary differential equation (ODE) d^3x/dt^3 + 2dx/dt = 0 using the substitution x = e^(λt). The characteristic polynomial derived is λ^3 + 2λ = 0, yielding one real root (λ1 = 0) and a pair of complex roots (λ2, λ3 = ±√2i). The general solution is expressed as x(t) = C1 + C2cos(√2t) + C3sin(√2t). A discrepancy arises when comparing this solution to the output from Wolfram Mathematica, which includes different constants, attributed to Mathematica's method of finding general solutions that may involve more complex forms. Ultimately, the differences in constants do not affect the validity of the solutions, as integration constants are arbitrary.
dperkovic
Messages
17
Reaction score
0
I have a problem with this homogenous ODE:
<br /> \frac{d^3x}{dt^3} + 2\frac{dx}{dt} = 0$<br />

With substitution x=e^\lambda t, I got characteristic polinom: \lambda^3+2\lambda=0.
So,
\lambda_1 = 0, \lambda_{1,2} = \pm\sqrt{2}i

Because we have one real,and one pair of complex roots of characteristic polinom, general solution of ODE is:

x(t) =C1e^{0\cdot t} + C2\cos(\sqrt{2}t) +C3\sin(\sqrt{2}t),
or:
x(t) =C1 + C2\cos(\sqrt{2}t) +C3\sin(\sqrt{2}t)

But, Wolfram Mathematica's command DSolve[x'''[t] + 2*x'[t] == 0, x[t], t],give to me this:

x(t) =C1 - \frac{C2\cos(\sqrt{2}t)}{\sqrt{2}} +\frac{C3\sin(\sqrt{2}t)}{\sqrt{2}}.

As we can see, there is difference between my and wolfram's solution. Why ?? Where is a mistake in my procedure ?
 
Last edited:
Physics news on Phys.org
There is no difference in the solutions, remember the constants of integration are arbitrary.
 
Thank's. I thought that. But still, I can't figure out why constants in Mathematica's solution is different than constants i my solution. From which algorithm comming this "complicated" constants ?
 
The Mathematica DSolve usually looks for general solutions that match your particular equation and then it uses the form of the general solution to solve your specific problem. Often the general solution has such "complicated" constants.
 
In this case the general form is probably something like:
a x&#039;&#039;&#039;(t)+b x&#039;(t)=0

Which has the general solution:
\frac{\text{C1} \sqrt{a} \sin \left(\frac{\sqrt{b}<br /> t}{\sqrt{a}}\right)}{\sqrt{b}}-\frac{\text{C2} \sqrt{a} \cos<br /> \left(\frac{\sqrt{b} t}{\sqrt{a}}\right)}{\sqrt{b}}+\text{C3}

Which turns into exactly the expression you have above for a=1 and b=2
 
This is the explanation which I looking for :smile: :smile: :smile:. Thank you !
 
Thread 'Question about pressure of a liquid'
I am looking at pressure in liquids and I am testing my idea. The vertical tube is 100m, the contraption is filled with water. The vertical tube is very thin(maybe 1mm^2 cross section). The area of the base is ~100m^2. Will he top half be launched in the air if suddenly it cracked?- assuming its light enough. I want to test my idea that if I had a thin long ruber tube that I lifted up, then the pressure at "red lines" will be high and that the $force = pressure * area$ would be massive...
I feel it should be solvable we just need to find a perfect pattern, and there will be a general pattern since the forces acting are based on a single function, so..... you can't actually say it is unsolvable right? Cause imaging 3 bodies actually existed somwhere in this universe then nature isn't gonna wait till we predict it! And yea I have checked in many places that tiny changes cause large changes so it becomes chaos........ but still I just can't accept that it is impossible to solve...
Back
Top