Kepler’s Impossible Equation

Why is such a simple equation so difficult to solve? Head to www.kiwico.com... and use code WELCHLABS for 50% off your first monthly club crate or 20% off your first Panda Crate.
Welch Labs Book Ships December 2024: www.welchlabs....
Welch Labs Posters: www.welchlabs....
How the Bizarre Path of Mars Reshaped Astronomy: • How the Bizarre Path o...
Support Welch Labs on Patreon! / welchlabs
Special thanks to Patrons: Juan Benet, Ross Hanson, Yan Babitski, AJ Englehardt, Alvin Khaled, Eduardo Barraza, Hitoshi Yamauchi, Jaewon Jung, Mrgoodlight, Shinichi Hayashi, Sid Sarasvati, Dominic Beaumont, Shannon Prater, Ubiquity Ventures, Matias Forti, Brian Henry, Tim Palade, Petar Vecutin, Nicolas baumann
Learn more about WelchLabs!
www.welchlabs.com
TikTok: / welchlabs
Instagram: / welchlabs
REFERENCES
Colwell, P. (1993). Solving Kepler's Equation Over Three Centuries. United Kingdom: Willmann-Bell.
Needham, T. (1997). Visual Complex Analysis. United Kingdom: Clarendon Press.
Bate, R. R., Mueller, D. D., White, J. E. (1971). Fundamentals of Astrodynamics. Egypt: Dover Publications.
Vallado, D. (2001). Fundamentals of Astrodynamics and Applications. Netherlands: Springer Netherlands.
Borghi R. On the Bessel Solution of Kepler’s Equation. *Mathematics*. 2024; 12(1):154. doi.org/10.339...
Tom Archibald, Craig Fraser, Ivor Grattan-Guinness, The History of Differential Equations, 1670-1950. Oberwolfach Rep. 1 (2004), no. 4, pp. 2729-2794
Francisco G. M. Orlando, C. Farina, Carlos A. D. Zarro, P. Terra; Kepler's equation and some of its pearls. Am. J. Phys. 1 November 2018; 86 (11): 849-858.
Arthur A. Rambaut, M.A. A Simple Method of obtaining an Approximate Solution of Kepler's Equation. *Monthly Notices of the Royal Astronomical Society*, Volume 50, Issue 5, March 1890, Pages 301-302.
Ben Coleman. How to Find the Taylor Series of an Inverse Function. randorithms.co...
17 Recent papers on Kepler’s equation can be found in references 2-16 here: www.mdpi.com/2...

[@WelchLabsVideo]

Now is a great time to have a fresh look at this video’s sponsor KiwiCo - head to www.kiwico.com/welchlabs and use code WELCHLABS for 50% off your first monthly club crate or 20% off your first Panda Crate.

[@stopper0203]

Will your book be shipping outside America?

[@Anon54387]

I just factor out the E, set 1-e*sin to zero and call it a day.

[@TheDavidlloydjones]

@@Anon54387 But be sure to make allowance for the fact that the kiwis are all upside down. Since -0 = 0, this is very often a wash. Also watch out for leap years. 2400 CE is going to be another weird one.

[@revenevan11]

As a space nerd who loves this sort of math, I'm astounded that I never realized how central Kepler's equation was to the historical development of calculus!

[@revenevan11]

And of complex analysis, more broadly I guess.

[@Christopher-e7o]

X,2x+5=8')

[@akaHarvesteR]

Oh man this brings back memories. KSP used three different methods to compute E. The Newton-Raphson method actually turned out to be very unreliable for several cases, and would often blow up. What worked pretty well was a binary partitioning method, where you start with a large correction and half the search space at each iteration. It was slower to converge, but almost always found a decent solution. There was also a dedicated solver for extreme eccentricities, and of course another one for hyperbolic cases.

[@Tuned_Rockets]

Woah it's the man himself!

[@VeganSemihCyprus33]

You don't need to solve it, and that's how it is solved 👉 The Connections (2021) [short documentary] ❤

[@ionre24]

I was wondering how you knew this and then I saw your username. I was working on a 2d orbital space game a few years ago and I wondered how KSP calculated orbits, and I guess now I know :)

[@Bizob2010]

Woah :D. Would you approach it any differently today?

[@jmr5125]

Do you remember under what circumstances Newton's method failed to converge for elliptical orbits? I know that hyperbolic orbits are problematic, especially with high e, but I'm surprised that Newton's method would ever fail for elliptical orbits.

[@marcorademan8433]

Keplers algorithm resembles a control system closed loop in discrete steps, where the error signal is fed back with a proportional gain of 1. Some gains can make the closed loop unstable. This is an early version of a P-controller! You can analyse stability in terms of the z-transform of the open and closed loops.

[@logickedmazimoon6001]

so you can make a block diagram for this equation interesting.

[@paramdandekar562]

yeah lol when i saw that part that's exactly what I thought of too! now this has me wondering if newton's method can be expressed in control system terms too

[@jaspervandenbosch3838]

​@paramdandekar562 It would be a PD method, since it involves both the output value at the time step (P) as well as its derivative (D)

[@paramdandekar562]

@@jaspervandenbosch3838 not exactly though right? PD is when U(s)/E(s) = Kd*s + Kp, whereas Newton's method has the derivative divided by the value?

[@jaspervandenbosch3838]

@paramdandekar562 Yeah, you're right it's not really a PD. I think it should be possible to express in terms of a control scheme though

[@QuantumHistorian]

3:00 That's slightly off: the orbit is a parabola only at exactly e=1. For greater values of e it's a hyperbola. Not that having the right name for the right conic section is all that important here.

[@PhoenixEditz69]

Well coincidentally, our teacher today began teaching about parabola and when he said for e=1, locus is parabola There i got confused because in the video he mentioned something else😅

[@isaacclark9825]

@@PhoenixEditz69 If I recall correctly, a completely different approach is used for hyperbolic orbits. Parabolic orbits, which require e is exactly 1 are as rare as perfect circular orbits. We don't use Kepler's equation for that either, because we don't have to solve for E (eccentric anomaly) to get the true anomaly. We get the true anomaly directly from a much simpler equation.

[@PhoenixEditz69]

@@isaacclark9825 yeah i also think in actual life parabolic trajectories are rare as nothing is perfect and hyperbolas are more relevant I dont know if im correct but this could also be the reason why we study hyperbolic geometry and spherical geometry at higher level I haven't seen many people talk about parabolic geometry specially being a high school student😅 Most of college stuff like multi variable calculus, differential equations, and astrophysics i have learnt through YT videos So finding right stuff sometimes become tricky Like why eccentricity can't be zero But rather it tends to zero🤔 Though visualy it seems correct but mathematically it's counter-intuitive

[@bobbyheffley4955]

​@@isaacclark9825e=0 for a circle

[@JustRob96]

@@isaacclark9825 ​​⁠​⁠​⁠”we get the true anomaly from a much simpler equation” - Yes, it’s on the screen at 09:07! Barker’s goes direct from true anomaly to time for parabolic orbits, with analytical solutions in both directions. Because the comets have such a high eccentricity, they can be approximated with a parabolic close to the Sun - hence Halley’s method mentioned at this point in the video For hyperbolic orbits you are right, we have a different equation, but it is similar in that it is also transcendental and requires similar numerical techniques to solve

[@andrewmccandliss399]

Using e for eccentricity is evil. Seeing a formula with e^n in it where e isn't euler's number shakes me to my core

[@VerrouSuo]

i do a lot of disk work in astrophysics and we use capital Σ as a parameter (surface mass density) and im still not used to it

[@andrewmccandliss399]

@@VerrouSuo That is weird, usually I see lowercase sigma for surface density. Using capital greek letters as parameters is very weird

[@VerrouSuo]

@ i always supposed its because lowercase σ is the stefan-boltzmann constant, but ive also always used σ_b for that anyway (incidentally its also the standard deviation, but while that comes up all the time when you analyze data, i never actually have to write it symbolically)

[@CepheusMappy]

Some use ε as eccentricity though...

[@andrewmccandliss399]

@@CepheusMappy Yes. That's what I've always seen used. It's an e without being the symbol for euler's number essentially.

[@TristanPopken]

I remember trying to solve this before knowing it was a famously difficuly/impossible problem, got to the exact form of M=E-esin(E), though with different variable names, and tried to use the taylor series to make an estimate. Worked well for low eccentricities, but lets not even talk about the high ones haha. Great video about this topic, so here is a comment for the algorithm :)

[@maxmn5821]

Similar thing happened to me with x+log(x)=0 For a while I refused to acknowledge the problem - this was just a boring heat exchanger, after all. Until wolfram alpha convinced me

[@andrewzhang8512]

@@maxmn5821 lambert w function moment

[@unflexian]

​@@maxmn5821similar experience, banging my head against a wall trying to invent the Lambert W function until i heard about it in a blackpenredpen video

[@Sluppie]

Same. I was working on a programming project that needed to model orbits. The idea was that you had a "time slider" that you could move back and forth like a scroll bar to make the planets move around the star. I also needed to have a rocket ship that could apply thrust at specific points in time to alter orbits and transfer between different planets and the star. I thought that E would be an easy solve. I thought... I eventually just ended up taking the idea of a binary search tree and applying it to E, so that I could find it out via brute force. It worked pretty fast since I was using a computer, so I didn't worry about it too much beyond that.

[@darrendrapkin4508]

@@Sluppie I think you have discovered that simple, stable, numerical methods are preferable to more advanced ones, on a computer. I know someone who took thousands of terms of the Taylor Series for sin(x) to examine it for convergence. He found that even the error on a sine wave is a sine wave!

[@cheesedmacaroni]

I have a beautiful solution to Kepler's equation, but it could not be contained within the margins of this comment section

[@thesmartstickguy1145]

I get that reference

[@sajjadakbar6649]

Good Now i will look for your proof

[@mrfarts5176]

Ou must understand that Einsteins equations but be run assuming a flat earth. Well done...

[@VeganSemihCyprus33]

You don't need to solve it, and that's how it is solved 👉 The Connections (2021) [short documentary] ❤

[@thesmartstickguy1145]

@@VeganSemihCyprus33 Lol. It's still an unsolved problem. Did you not watch the video?

[@-_Nuke_-]

I honestly think that this is one of the best channels in all of the internet, not just KZbin. I feel blessed everytime I venture in the minds of the geniuses of our past.

[@VeganSemihCyprus33]

You don't need to solve it, and that's how it is solved 👉 The Connections (2021) [short documentary] ❤

[@OrangeDurito]

There is one guy named Michael MacKelvie who makes videos about sports which is at par with with the research, narrative style, and animation of this channel. That is also brilliant. Check it out!

[@qadirtimerghazin]

While watching the part about Kepler’s numerical approach, I thought “That’s what the Newton(-Raphson) method is for…” Turns out Newton came up with it specifically for this equation!

[@VeganSemihCyprus33]

You don't need to solve it, and that's how it is solved 👉 The Connections (2021) [short documentary] ❤

[@oniondeluxe9942]

I experienced exactly the same thoughts. 😮 It’s so sad when I realize I have learned so many things through life, and still, my life is a failure.

[@MathwithMing]

(Edit: the last paragraph) Superb work!! It always baffles me that Lagrange inversion isn't a standard topic in Calculus courses in universities. Most equations in science and engineering are "impossible" to solve in the sense that no closed-form solutions exist. Yes, we have powerful computers to quickly obtain highly accurate numerical solutions, but an approximate, analytic solution gives a much better qualitative picture. (perhaps it is difficult to justify convergence at elementary level, i.e., without complex analysis, but the inversion technique is just too good to miss out!) And thank you very much for providing reference btw. I'd like to add one more: Theory of Functions of a Complex Variables, by Soviet mathematician Alexei. I. Markushevich, translated by Richard A. Silverman. It has an extensive treatment on the technique of inversion of power series. In particular, the power series solution to Kepler's equation is covered in Example 3, page 99, Volume II. Markushevich's treatise is also extremely thorough (more so than western classic such as Ahlfors) and is the standard reference on complex analysis in the Russian-speaking world. PS. For those who wish to learn about Lagrange inversion but without a background in complex analysis, see pp151-155 Art.125,126 of the book "an Elementary Treatise on the Differential Calculus; Containing the Theory of Plane Curves, with Numerous Examples" by Benjamin Williamson. It is an old book and is available on the Internet Archive.

[@tomerpeled8922]

I solved it actually. Will post proof in a sec

[@platymusPrime]

It’s true, he showed it to me last night

[@DevangPatil-i9f]

It's true my professor told me about him

[@Bombito_]

​@@platymusPrime¿In a dream with a Indian god?

[@ryanschneer]

Ok Fermat

[@maloxi1472]

It's been 676 seconds 😢

[@Rubrickety]

This Kepler series has been superb. (And let’s all take a moment to thank Grant Sanderson for making Manim open source! Its signature look and power is all over KZbin these days.)

[@SpectatorAlius]

There is no mystery here.. Th\ss equation, like many others in physics, has no *closed form* solution. But we live in the era of computers, so we can approximate is numerically, something Kepler could do only in simple cases.

[@janisir4529]

A function not having analytic solution is indeed a good reason why it's hard to solve.

[@VaraNiN]

16:52 Laplace such a good physicist he still attends lectures 4 years after his death :P

[@lethargogpeterson4083]

Now THAT is dedication to accuracy. Your ghost is running out of academy meetings to issue errata.

[@jannegrey]

And lot of dedication from Cauchy, given that he was abroad - banished from France at the time.

[@siddharth-gandhi]

Bro is single handedly bringing back the science in science education channels on YT, nothing but bravos to you! Always excited to see what you teach nexy

[@darrenstensland5301]

Well, with some help from Grant Sanderson and Maxim.

[@giannisr.7733]

WAIT A MINUTE. I remember watching your whole series on imaginary numbers back when i was 16 and I knew nothing about them. You sparked my joy for maths and now here i am 6 years later struggling with my elecrtical engineering degree :(

[@thai-cheese]

Watching this two weeks into a complex analysis course would add do much motivation

[@BernardLowe-v1w]

Oh, man! Believe it or not, I just started taking a course on Numerical methods recently and only 2 days before, I had a lecture on the Newton-Raphson method in class. AND, I just started learning python animations with Manim! That's how you make your videos, right? You chose the perfect time to drop this, lol. I am starting to see just how much code and time you need to devote to get these beautiful animations. Thanks for the great vid!

[@Altair705]

My favorite method to solve it is the Halley method. It's similar to the Newton-Raphson method but also takes into account the second derivative, which makes the convergence faster. Generally speaking, that method is not very popular because the extra calculations needed to get the second derivative tend to offset the faster convergence, but in this case the second derivative of "E-e.sin(E)" is just "e.sin(E)", so no extra calculation is needed. In the end, the Halley method appears to be more efficient in the case of the Kepler equation, which makes me wonder why it's not more popular.

[@benruniko]

I love your video on imaginary numbers. It helped me gain the intuition needed to work with them. I think one reason we may never find a perfect solution to Kepler’s equation is the fact there is no perfect formula to draw an ellipse. You can easily plot a circle with a single function but true ellipses… not to much. Not without multiple steps, iteration or some approximations. It would be very cool to discover one some day though

[@yan.weather]

So good to have these content… keeps studies exciting and broadens the context of material

[@markuswx1322]

Years ago when learning how to compute implied interest rates and future cash flows based on the present value of an investment, I had no idea how closely these were related to fundamental problems in mechanics and other fields. Thanks for a very informative video.

[@euromicelli5970]

Not explicitly mentioned I presume for simplicity reasons, the Newton-Raphson method relies on the derivative of the function (see if you can spot it 7:03), which of course Kepler did not have but Newton very much did. Newton-Raphson remains to this day a very practical and useful numerical method (although it’s not the only one and it doesn’t work well for every kind of equation).

[@isaacclark9825]

Wow. I have worked with this equation since I was a high school lad, even up until last year when I implemented a solution (for low eccentricitry orbits) in software last year I expected that the subject would interest only about four KZbin viewers, but you have done much better than that!!

[@SolarAlmanac]

Ahhh, thank you so much. If only I had known all this when I started writing my astronomy code in high school :) PS. Last year I found Kepler's Equation in a strange place… in a _gardening_ problem. Given a length of chicken wire: Form the chicken wire into a circular arc (less than 2π radians) and place the ends against a straight exterior wall of the house. How large a sector of the circle (call it 2θ) should I use to maximize the enclosed area? … Well, when θ − sin θ appeared in the derivative, I said "No no no! That's Kepler's equation and can't be solved explicitly." By this time I knew better than to try! So I graphed it, and π/2 looked good, so I plugged it in, and yada yada yada the answer is: Semicircle.

[@carultch]

I stumbled upon this equation when trying to solve the exact same thing, for a recent Walter Lewin problem about the transfer orbit to Mars. I wanted to make a Geogebra animation that showed the path of the rocket as a function of time, with Earth and Mars taking circular orbits at the same time. I derived the equation of this form, and realized there was no analytic solution, and implemented a first order use case of Newton's method to solve it.

[@zimriel]

For my part I was trying to backtrack Mars' orbit (which is elliptic) from a point where the planet and my spacecraft would meet. I succeeded... some of the time.

[@Calcprof]

Ah, playing with Gil Strang's Introduction to Applied Mathematics. I taught from this book many times in the 1990s.

[@aljawad]

Ahh, this brings back so many fond memories of my first experience in programming. Some 50 years ago when I started to plot cometary positions based on their orbital elements, I used a programmable calculator to hunt for the value of E, basically using iteration (my first Hewlett Packard calculator had only 49 steps). Eventually when I acquired my first computer (an Apple II), I wrote dedicated programs (initially in Basic, then Pascal) to tabulate an ephemeris of the positions. A fascinating book on the subject is Peter Colwell’s “Solving Kepler’s Equation”, which I highly recommend 👍🏼

[@liamhagan4434]

Hey man, really excited about the imaginary numbers book. Your series on imaginary numbers in 2016 was extremely impactful to me then and now. I still think about those videos often, and can’t wait to have this book in my possession!!

[@deer001]

After watching your complex number series, I felt heartbroken that you stopped posting.I didn't even know you were back until a short of yours popped up! So glad you are back! Take love❤

[@cameronspalding9792]

If we let f(E) = E-e*sin(E), then f(M-1)= (M-1)-e*sin(M)

[@maitland1007]

This was really great. The only thing Id want more of is more explanation of the history/derivation of the equation. Maybe another video?

[@WelchLabsVideo]

I cover some of this here! kzbin.info/www/bejne/hpnWlJ2iZdppiqs

[@maitland1007]

@@WelchLabsVideo Thanks.. just watched the first one.. my reaction is.. WOW!.. so amazingly well done! thanks!

[@kellymoses8566]

Veritasium has a good video about this.

[@swamihuman9395]

- I can't love this enough! :)

[@Nabla_Squared]

I love so much your videos, your style, explanations and the topics you choose are simply perfect. Excelent work!

[@ivoryas1696]

10:20 That... was pretty nice! Good point! I've honestly want to someday construct a mechanical integrator for similar purposes. 😅

[@ivicasmolic2221]

Small remark: in order to get series in e via Lagrange inversion theorem, mentioned in the video cca 13:00, it is convenient to choose f(E) = (E-M)/sin(E) (note that f is analytic around M and f'(M) = 1/sin(E) is not zero) and look at the equation f(E) = e.

[@fg786]

5:13 Something is off with your line for M = 77,7°. It intersects the white curve at E >120° but when you zoom in the whole graph is shifting. It looks too low to be the value for E close to 90°. So it's up to you to figure out what's wrong here.

[@ErwanLeGac]

You're reminding me my first steps in programmation (in BASIC!), when I tried to make a program to calculate the position of the planets...

[@astronomianova797]

Prussing and Conway's Orbital Mechanics book suggests using the Laguerre algorithm. I used it years ago when I wrote programs for astrodynamics. The Laguerre method was meant for solving polynomials of degree n. Although Kepler's Equation is not a polynomial, you can just use n=4, 5 or 6 and run the iterations. I verified what was said in the book; it converges quicker than Newton's method and converges for any initial guess (a bad first guess in Newton's method may lead to no convergence).

[@officialjfendi]

Amazing video! I'm actually taking a grad-level orbital mechanics class right now and had to implement a bunch of these algorithms that you mentioned in the video, this explains it in a really nice and conceptually palatable way so nice job :D. Such a fun topic haha

[@phenixorbitall3917]

Tbh...Cauchy is a real beast! This guy was disrespectacularly smart - no kidding

[@yan.weather]

Disrespectacularly… I learnt a new word

[@phenixorbitall3917]

I made it up :) you're welcome

[@jedswift]

Fascinating history of a real problem being investigated by some truely awesome intellects. As an engineer, I just used dA/dt = C to estimate the position vs time on an orbital trajectory of a spacecraft in a geostationary transfer orbit from low Earth orbit. I just divided the orbit into 1º slices assuming straight lines for every slice resulting in an ellipse with 360 sides, made of triangles all with an easily calculatable areas summing to the total area of the eclipse. Linear interpolation on the little line segments worked just fine for positions not on even degrees. Turned out to be accurate within a second of the actual ~12 hr orbital period. All hail Excel :-).

[@dendarius9906]

Kepler is one of my favorite. His equations are remarkable. Dark matter seem to enforce both his law of equal distance over equal time and newton law.

[@zeFresk]

This was so beautiful I almost teared up. Thank you for posting this video.

[@EtotheFnD]

Is okay to admit i enjoyed watching this entire video, shared the video to a discord I am in, but have NO IDEA HOW THE MATHEMATICS WORK?

[@costbubbles8336]

ive been waiting for this video for so long after seeing the yt short

[@cparks1000000]

The Lagrange Inversion Formula looks super complicated but it's just Taylors Theorem at f(a) plus "calculating" the (higher) derivative(s) of the inverse function (i.e., prove that $(f^{-1})^{(n)}(f(a))= \lim_{x \to a} (d/dx)^{n-1}[(x-a)/(f(x)-f(a)]$)

[@Locateson]

Bessel functions bring back happy memories, I wrote a planetary system simulator some time in 2002 and used that stuff extensively.

[@iccuwarn1781]

Fantastic video! You made the convergence of Taylor series so clear and easy to understand! =)

[@StuartJdAntoine]

E simplify in a scalar matrix to -sin(e) since that returns standard deviation of a log variance for any index in E. If E is an empty set, it returns -1.

[@tune490]

This is an amazing video! Thank you Stephen :D

[@exponentmantissa5598]

Well done. I always like your stuff. - Retired Engineer

[@0cgw]

If confronted with this equation I would rewrite it as a fixed point equation: E=f(E) where f(E)=M+e sin E. For e1 (hyperbolic orbits), we have E=g(E) where g(E)= sin⁻¹ [(E-M)/e] and the same approach should work at least if e is sufficiently large.

[@butthardley5160]

I like your thinking

[@JustRob96]

Here’s a fun idea. Imagine taking the Newton-Raphson method, and upgrading it to use higher order derivatives too. In general I believe this is known as Householder iteration. The trouble with this is that in order to use it, you have to compute those higher order derivatives. The rate of convergence for each iteration is higher, but the workload per iteration increases as well. Sometimes it can be worth it. For example, there’s a 2014 paper by a man named Izzo who uses this to solve Lambert’s problem from the field of… orbital mechanics! Suppose you have two position vectors of an asteroid and the times of measurement, can you calculate what the asteroid’s orbit must be? That’s Lambert’s problem

[@zimriel]

Altair 705 here has pointed out that the second derivative is Halley's method, which here is just... e sinE.

[@qadirtimerghazin]

8:57 “high curvature causes Newton's method to overshoot the correct answer” I assume this also refers to 6 iterations, as it was for the Kepler’s method? I wonder how many more NR iterations are actually needed to converge? In fact, if it’s done manually without a computer, it seems to me that it would be easier and faster to do a few more NR steps than compute all those really complicated series…

[@jmr5125]

@qadirtimerghazin Yes, Newton's method will generally converge -- sometimes it requires extra iterations. However, based on other comments in this video, there are specific enough values for which Newton's method *won't* converge, no matter how many iterations you use. When / how this occurs I do not know although I strongly suspect it has to do with eccentricity values close to (but below) 1 and a very poor initial guesses for E.

[@qadirtimerghazin]

@@jmr5125 I didn’t have time to investigate this in detail, but our future overlord GPT 4o said that “Newton-Raphson method is widely used for solving Kepler’s equation because it converges very quickly for small and moderate eccentricities. For high eccentricities and problematic initial guesses, however, alternative or hybrid approaches might be necessary to ensure reliable convergence.” As a computational chemist, I’ve dealt with many cases when the NR method frustratingly just couldn’t converge trying to find a critical point on a highly multi-dimensional molecular potential surfaces, but it’s quite interesting that even this seemingly simple function can be challenging for the NR method…

[@jonquil3015]

Unbelievably good presentation!

[@pizza8725]

I might of done a mistake that simplified the problem bc I just solved E based on M The formula is E=-i*ln(-(2W(1÷2-e^(iM)÷2)-1) Where W is the W lambart function and i is the imaginary unit

[@patriciodasilva7902]

Using Newton's method, we solved Kepler's equation for the given mean anomaly and eccentricity . The resulting value for the eccentric anomaly is approximately . This iterative method finds by gradually refining an initial estimate until it converges within a set tolerance, giving us an accurate solution for the position of the planet in its orbit.

[@josephyoung6749]

Didn't understand this the first time through but will watch again

[@Elektrolite111]

It would be nice to see a deep dive into the radius of convergence

[@pikkisir4664]

I love your videos a lot

[@unflexian]

zetamath video about analytic continuation is focused on the subject of the radius of convergence in the complex plane, and it's fantastic!

[@bobnobrain]

Yay it's finally out!!

[@regulus8518]

please do a video like this on the history behind runge kutta methods ode45 etc for differential equation solution ... it was very interesting to learn context behind newton raphson and how and why it was developed, it is used in electrical engineering for something called load flow analysis and allows for things like economic dispatch and maintainence on the grid without swictching the entire thing off

[@xanterrx9741]

Great video , thanks for making it

[@icenarsin5283]

Awesome video.... Thank you!

[@itzmetanjim]

4:15 it does work, the absolute value of the error is reducing , even though the sign is flipping. this is still converging (much more slowly though) edit: he said 'does not work *at all* "

[@sage5296]

I remember I once tried to program in planetary motion based primarily on Kepler's law and subsequently gave up after several hours of trying. Fun times lol

[@witsued]

Edmond Halley, like rally. J.K. Rowling, like bowling.

[@RickyMud]

I swear the algo has been hyping up this equation recently, unless that you guys doing it

[@lorenwilson8128]

Hailey expanded Newton's method to use the second derivative to correct for under or overshoot. It will converge significantly faster.

[@Magicraft13]

Hi! Amazing video, I was wondering if an EU-shipping delivery will be available for your shop in a close future, I really want to buy some items :). Thanks for your astonishing work!

[@michaelpieters1844]

Amazing video! Where did you get all this information? Which sources did you consult or did you go to the original works of Kepler, Newton, Laplace ... themselves?

[@WelchLabsVideo]

Thanks! Sources in description.

[@harryh4398]

I thought this would be impossible for me to follow but apart from the first 22 mins and 41secs, it really wasn't too bad.

[@wjalp]

Kepler's method is very smart! 6:00 :))

[@amorphant]

The music really adds gravitas.

[@ruperterskin2117]

Cool. Thanks for sharing.

[@rb8049]

A Pade expansion including both poles and zeros is more general than a Taylor series expansion.

[@mrhatman675]

It isn t really

[@FadkinsDiet]

Padé approximants are cool but not really more applicable except in special cases.

[@mumujibirb]

4:10 I've seen people use a guess and check approah, i.e. they try a value and see if it's too small or large, then guess again So i think this is really just a natural extension, after all, -since sin(x) = x for small values-

[@PrimordialOracleOfManyWorlds]

fantastically fascinating.

[@alliknowissuperposition2848]

Finally, after a lot of teaser 😅

[@RussellJones1961]

I think you’ll find that Sir Edmund Halley’s name was actually …. Halley.

[@VictorOrji-mz1jw]

Is it me or does the graph at 4:58 like like a cumulative frequency curve

[@-Yousof-]

Great content!

[@rgerk]

Maybe the solution will come if you use rational trigonometry, changing angles for spreads and distances for quadrances, avoiding square roots and transcendental functions.

[@saxtant]

It's equivalent to conservation of angular momentum, ellipse or circle, the energy is the same.

[@tmst2199]

He doesn't say that the reason for determining the area swept out in the circle is because we don't know the elepticity of the ellipse.

[@Kaviranghari]

ok this is crazy i was talking about the newton raphson method all day and my teacher was and this video has popped up (by the way i was talking about it because i noticed it converges very quickly )

[@davidrandell2224]

“The Final Theory: Rethinking Our Scientific Legacy “, Mark McCutcheon for proper physics including the CAUSE of gravity, electricity, magnetism, light and well.... everything. Without the cause of gravity these models will remain incomplete. Newton/Einstein etc become obsolete mathematical projections. The Geometric Orbital Equation is: v^2R=K. K is the constant for different planets/sun’s etc.R=radius; v=velocity.

[@Emry11]

In your imaginary numbers series you start with a false premise that is so rampant among academics and it is very sad and disappointing. Basically, there exist nothing called i=sqrt(-1) That is a non-existent definition and nowhere in math history any reliable source has ever defined it. The true and unique definition is i^2 = -1 and from that no one can possibly deduce i= sqrt(-1). To make this clear once for all, let's (just for the sake of argument!) assume that i=sqrt(-1) exist (again, it doesn't but we pretend it does!). Then we will have: -1 = i^2 = i*i = sqrt(-1) * sqrt(-1) = sqrt( -1 * -1) = sqrt (+1) = +1 which wrongly implies -1 = +1 Besides, try the Euler's equation itself: e^(ix) = Cos(x) + i Sin(x) and try to substitute i with sqrt(-1) and see if you get the same result which of course you don't because i is the complex number and can't be substituted by any other definition. Another obvious example would be e^(i * Pi) = -1. Try to substitute i with sqrt(-1) in that and see if e^(sqrt(-1) * Pi) will be equal to -1 which you won't be able to show. I hope that you and all others are now clear about this very common and unfortunate mistake and actually mis-definition of the complex number "i" and won't use that false definition from now on. Good Luck!

[@alexbennie]

All the OGs of Real Analysis and Complex Analysis getting together on a question about horoscopes! Being able to tell precisely when Uranus aligns with Venus is useful. Especially if you are a Certified Genius (together or without the associated stereotypes) in a club... With your crush that's into, for example: Horoscopes. I do see a lot of mentions of Kepler's math in the more "Sciencey" Horoscope videos though... So it's plausible that he and they were not unfamiliar. Sooo, yah... He tried to measure the orbits of Mars and whatever, to impress a chick. You can't tell me that at least one of these Greats: Cauchy, Laplace, Lagrange, Bessel, et al... tried to pick up a chick by knowing the science behind her 'hobby'?

[@souldrip2000]

at 16:47... it would have been quite difficult for Laplace to hear Cauchy's 1831 lecture at the French Academy of Sciences. (Laplace died in 1827.)

[@primenumberbuster404]

Cauchy is like that one op anime character that basically does everything with style.

[@thamiordragonheart8682]

He probably got the decade wrong while he was talking or something simple like that.

[@manhhoanguc837]

Just use Bisected Direct Quadratic Regula Falsi. For M in range 0 to π, set the search range [a,b] of E to [0,π] and the start point to M. The same apply for the range [π,2π]. This has the same quadratic convergence rate as the N-R method, but will always converge for e very close to 1.(Even if e is 0.9998 or closer depends on the floating point precision.) If you tweak somewhat, it will also works for the hyperbolic case e>1, and for very large e > 50, and M > 10π or M

[@FadkinsDiet]

That requires many more steps to get the error down to acceptable levels, even if it is theoretically quadratic in the limit

[@manhhoanguc837]

@@FadkinsDiet Nope, each step only requires one additional evaluation of the current function, without having to evaluate any of the derivatives like the harley families. For accuracy up to 12-13 decimal points, this method only requires 5-7 func evaluations.

[@einfisch-z9f]

I understood nearly nothing but this sounds very interesting and important

[@ExistenceUniversity]

I don't understand why people are so proud of not understanding such that it needs to be shared. Did you try pausing and following along with a pen and paper? Why openly admit you don't understand?

[@DavyCDiamondback]

It's always aggravating when classical physics problems are unsolvable, because if classical is unsolvable, what hope is there for relativistic and quantum physics?

[@Sanchuniathon384]

HOLD UP, where do we find how to use Bessel functions to model light diffracting through DNA strands? I am so down for this.

[@RyanLynch1]

2:13 hey that's my birthday!

[@elektr1x7881]

mine too!

[@tgrey_shift..mp334]

The error correction portion reminds me of Bayes Theorem! Maybe its related to that some way?