Ah, ****, here we go again.
The fundamental theorem of engineering.
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=http%3A%2F%2Fpm1.aminoapps.vertvonline.info%2F7617%2F9e548d5623ce35969a9c6098fdcd56d6ee847b4br1-1672-945v2_hq.jpg)
Please don't use it, guys. I mean, no, but...
------------------------
sin x = x, really?
Well, before I start ranting again...
Let's take a look at what the second result of the fundamental theorem of engineering actually points at.
sin x = x
Let f: R -> R be a function. Let's suppose that f is smooth enough, so it's continuous and differentiable as many times as you want on R.
So, if we take a point t, and h small enough so that t+h is really close to x, it's possible to approximate f(t+h) in function of f(t).
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=http%3A%2F%2Fpm1.aminoapps.vertvonline.info%2F7617%2Fa272bee28c406ac4b4008fd31264a2da54298c17r1-1189-823v2_hq.jpg)
As you can see, there's still some error in that approximation. But in order to lessen the error, you need to understand this..
We can try to approximate a function around a point by using some polynomial expression., which coefficients include f and its successive derivatives. It works the same way as this.
Of course, the more you haven, the less error you get.
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=http%3A%2F%2Fpm1.aminoapps.vertvonline.info%2F7617%2F6b3f920f40f25feb2958fee2ca4fe37dfc48deb4r1-1019-126v2_hq.jpg)
This is called Taylor expanding a function at a point t0.
Now, how does that apply to sin t = t?
Well, by expanding with 1 extra term near t0 = 0, as h -->0:
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=http%3A%2F%2Fpm1.aminoapps.vertvonline.info%2F7617%2Fd95b2ef4fe5b8375d0e5e4e9d07e69a9f3309d23r1-435-217v2_hq.jpg)
There, a rather minimalist approach to the fundamental theorem of engineering.
---------------------------------------------------------
Application: Euler's method
Taylor expanding can come very handy if you want to approximate solutions within the bounds in which the Taylor expansion remains accurate enough.
Let's see an application of that: falling down with air friction.
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=https%3A%2F%2Fpa1.aminoapps.vertvonline.info%2F7617%2F8efa71a6fd58e5b0e10879b8ec6c0783b4d2d3d0r1-359-202_hq.gif)
After having been summoned to Disboard, Sora and Shiro (modeled by a point M at height z(t)) fall down from the sky without initial speed.
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=http%3A%2F%2Fpm1.aminoapps.vertvonline.info%2F7617%2Fe4d9d66e12b26533214b1bbba9a125aea0603fdar1-736-550v2_hq.jpg)
System: M
Frame of referential: Terrestrial, supposed to be galilean
Forces: The weight P and the friction force F, as shown above.
Note: f is a constant here. F is proportional to the square of the speed for high speed falls. For slower falls, it's better to model F as proportional to the speed.
Furthermore, we should expect the speed to converge because the fall gets slown down by the air friction.
Newton's second law, projected on ez, gives, for any t:
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=http%3A%2F%2Fpm1.aminoapps.vertvonline.info%2F7617%2Fad2ca278ad491c1bdbcfeb759b7b6b4eb5f7f1b7r1-340-222v2_hq.jpg)
We don't know how to solve non linear differential equations.
But we can at least approximate the solutions more or less accurately.
-------------------------------------------
Now, assuming that during the fall v can be modeled as a smooth function, Taylor expanding to the first order for a small turation h would give:
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=http%3A%2F%2Fpm1.aminoapps.vertvonline.info%2F7617%2Fdbc1ea03302381e473f2eed0b8761dd69f246211r1-416-199v2_hq.jpg)
So, what does this mean here?
Suppose that you know what at time t, what v(t) and dv/dt(t) are.
Then, you can approximate v(t+h) using the formula above, then approximate dv/dt(t+h) by using the differential equation:
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=http%3A%2F%2Fpm1.aminoapps.vertvonline.info%2F7617%2F5b4c1d04c8655362e15b8e6109e23a31ea73309dr1-358-108v2_hq.jpg)
Now, doing this step by step, you can actually more or less find the numeric values of v and v' at the times t = nh, for n an integer. You just have to guess the value of a point from the value of the previous one.
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=https%3A%2F%2Fpa1.aminoapps.vertvonline.info%2F7617%2Fb65b15f669fa1ec363beb467d6360622b487a8e2r1-478-437_hq.gif)
Now, let's punch that onto Python. To ease the problem, we'll name K = f/m.
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=http%3A%2F%2Fpm1.aminoapps.vertvonline.info%2F7617%2F86d035a68608710052ee51c5be610b46c11052a9r1-822-611v2_hq.jpg)
Now. Punching that onto the interpreter and plotting a graph for you:
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=http%3A%2F%2Fpm1.aminoapps.vertvonline.info%2F7617%2F347805917515b27bf7f1bcd5ab6dc1e493c4fdf5r1-1592-839v2_hq.jpg)
As expected, the falling speed actually stabilizes around (converges to) a certain value:... Roughly, 35m/s, which is around 126km/h.
----------------------------------------------------
Conclusion
Euler's method is a very powerful algorithm that can numerically solve most differential equations, linear or not. Taylor expand any derivative of your function until the last derivative, then find the latter by using the differential equation.
As you can see, there is a certain domain of validity to the fundamental theorem of engineering. You don't have to discard it right away, but PLEASE do try to stay as close as possible to the domain of validity if you want to have consistent results.
![An improved version of the fundamental theorem of engineering?-[C]Ah, ****, here we go again.
[BC]The fundamental theorem of](https://image.staticox.com/?url=https%3A%2F%2Fpa1.aminoapps.vertvonline.info%2F7617%2F7b195f820b8d28d924fa5313a8c2b72251d4d17fr1-480-360_hq.gif)
Thank you for listening to my long ass rant, and see you next time.
===============================
Post scriptum
The non differential equation here is actually solvable. I'll discuss the proof in another post.
Comments (5)
Fantastic post as always :yum:
Thank you!
I really enjoyed this post :blush: what an amazing post
Thanks! Since it's summer for me, I'll be posting more open problems, maths, physics and computer science. I hope you'll look forward to it!
Reply to: Rion
Yes i surely look forward to it. to also relax since it's summer!