Rewriting ODE's into lower orders

  • Thread starter Thread starter Euler2718
  • Start date Start date
  • Tags Tags
    matlab ode
Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
2 replies · 1K views
Euler2718
Messages
90
Reaction score
3

Homework Statement



Express

[tex]\frac{d^{2}x}{dt^{2}} + \sin(x) = 0[/tex]

In a system in terms of [itex]x'[/itex] and [itex]y'[/itex].

Homework Equations

The Attempt at a Solution


[/B]
I seen this example:

[tex]x^{\prime\prime\prime} = x^{\prime}(t)\cdot x(t) - 2t(x^{\prime\prime}(t))^{2}[/tex]

Where they then wrote:

[tex]x^{\prime} = y[/tex]
[tex]y^{\prime} = z[/tex]
[tex]z^{\prime} = y^{\prime\prime} = x^{\prime\prime\prime} = x^{\prime}(t)\cdot x(t) - 2t(x^{\prime\prime}(t))^{2}[/tex]
So:
[tex]z^{\prime} = xy - 2tz^{2}[/tex]
Since [itex]y^{\prime\prime} = x^{\prime} = z[/itex]

And thus an Euler's method can be devised in MATLAB, given some IVP's of course. Is this then the correct approach for my problem:

[tex]x^{\prime} = y[/tex]
[tex]y^{\prime} = x^{\prime\prime} = -\sin(x)[/tex]

Not really familiar with ODE's and such processes, but I need to apply this correctly to proceed with my work.
 
Physics news on Phys.org
Morgan Chafe said:

Homework Statement



Express

[tex]\frac{d^{2}x}{dt^{2}} + \sin(x) = 0[/tex]

In a system in terms of [itex]x'[/itex] and [itex]y'[/itex].

Homework Equations

The Attempt at a Solution


[/B]
I seen this example:

[tex]x^{\prime\prime\prime} = x^{\prime}(t)\cdot x(t) - 2t(x^{\prime\prime}(t))^{2}[/tex]

Where they then wrote:

[tex]x^{\prime} = y[/tex]
[tex]y^{\prime} = z[/tex]
[tex]z^{\prime} = y^{\prime\prime} = x^{\prime\prime\prime} = x^{\prime}(t)\cdot x(t) - 2t(x^{\prime\prime}(t))^{2}[/tex]
So:
[tex]z^{\prime} = xy - 2tz^{2}[/tex]
Since [itex]y^{\prime\prime} = x^{\prime} = z[/itex]

And thus an Euler's method can be devised in MATLAB, given some IVP's of course. Is this then the correct approach for my problem:

[tex]x^{\prime} = y[/tex]
[tex]y^{\prime} = x^{\prime\prime} = -\sin(x)[/tex]
This is it right here.
 
  • Like
Likes   Reactions: Euler2718
Mark44 said:
This is it right here.

Thanks. I really had to be sure before continuing.