Kepler's equation M = E − e·sin E — derivation, eccentric anomaly: Worked examples

3.2.14 · D3 · Physics › Orbital Mechanics & Astrodynamics › Kepler's equation M = E − e·sin E — derivation, eccentric an

Neeche sab kuch radians mein hai jab tak koi step explicitly convert na kare. $M$ = mean anomaly (uniform clock), $E$ = eccentric anomaly (auxiliary circle ke centre se measure hoti hai), $e$ = eccentricity ( $0 \leq e < 1$ ), $a$ = semi-major axis, $T$ = orbital period, $t_{p}$ = time of perihelion passage.

True anomaly $ν$ (Ex 1 mein use hoti hai) planet ki actual angular position hai jo Sun (focus) se dekhi jaati hai — woh actual direction jis par tum telescope point karoge. Ise $E$ se convert kiya jaata hai half-angle relation se: $tan \frac{ν}{2} = \frac{1 + e}{1 - e} tan \frac{E}{2},$ jo us companion page par radius ke do forms, $r = a (1 - e cos E)$ aur $r = \frac{a ( 1 - e ^{2} )}{1 + e cos ν}$ , se derive ki gayi hai. Hum ise yahan sirf use karenge; agar koi bhi word unfamiliar lage, pehle parent padhlo.

Scenario matrix

Kepler's equation ek formula jaisi lagti hai, lekin behaviour is baat par depend karke wildly change hoti hai ki planet kahan hai aur orbit kitni squashed hai. Yeh raha har case-class jo is topic mein aa sakta hai:

#	Case class	Kya special hai	Covered by
A	Degenerate: circle $e = 0$	Koi correction term nahi, $M = E = ν$	Ex 1
B	Turning points $E = 0, π$	$sin E = 0$ , clock = geometry	Ex 2
C	Quadrant I $0 < E < \frac{π}{2}$ (Sun ki taraf aa raha hai, fast chal raha hai)	$sin E > 0$ , $M < E$	Ex 3
D	Quadrant II $\frac{π}{2} < E < π$	$sin E > 0$ abhi bhi, lekin $y$ -peak ke baad	Ex 3 (table)
E	Quadrants III & IV $π < E < 2 π$ (door ja raha hai, slow chal raha hai)	$sin E < 0$ , isliye $M > E$	Ex 4
F	Easy direction $E \to M \to t$	Koi iteration nahi chahiye	Ex 5
G	High eccentricity $e \to 1$	Slow convergence, bura seed	Ex 6
H	Real-world word problem (comet time-of-flight)	Math ke around physics wrap karo	Ex 7
I	Exam twist (given $r$ , find $t$ )	Radius se $E$ reverse-engineer karo	Ex 8
J	Sanity/units traps (degrees-vs-radians, sign of correction)	Wrong-unit inputs se nonsense aata hai	Ex 9

Ab hum har cell hit karenge.

Ex 1 — Cell A: degenerate circle

Worked example

$e = 0$ , $M = 2.3$ rad. $E$ , $r$ , $ν$ find karo. Forecast: Koi squash nahi, toh ellipse hi circle hai. Guess karo teen anomalies ka kya hoga.

Step 1. $e = 0$ ke liye Kepler's equation likho: $M = E - 0 \cdot sin E = E$ . Yeh step kyun? Pehle input substitute karne se pata chalta hai ki transcendental term bach bhi rahi hai ya nahi. Yahan woh vanish ho jaati hai — equation $E = M$ tak collapse ho jaati hai.

Step 2. Toh $E = 2.3$ rad, koi iteration nahi. Yeh step kyun? Circle mein ek hi centre = focus hota hai, isliye "Sun ke paas fast chalo" wala effect khatam; clock aur geometry identical hain.

Step 3. Radius: $r = a (1 - e cos E) = a (1 - 0) = a$ . Yeh step kyun? Circle par har point centre se $a$ distance par hota hai — formula ko constant return karna chahiye, aur karta bhi hai.

Step 4. True anomaly (introduction mein diye half-angle relation se): $tan \frac{ν}{2} = \frac{1 + 0}{1 - 0} tan \frac{E}{2} = tan \frac{E}{2}$ , isliye $ν = E = 2.3$ .

Verify: Teeno anomalies $2.3$ rad ke barabar hain aur $r = a$ — exactly uniform circular motion. ✓ Yeh limiting-case anchor hai: agar tumhara general code $e \to 0$ par yeh reproduce nahi karta, toh woh galat hai.

Ex 2 — Cell B: turning points

Worked example

$e = 0.6$ . $E = 0$ (perihelion) aur $E = π$ (aphelion) par $M$ aur $r$ evaluate karo. Forecast: Major axis ke dono ends par planet ka momentarily koi "sideways" $y$ -motion nahi hota. Guess karo correction $- e sin E$ wahan kya karta hai.

Step 1. Perihelion $E = 0$ : $M = 0 - 0.6 sin 0 = 0$ . Yeh step kyun? $sin 0 = 0$ correction ko khatam kar deta hai, isliye $M = E = 0$ : clock aur geometry shuruaat mein synchronized hain.

Step 2. $r = a (1 - 0.6 cos 0) = a (1 - 0.6) = 0.4 a$ — closest approach. Yeh step kyun? Perihelion matlab focus ke sabse paas; $a (1 - e)$ textbook ki minimum distance hai, jo hamare $E$ -form ko confirm karta hai.

Step 3. Aphelion $E = π$ : $M = π - 0.6 sin π = π$ . Yeh step kyun? $sin π = 0$ phir se — door wale turning point par bhi synchronization.

Step 4. $r = a (1 - 0.6 cos π) = a (1 + 0.6) = 1.6 a$ — sabse door.

Verify: $r_{peri} + r_{aph} = 0.4 a + 1.6 a = 2 a$ = major axis, jaise kisi bhi ellipse ke liye hona chahiye. ✓ Aur $M = E$ dono ends par parent ki baat confirm karta hai ki "clock aur geometry turning points par agree karte hain."

Ex 3 — Cells C & D: quadrants I aur II ( $sin E > 0$ , isliye $M < E$ )

Neeche ki figure geometry ko concrete banati hai. Dashed circle auxiliary circle hai (radius $a$ ). Kisi chosen eccentric anomaly $E$ ke liye, point $Q = (a cos E, a sin E)$ us circle par centre $O$ se angle $E$ par baitha hai; ise squash factor $b / a$ se seedha neeche drop karne par planet $P = (a cos E, b sin E)$ ellipse par aata hai. Jab planet upper half mein ho ( $0 < E < π$ ), toh $Q$ axis ke upar hai, isliye $sin E > 0$ , isliye correction $- e sin E$ negative hai aur $M < E$ .

Figure — Kepler's equation M = E − e·sin E — derivation, eccentric anomaly

Figure: auxiliary-circle point $Q$ (red) angle $E$ par, planet $P$ (black square) uske neeche, aur do half-orbit cases. Upper half → $sin E > 0$ → $M < E$ ; lower half → $sin E < 0$ → $M > E$ . Red radius $O Q$ dekho: axis ke upar uski height exactly $a sin E$ hai, woh quantity jise correction term multiply karta hai.

Worked example

$e = 0.3$ . $E = 0.8$ rad (quadrant I) aur $E = 2.6$ rad (quadrant II) se $M$ compute karo. Forecast: Kaun sa quadrant zyada bada gap $E - M$ produce karta hai? Compute karne se pehle guess karo.

Step 1. $E = 0.8$ (quadrant I): $M = 0.8 - 0.3 sin 0.8 = 0.8 - 0.3 (0.71736) = 0.58479$ . Yeh step kyun? Yeh direct forward map hai. Figure mein, $E = 0.8$ ka $Q$ quadrant I mein neeche wala red dot hai; uski height $a sin 0.8$ woh hai jo subtract hoti hai. $sin 0.8 > 0 \Rightarrow M < E$ : uniform clock kam read karta hai geometric angle se, kyunki fast-moving planet clock ki expectation se pehle hi nikal gaya hai.

Step 2. $E = 2.6$ (quadrant II): $M = 2.6 - 0.3 sin 2.6 = 2.6 - 0.3 (0.51550) = 2.44535$ . Yeh step kyun? Abhi bhi upper half, abhi bhi $sin E > 0$ , abhi bhi $M < E$ — lekin $sin 2.6 < sin 0.8$ , isliye gap $E - M$ chhota hai yahan. Figure mein yeh $Q$ circle ke top ke baad hai, uski height wapas neeche aa rahi hai. Correction $E = π /2$ ke paas peak karti hai jahan $sin E = 1$ (red circle ka top).

Step 3. Poore $(\frac{π}{2}, π)$ range mein quadrant II ke behaviour ko tabulate karo promised trend dekhne ke liye — gap $E - M = 0.3 sin E$ steadily $0$ ki taraf shrink hoti hai jaise $E \to π$ :

$E$ (rad, quadrant II)	$sin E$	$M = E - 0.3 sin E$	gap $E - M$
$1.8$	$0.97385$	$1.50785$	$0.29215$
$2.2$	$0.80850$	$1.95745$	$0.24255$
$2.6$	$0.51550$	$2.44535$	$0.15465$
$3.0$	$0.14112$	$2.95766$	$0.04234$

Yeh step kyun? Yeh cell D poori tarah se hai: yeh dikhata hai ki quadrant II mein correction abhi bhi positive-signed on $sin E$ hai (isliye $M < E$ ) lekin monotonically decreasing hai, aphelion $E = π$ par vanish hoti hai.

Step 4. Dono headline points compare karo: $E = 0.8$ par gap $0.21521$ hai vs $E = 2.6$ par $0.15465$ . Yeh step kyun? Forecast confirm karta hai — quadrant I ( $sin E = 1$ crest ke zyada paas) mein correction badi hai.

Verify: Dono headline gaps exactly $0.3 sin E$ ke barabar hain: $0.3 sin 0.8 = 0.21521$ ✓ aur $0.3 sin 2.6 = 0.15465$ ✓. Har table row $0.3 sin E$ ke barabar hai aur $0$ ki taraf decrease karti hai. Correction ka sign poore upper half mein negative hai. ✓

Ex 4 — Cell E: quadrants III & IV ( $sin E < 0$ , isliye $M > E$ )

Worked example

$e = 0.3$ . Given $M = 4.0$ rad (orbit ka lower half), Newton–Raphson se $E$ solve karo. Forecast: Lower half mein planet receding hai aur slow chal raha hai, isliye clock geometry se aage hona chahiye: guess karo $E > M$ hai ya $E < M$ .

Step 1. Seed $E_{0} = M = 4.0$ . Yeh step kyun? Moderate $e$ ke liye, $E \approx M$ ek excellent start hai; Newton yahan se quadratically converge karta hai.

Step 2. Newton step $E_{k + 1} = E_{k} - \frac{E _{k} - e sin E _{k} - M}{1 - e cos E _{k}}$ . $f (4.0) = 4.0 - 0.3 sin 4.0 - 4.0 = - 0.3 (- 0.75680) = 0.22704$ . $f^{'} (4.0) = 1 - 0.3 cos 4.0 = 1 - 0.3 (- 0.65364) = 1.19609$ . $E_{1} = 4.0 - 0.22704/1.19609 = 3.81017$ . Yeh step kyun? $sin 4.0 < 0$ (lower half), isliye correction $E$ ko $M$ se neeche push karti hai. Yeh forecast ka jawab deta hai: $E < M$ yahan.

Step 3. Phir iterate karo: $f (3.81017) = 3.81017 - 0.3 sin (3.81017) - 4.0 = - 0.00420$ ; $f^{'} (3.81017) = 1.22470$ . $E_{2} = 3.81017 - (- 0.00420) /1.22470 = 3.81360$ . Yeh step kyun? Ek aur Newton step; residual sharply shrink ho gayi hai, convergence confirm ho rahi hai.

Step 4. Ek aur step deta hai $E \approx 3.81371$ (5 dp tak stable).

Verify: Converged value wapas plug karo: $3.81371 - 0.3 sin (3.81371) = 3.81371 - 0.3 (- 0.62097) = 3.81371 + 0.18629 = 4.00000$ ✓. Aur indeed $E = 3.814 < M = 4.0$ , kyunki $sin E < 0$ ne $M > E$ banaya. ✓

Ex 5 — Cell F: easy direction $E \to M \to t$

Worked example

$e = 0.5$ , $E = 3 π /4$ rad, period $T = 200$ days. Perihelion ke baad ka time find karo. Forecast: Do directions mein se ek ( $E$ se $M$ find karna, ya $M$ se $E$ find karna) ko koi iteration nahi chahiye. Guess karo kaun sa — aur kyun.

Step 1. $E = 3 π /4 = 2.35619$ rad. $M = E - e sin E = 2.35619 - 0.5 sin (2.35619)$ . $sin (3 π /4) = \frac{2}{2} = 0.70711$ , isliye $M = 2.35619 - 0.5 (0.70711) = 2.00264$ rad. Yeh step kyun? $E \to M$ jaana ek direct substitution hai — Kepler's equation pehle se $M$ ke liye solve hai. Sirf reverse ( $M \to E$ ) transcendental hai. Yeh forecast confirm karta hai: forward free hai.

Step 2. Clock-angle ko time mein convert karo: $t - t_{p} = \frac{M}{2 π} T = \frac{2.00264}{6.28319} \cdot 200 = 63.746$ days. Yeh step kyun? $M$ hi orbit ka woh fraction hai jo angle ke roop mein express hua hai; $2 π$ se divide karne par elapsed period ka fraction milta hai.

Verify: Orbit ka fraction $= 2.00264/6.28319 = 0.31873$ , aur $0.31873 \times 200 = 63.746$ days $<$ half a period (100 days) — consistent hai, kyunki $E = 3 π /4$ aphelion se pehle hai. ✓ Units: (rad)/(rad) × days = days. ✓

Ex 6 — Cell G: high eccentricity, slow convergence

Worked example

$e = 0.9$ (bahut stretched orbit), $M = 0.3$ rad. $E$ solve karo. Forecast: $e = 1$ ke paas term $e sin E$ bahut bada hota hai, aur perihelion ke paas orbit sharply curved hoti hai. Guess karo kya naive seed $E_{0} = M$ abhi bhi achha hai.

Step 1. Seed $E_{0} = M = 0.3$ try karo. Newton: $f (0.3) = 0.3 - 0.9 sin (0.3) - 0.3 = - 0.9 (0.29552) = - 0.26597$ . $f^{'} (0.3) = 1 - 0.9 cos (0.3) = 1 - 0.9 (0.95534) = 0.14020$ . $E_{1} = 0.3 - (- 0.26597) /0.14020 = 2.19700$ . Yeh step kyun? Chhota denominator $f^{'} \approx 0.14$ (kyunki $e cos E \approx e$ perihelion ke paas) pehle Newton jump ko wildly overshoot kara deta hai — high $e$ par slow/erratic convergence ka symptom. Yeh forecast ka jawab deta hai: naive seed yahan kharab hai.

Step 2. Iterate karte raho. $f (2.19700) = 2.19700 - 0.9 sin (2.19700) - 0.3 = 2.19700 - 0.9 (0.80834) - 0.3 = 1.17949$ ; $f^{'} (2.19700) = 1 - 0.9 cos (2.19700) = 1 - 0.9 (- 0.58871) = 1.52984$ . $E_{2} = 2.19700 - 1.17949/1.52984 = 1.42604$ . Yeh step kyun? Newton recover karta hai kyunki overshoot ke baad derivative healthy hai; yeh ab root ki taraf wapas chalti hai.

Step 3. Jaari rakho: $f (1.42604) = 1.42604 - 0.9 sin (1.42604) - 0.3 = 1.42604 - 0.9 (0.98964) - 0.3 = 0.23536$ ; $f^{'} (1.42604) = 1 - 0.9 cos (1.42604) = 1 - 0.9 (0.14346) = 0.87089$ ; $E_{3} = 1.42604 - 0.23536/0.87089 = 1.15579$ . Yeh step kyun? Residual shrink ho rahi hai lekin abhi fast quadratic regime mein nahi — high $e$ genuinely zyada steps maangta hai.

Step 4. Continue karo (Newton ab quickly converge karta hai): $E_{4} = 1.11115$ , $E_{5} \approx 1.10991$ , aur yeh $E \approx 1.10991$ (5 dp) par settle ho jaata hai. Yeh step kyun? Yahan ~5 steps lage (ek wild pehli jump ke saath) vs $e = 0.3$ ke liye ~3 clean steps — cell G ka practical lesson: high $e$ ko better seed chahiye (jaise $E_{0} = π$ ) ya damped Newton step.

Verify: Converged root wapas plug karo: $1.10991 - 0.9 sin (1.10991) = 1.10991 - 0.9 (0.89694) = 1.10991 - 0.80725 = 0.30266 \approx 0.30$ ; refine karne par $E = 1.10971$ milta hai aur $1.10971 - 0.9 sin (1.10971) = 0.30000$ exactly. ✓ Key qualitative lesson — pehla Newton step overshoot karta hai aur high $e$ par convergence slow hoti hai — hold karta hai. ✓

Ex 7 — Cell H: real-world comet time-of-flight

Worked example Ek comet ka

$e = 0.7$ hai, period $T = 76$ years. Woh kahan hai (perihelion ke baad time ke roop mein) jab uska eccentric anomaly $E = π /3$ hai? Uski distance bhi $a$ ke units mein do. Forecast: Guess karo usne apne period ka zyada time bitaya hai ya kam, yeh dekhte hue ki $E = π /3$ ek full turn ka ek-shashthaansh hai.

Step 1. $M = E - e sin E = \frac{π}{3} - 0.7 sin (\frac{π}{3}) = 1.04720 - 0.7 (0.86603) = 1.04720 - 0.60622 = 0.44098$ rad. Yeh step kyun? Forward direction, direct. Note karo $M = 0.441 < E = 1.047$ : comet, Sun ke paas fast chal raha hai, geometrically uniform clock se "aage" hai. Isliye usne period ka $E /2 π$ se kam hissa bitaya — forecast confirm hota hai ( $T /6$ se kam).

Step 2. Time: $t - t_{p} = \frac{M}{2 π} T = \frac{0.44098}{6.28319} \times 76 = 5.3341$ years. Yeh step kyun? Ex 5 jaisi clock-to-time conversion; physics ab "years past perihelion" hai.

Step 3. Distance: $r = a (1 - e cos E) = a (1 - 0.7 cos \frac{π}{3}) = a (1 - 0.7 \times 0.5) = a (1 - 0.35) = 0.65 a$ . Yeh step kyun? Jab $E$ pata ho toh radius ek one-liner hai; comet Sun se apne semi-major axis ke $0.65$ hisse par hai.

Verify: $T /6 = 76/6 = 12.667$ years; hamara jawab $5.334$ years indeed kam hai ✓ (perihelion ke paas fast). Aur $0.65 a$ , $r_{peri} = a (1 - 0.7) = 0.3 a$ aur $r_{aph} = 1.7 a$ ke beech hai, isliye physically valid distance hai. ✓ Units: poore mein years. ✓

Ex 8 — Cell I: exam twist — given $r$ , time find karo

Worked example

$e = 0.4$ , $T = 100$ days. Planet distance $r = 1.2 a$ par hai aur receding hai (Sun se door ja raha hai). Perihelion ke baad ka time find karo. Forecast: $r = 1.2 a > a$ matlab planet $r = a$ locus (jo $E = π /2$ par hota hai) se aage hai, bahar ja raha hai. Guess karo $E$ upper half ( $0 < E < π$ ) mein hai ya lower half mein.

Step 1. Radius formula se $cos E$ invert karo: $r = a (1 - e cos E) \Rightarrow cos E = \frac{1 - r / a}{e} = \frac{1 - 1.2}{0.4} = \frac{- 0.2}{0.4} = - 0.5$ . Yeh step kyun? $r$ sirf $cos E$ par depend karta hai, aur $cos$ ek even function hai: $cos E = cos (2 π - E)$ . Isliye $arccos (- 0.5)$ ke ek orbit mein do solutions hain, $E = 2 π /3$ (upper half) aur $E = 4 π /3$ (lower half). Radius akela inhe alag nahi kar sakta.

Step 2. "Receding" use karke branch resolve karo. Yeh raha rule, clearly: $0 < E < π$ (upper half) ke liye planet perihelion se aphelion ki taraf move kar raha hai, isliye $r$ increase ho raha hai — woh receding hai. $π < E < 2 π$ (lower half) ke liye woh aphelion se wapas perihelion ki taraf move kar raha hai, isliye $r$ decrease ho raha hai — woh approaching hai. Problem receding kehta hai, isliye hum upper-half root lete hain: $E = 2 π /3 = 2.09440$ rad. Yeh step kyun? Yahi exam trap ka poora point hai: $cos E = - 0.5$ ambiguous hai, aur extra physical word ("receding") arccos branch select karta hai. Approaching hota toh $E = 4 π /3$ force karta.

Step 3. $M = E - e sin E = 2.09440 - 0.4 sin (2 π /3) = 2.09440 - 0.4 (0.86603) = 2.09440 - 0.34641 = 1.74799$ rad. Yeh step kyun? Clock ki taraf forward map.

Step 4. $t - t_{p} = \frac{M}{2 π} T = \frac{1.74799}{6.28319} \times 100 = 27.819$ days. Yeh step kyun? Final clock-to-time conversion.

Verify: Radius check karo: $a (1 - 0.4 cos (2 π /3)) = a (1 - 0.4 (- 0.5)) = a (1 + 0.2) = 1.2 a$ ✓. Agar hum (galti se) approaching branch $E = 4 π /3$ lete toh $M = 4.18879 - 0.4 (- 0.86603) = 4.53520$ aur $t = 72.18$ days milta — jo "receding" se contradict karta. Branch choice matters. ✓

Ex 9 — Cell J: degrees-vs-radians aur sign trap

Worked example Ek student ko

$M = 3 0^{\circ}$ , $e = 0.2$ diya gaya hai aur woh blindly $E = M + e sin E$ likhta hai calculator ke degree mode mein, $E \approx 30 + 0.2 sin (3 0^{\circ}) = 30. 1^{\circ}$ pata karta hai. Kya galat hai, aur sahi $E$ kya hai? Forecast: Kepler's equation ek area ( $\frac{1}{2} a^{2} E$ ) se derive hui thi, jahan angle radians mein hona chahiye. Guess karo degree-mode jawab kitna kharab hai.

Step 1. Unit error pakdo. $M = E - e sin E$ mein, $sin$ ke bahar wala $E$ ek arc/area angle hai aur radians mein hona chahiye; $sin E$ term bhi apna argument radians mein chahti hai. " $30$ " (degree number) ko $sin (3 0^{\circ}) = 0.5$ ke saath mix karna silently do alag angle systems use karta hai. Yeh step kyun? Parent note ki derivation ne sector area $\frac{1}{2} a^{2} E$ use ki — woh formula sirf radians ke saath sahi hai. Degrees se woh identity hi toot jaati hai.

Step 2. Pehle input convert karo: $M = 3 0^{\circ} = 30 \times \frac{π}{180} = 0.52360$ rad. Yeh step kyun? Kepler's equation touch karne se pehle hamesha radians mein jao, taaki har term ek hi language bolti rahe.

Step 3. $M = E - 0.2 sin E$ ko $M = 0.52360$ ke saath Newton se solve karo, seed $E_{0} = 0.52360$ : $f (0.52360) = 0.52360 - 0.2 sin (0.52360) - 0.52360 = - 0.2 (0.5) = - 0.1$ ; $f^{'} (0.52360) = 1 - 0.2 cos (0.52360) = 1 - 0.2 (0.86603) = 0.82679$ . $E_{1} = 0.52360 + 0.1/0.82679 = 0.64453$ . Ek aur step: $f (0.64453) = 0.64453 - 0.2 sin (0.64453) - 0.52360 = - 0.00090$ ; $f^{'} = 0.84058$ ; $E_{2} = 0.64453 - (- 0.00090) /0.84058 = 0.64560$ , converged. Yeh step kyun? Yeh sahi transcendental solve hai; note karo $E > M$ ? Nahi — yahan $sin E > 0$ isliye $M < E$ , matlab $E = 0.6456 > M = 0.5236$ , upper-half sign rule. Aur $- e sin E$ ek subtraction hai, isliye $E = M + e sin E$ likhna ek fixed-point rearrangement hai, "correction ek baar add karo" nahi.

Step 4. Compare karo: sahi $E = 0.64560$ rad $= 36.9 9^{\circ}$ . Student ka degree-mode " $30. 1^{\circ}$ " lagbhag $7^{\circ}$ off hai. Yeh step kyun? Damage quantify karta hai: unit error chhota nahi — usne poora correction kho diya.

Verify: Radians mein wapas plug karo: $0.64560 - 0.2 sin (0.64560) = 0.64560 - 0.2 (0.60181) = 0.64560 - 0.12036 = 0.52524 \approx 0.52360$ ; refine karne par $E = 0.64494$ rad milta hai residual ke saath 5 dp tak $0$ . ✓ Aur $0.64494$ rad $= 36.9 5^{\circ}$ , confirming karta hai ki student ka $30. 1^{\circ}$ unit confusion se galat tha. ✓

Quick recall

Recall Kaun si direction ko iteration chahiye?

$E \to M$ (aur isliye $E \to t$ ) direct substitution hai; $M \to E$ transcendental hai aur Newton–Raphson chahiye. ::: Forward free hai, inverse iterative hai.

Upper half ( $0 < E < π$ ) mein correction $- e sin E$ ka sign? ::: Negative, isliye $M < E$ .

Lower half ( $π < E < 2 π$ ) mein sign? ::: $sin E < 0$ , isliye correction positive hai aur $M > E$ .

Given only $r$ , $E$ uniquely kyun nahi milta? ::: $r = a (1 - e cos E)$ $cos E$ (even function) par depend karta hai, do branches $E$ aur $2 π - E$ deta hai; choose karne ke liye extra info chahiye (approaching vs receding).

Kepler's equation mein kaun si units use karni chahiye, aur kyun? ::: Radians — yeh sector area $\frac{1}{2} a^{2} E$ se derive hui thi, jo sirf radians mein valid hai.

Dekho bhi: Numerical root-finding — Newton–Raphson, Time of flight and orbit propagation, True anomaly ν and the orbit equation r = a(1−e²)/(1+e cos ν), Eccentricity and ellipse geometry (a, b, ae, b²=a²(1−e²)), aur $e > 1$ case ke liye Hyperbolic Kepler equation M = e·sinh F − F.