![Klada10000.png](https://cdn.steemitimages.com/DQmZo1uTZnk1n9bRfDpyWPuBqoX5RwSpLCmMRqxMJzHqFQZ/Klada10000.png)

Hi guys. Today I am back with a whole new topic I have never talked about. I intend to do a lot of posts like this as I think it is very educational and I personally love solving problems in classical mechanics. A lot can be done analytically which is a true joy!

So let me answer the question I posted earlier. There is a simpler more systematic way of describing systems of classical objects and it is called [The Lagrange Formalism](https://en.wikipedia.org/wiki/Lagrangian_mechanics). I wont talk about what it is as it is a very well described topic that you can find on the internet, but I do intend to showcase its strength in this post and some more posts coming in the future.

The central object of this formalism is the Lagrange function that is a generalized function of the coordinates of the objects in the system - q<sub>i</sub> - and its first ***time*** derivatives -  q<sub>i</sub><sup>,</sup>. It is defined as:

 <center> L(q<sub>i</sub>, q<sub>i</sub><sup>,</sup>) = T(q<sub>i</sub>, q<sub>i</sub><sup>,</sup>) - V(q<sub>i</sub>, q<sub>i</sub><sup>,</sup>)</center>

The **T(q<sub>i</sub>, q<sub>i</sub><sup>,</sup>)** is the kinetic energy of the entire system and **V(q<sub>i</sub>, q<sub>i</sub><sup>,</sup>)** is the potential energy of the system. Once you defined your **L** for your system, you can very easily get your complete set of equations of motion using:

<center>d/dt(&delta;L/&delta; q<sub>i</sub><sup>,</sup>) = &delta;L/&delta; q<sub>i</sub></center>

I am sorry about the weird notation but Markdown is not very math friendy. &delta; stands for a partial derivative, while d/dt is the total time derivative.

Lets look at a simple example of a mass tethered to a wall with a spring with its coefficient being ***k*** and its mass ***m***. You can look at the sketch in the first picture above. The spring has a natural length of ***l***. The kinetic and potential energies of such a system is:

<center> T = 1/2 m (x<sup>,</sup>)<sup>2</sup></center>

<center> V = 1/2 k (x - l)<sup>2</sup></center>

<center> L  = 1/2 m (x<sup>,</sup>)<sup>2</sup> -1/2 k (x - l)<sup>2</sup> </center>

Here we remembered that the velocity is defined as the first time derivative of the position  - v = x<sup>,</sup>. Usually dots are used as time derivatives but it is impossible to use them here. For the next post I will think of something else. If you now use the recipe of how one derives the equations of motion you get:

<center> m x<sup>,,</sup> = -k (x - l)</center>

You get the well known second law of motion derived by Newton. The Lagrangian formalism works. Yaaaay. Now we can use it for a more interesting case.

![Klada2nacin0000.png](https://cdn.steemitimages.com/DQme9WExisHnyy8MsiVUXowhVXbgRvh8ohesiSaPM6kCWbF/Klada2nacin0000.png)

Lets look at a more complex system that contains two bodies with masses that are coupled together. Similarly as before we can easily write the Lagrange function for the system as:

<center> T = 1/2 m<sub>1</sub> (x<sub>1</sub><sup>,</sup>)<sup>2</sup> + 1/2 m<sub>2</sub> (x<sub>2</sub><sup>,</sup>)<sup>2</sup></center>

<center> V = 1/2 k<sub>1</sub> (x<sub>1</sub>)<sup>2</sup> + 1/2 k<sub>2</sub> (x<sub>2</sub>)<sup>2</sup> + 1/2 K (x<sub>1</sub> - x<sub>2</sub>)<sup>2</sup></center>

<center> L  = 1/2 m<sub>1</sub> (x<sub>1</sub><sup>,</sup>)<sup>2</sup> + 1/2 m<sub>2</sub> (x<sub>2</sub><sup>,</sup>)<sup>2</sup> - (1/2 k<sub>1</sub> (x<sub>1</sub>)<sup>2</sup> + 1/2 k<sub>2</sub> (x<sub>2</sub>)<sup>2</sup> + 1/2 K (x<sub>1</sub> - x<sub>2</sub>)<sup>2</sup>) </center>

The kinetic energy part is pretty self-explanatory. The potential energy is also not overly complicated but still a bit more thinking goes into it. The first two terms are the potential spring energies of the wall-mass couplings while the last term is the mass-mass coupling. If x<sub>1</sub> = x<sub>2</sub> then the middle spring didnt stretch. This is why there is a minus in the last term. If we again use the recipe we get the equations of motion:

<center>x<sub>1</sub><sup>,,</sup> = (-k<sub>1</sub>/m<sub>1</sub> - K/m<sub>1</sub>) x<sub>1</sub><sup>,,</sup> + K/m<sub>1</sub> x<sub>2</sub><sup>,,</sup></center>

<center>x<sub>2</sub><sup>,,</sup> = K/m<sub>2</sub> x<sub>1</sub><sup>,,</sup> + (-k<sub>2</sub>/m<sub>2</sub> - K/m<sub>2</sub>) x<sub>2</sub><sup>,,</sup></center>

Although there is a formalism that easily deals with such equations I will solve these equations using some magic. But first lets make a simplification just for calculation (and animation) sake. Lets say that:

<center>  k<sub>1</sub> =  k<sub>2</sub> = k</center>
<center> K = 2 k </center>
<center>  m<sub>1</sub> =  m<sub>2</sub> =m</center>

This means we can define a characteristic frequency &omega;<sub>0</sub><sup>2</sup> = k/m and we can write it to write a set of beautiful equations:

<center>x<sub>1</sub><sup>,,</sup> = -3 &omega;<sub>0</sub><sup>2</sup> x<sub>1</sub>  + 2 &omega;<sub>0</sub><sup>2</sup>x<sub>2</sub></center>

<center>x<sub>2</sub><sup>,,</sup> = 2 &omega;<sub>0</sub><sup>2</sup> x<sub>1</sub>  - 3 &omega;<sub>0</sub><sup>2</sup>x<sub>2</sub></center>

So here comes the magic so pay attention. We notice that there is a symmetry to the equations, which is why we define two new coordinates (x<sub>1</sub> - x<sub>2</sub>) and (x<sub>1</sub> + x<sub>2</sub>). We can write the equations of motion for these new coordinates by subtracting the old differential equations for the first of the new coordinates and by summing up the old equations for the second case. If we do this we get:

<center>(x<sub>1</sub> + x<sub>2</sub>)<sup>,,</sup> = - &omega;<sub>0</sub><sup>2</sup>(x<sub>1</sub> + x<sub>2</sub>)</center>

<center>(x<sub>1</sub> - x<sub>2</sub>)<sup>,,</sup> = - 5 &omega;<sub>0</sub><sup>2</sup>(x<sub>1</sub> - x<sub>2</sub>)</center>

We recognize these equations of motion as equations of motion describing pendulums and similar ***harmonic oscillator*** things. We see that the sum has will oscillate with a frequency of &omega;<sub>0</sub> while the other will oscillate with a frequency of about 2.23 &omega;<sub>0</sub>. But what are these coordinates actually representing.

Well let my gifs illustrate.

<div class="pull-right">

<center>![1nacinKlada.gif](https://cdn.steemitimages.com/DQmQALo4cK41TdRVVkyxm2zMFcYSXLgY6fv1iUi8jCGiKzQ/1nacinKlada.gif)
</center>

</div>
<div class="pull-left">

<center>![2nacinklada.gif](https://cdn.steemitimages.com/DQmbkuF3ZDgfRNow8qWDW6TUM9jdHnVLUr5KLNHcwGThnma/2nacinklada.gif)
</center>

</div>

This is what we called the eigen-motions of a harmonic oscillator, which means that they are stable. If nothing disturbs them they will go on forever. All the other possible motions of this two-body coupled harmonic oscillator can be written as a linear sum of these two eigen-motions. The eigen-frequency  &omega;<sub>0</sub> is the frequency of the system when only the center of mass is moving (right figure) while the eigen-frequency  2.23 &omega;<sub>0</sub>  is the frequency of the system when the masses are moving towards each other (left figure).

# <center> Conclusion </center>

I hope I was able to convince you that this is an easier way of deriving the equations of motion than using the classic Newtonian kinematics. Of course the results are the same and you still have to solve some differential equations, but the way of constructing the differential equations is more systematic. I suggest you do the calculations alone the old way and the Lagrangian formalism way I introduced here for the coupled oscillator problem. You will see it is much easier this way.

I intend to bring you more interesting systems of coupled "stuff". I think it is really fun solving such easy problems that are really rich with physics.

<center>![stemq-concept2.png](https://cdn.steemitimages.com/DQmZLLVZyQ23GWzs3FNU6hzm1Ai5frMnbuRDYQ8C7ZBABGr/stemq-concept2.png)</center>

![Steem1.gif](https://steemitimages.com/DQmf18V3PQenbUzYfkZcWSVGgmV5CtZNymi824SYwAQ2GTp/Steem1.gif)

**StemQ Notice:** *This post was originally submitted on [StemQ.io](http://www.stemq.io), a Q&A application for STEM subjects powered by the Steem blockchain.*