Solving time dependent differential equation and plotting

pleasehelpmeno
Messages
154
Reaction score
0
Hi i have the differential equation

\frac{d^{2}}{dt^{2}}X(t) +(A+B\frac{sin^{2}(mt)}{mt})X(t) I have tried by hnd to solve this and am getting knowhere does anyone know how to solve it and then plot X against t (where the constants A, B and m will be arbitrarily added), possibly using maple? I know that it is related to the mathieua equations and have found this.

Thanks
 
Last edited by a moderator:
Physics news on Phys.org
Do you need an "analytic" solution? If you just want to graph the solution, I would recommend a "Runge-Kutta" numerical solution.
 
I just want to get a X= f(t) equation that i can then plot, for t values
 
I guess you mean:

\frac{d^{2}}{dt^{2}}X(t) +(A+B\frac{sin^{2}(mt)}{mt})X(t)=0

In which case as always multiply by dX/dt

1/2 \frac{d}{dt}(\frac{dX(t)}{dt})^2 + (A+B\frac{sin^{2}(mt)}{mt})/2 \frac{d}{dt}(X^2(t)) = 0

Now there's a problem, you need to integrate by parts to get the next equation:

1/2 (\frac{dX(t)}{dt})^2 + (A+B\frac{sin^{2}(mt)}{mt})/2 X^2(t) - \int X^2(t)/2 (B\frac{d}{dt}( \frac{sin^{2}(mt)}{mt}) dt = E

So now we get a mixed differentail-integral equation, which I am not sure how to approxiamte.
 
pleasehelpmeno said:
I just want to get a X= f(t) equation that i can then plot, for t values

You can't on the way you are thinking of.
The solutions cannot be expressed as a combination of a finite number of usual functions. In similar cases, some special functions were defined in order to represent the solutions on a closed form (see the example below, in case of a simpler equation). As far as I know, it was not done in the case of your equation. So, do not expect to find an exact analytical equation X=f(t).
Of course, this ODE can be solved and the result can be plot, thanks to numerical computation methods.
Another way for Physicians is to replace the equation by another simpler one which can be analyticaly solved, if the physical phenomena can be modeled with some approximations on some limited ranges of the parameters.
 

Attachments

  • Mathieu functions.JPG
    Mathieu functions.JPG
    12 KB · Views: 515
Thread 'Direction Fields and Isoclines'
I sketched the isoclines for $$ m=-1,0,1,2 $$. Since both $$ \frac{dy}{dx} $$ and $$ D_{y} \frac{dy}{dx} $$ are continuous on the square region R defined by $$ -4\leq x \leq 4, -4 \leq y \leq 4 $$ the existence and uniqueness theorem guarantees that if we pick a point in the interior that lies on an isocline there will be a unique differentiable function (solution) passing through that point. I understand that a solution exists but I unsure how to actually sketch it. For example, consider a...
Back
Top