diff --git a/_sources/intro.ipynb b/_sources/intro.ipynb
index 2a5332a..de7f47a 100644
--- a/_sources/intro.ipynb
+++ b/_sources/intro.ipynb
@@ -45,7 +45,7 @@
     "  + \\int_{S_0} \\rho \\mathbf{\\big[p \\times (\\boldsymbol{\\omega} \\times p)\\big](v_r \\cdot {n})} \\, dS\n",
     "```\n",
     "\n",
-    "The remainder of this paper concerns the torque-free motions of axisymmetric variable mass systems so $\\mathbf{M^* = 0}$. Note that though this non-whirling flow assumption is not realistic, it has been shown that the magnitude of the transverse angular rates are unaffected by such a flow assumption[2]. As stability information in this paper is derived from the magnitude of the angular rates, the non-whirling flow assumption is sufficient.\n",
+    "The remainder of this paper concerns the torque-free motions of axisymmetric variable mass systems so $\\mathbf{M^* = 0}$. Note that though this non-whirling flow assumption is not realistic, it has been shown that the magnitude of the transverse angular rates are unaffected by such a flow assumption {cite}`ekewang3`. As stability information in this paper is derived from the magnitude of the angular rates, the non-whirling flow assumption is sufficient.\n",
     "\n",
     "Referring back to Figure {numref}`fig-general-variable-mass-system`, $\\mathbf{b}_1$, $\\mathbf{b}_2$, $\\mathbf{b}_3$ are a dextral set of unit vectors affixed in B and parallel to its principal directions. The instantaneous central inertia dyadic of this axisymmetric system is:\n",
     "\n",
diff --git a/intro.html b/intro.html
index 9592037..9e1da93 100644
--- a/intro.html
+++ b/intro.html
@@ -490,7 +490,7 @@ <h1>Equations of Motion and the Angular Velocity of an Axisymmetric Systems<a cl
   + \boldsymbol{\omega} \times (\mathbf{I^*} \cdot \boldsymbol{\omega}) 
   + \frac{^B d \mathbf{I}^*}{dt} \cdot \boldsymbol{\omega} \\
   + \int_{S_0} \rho \mathbf{\big[p \times (\boldsymbol{\omega} \times p)\big](v_r \cdot {n})} \, dS\end{split}\]</div>
-<p>The remainder of this paper concerns the torque-free motions of axisymmetric variable mass systems so <span class="math notranslate nohighlight">\(\mathbf{M^* = 0}\)</span>. Note that though this non-whirling flow assumption is not realistic, it has been shown that the magnitude of the transverse angular rates are unaffected by such a flow assumption[2]. As stability information in this paper is derived from the magnitude of the angular rates, the non-whirling flow assumption is sufficient.</p>
+<p>The remainder of this paper concerns the torque-free motions of axisymmetric variable mass systems so <span class="math notranslate nohighlight">\(\mathbf{M^* = 0}\)</span>. Note that though this non-whirling flow assumption is not realistic, it has been shown that the magnitude of the transverse angular rates are unaffected by such a flow assumption <span id="id2">[<a class="reference internal" href="refs.html#id13" title="S.-M. Wang and F.O. Eke. Rotational dynamics of axisymmetric variable mass systems. J. Appl. Mech., 62(4):970-974, 1995.">Wang and Eke, 1995</a>]</span>. As stability information in this paper is derived from the magnitude of the angular rates, the non-whirling flow assumption is sufficient.</p>
 <p>Referring back to Figure <a class="reference internal" href="#fig-general-variable-mass-system"><span class="std std-numref">Fig. 1</span></a>, <span class="math notranslate nohighlight">\(\mathbf{b}_1\)</span>, <span class="math notranslate nohighlight">\(\mathbf{b}_2\)</span>, <span class="math notranslate nohighlight">\(\mathbf{b}_3\)</span> are a dextral set of unit vectors affixed in B and parallel to its principal directions. The instantaneous central inertia dyadic of this axisymmetric system is:</p>
 <div class="math notranslate nohighlight" id="equation-eq3">
 <span class="eqno">(3)<a class="headerlink" href="#equation-eq3" title="Permalink to this equation">#</a></span>\[\mathbf{I}^* \triangleq I \, \mathbf{b}_1 \mathbf{b}_1 + I \, \mathbf{b}_2 \mathbf{b}_2 + J \, \mathbf{b}_3 \mathbf{b}_3\]</div>
@@ -539,7 +539,7 @@ <h2>Spin rate<a class="headerlink" href="#spin-rate" title="Permalink to this he
 <p>Examining Equation <a class="reference internal" href="#equation-eq17">(17)</a>, the spin rate may decay, stay constant, grow, or fluctuate. Another observation that can be added to this is that the spin rate retains the polarity of its initial condition; if the initial condition is positive, then the spin rate is always non-negative, and vice versa. Further analysis in this chapter assumes a positive value for the initial spin rate. Additionally, if <span class="math notranslate nohighlight">\(k_3^2\)</span> is assumed to be a constant in Equation <a class="reference internal" href="#equation-eq17">(17)</a>, then</p>
 <div class="math notranslate nohighlight" id="equation-eq18">
 <span class="eqno">(18)<a class="headerlink" href="#equation-eq18" title="Permalink to this equation">#</a></span>\[\omega_3 = \omega_{30} \bigg( \frac{m(0)}{m}\bigg) ^{1 - \frac{R^2}{2k_3^2}}.\]</div>
-<p>Equation <a class="reference internal" href="#equation-eq18">(18)</a> asserts that the spin rate cannot fluctuate for a system with constant axial radius of gyration; it can only grow or decay monotonically, or stay constant. It can then also be inferred that a time-varying axial radius of gyration can lead to fluctuations in the spin rate. Thus, the radius of gyration has a crucial effect on the spin rate. These comments on the spin rate are in agreement with the work of Snyder and Warner <span id="id2">[<a class="reference internal" href="refs.html#id32" title="V.W. Snyder and G.G. Warner. A re-evaluation of jet damping. J. Spacecr. Rockets, 5:364-366, 1968.">Snyder and Warner, 1968</a>]</span>, and Wang and Eke <span id="id3">[<a class="reference internal" href="refs.html#id13" title="S.-M. Wang and F.O. Eke. Rotational dynamics of axisymmetric variable mass systems. J. Appl. Mech., 62(4):970-974, 1995.">Wang and Eke, 1995</a>]</span>.</p>
+<p>Equation <a class="reference internal" href="#equation-eq18">(18)</a> asserts that the spin rate cannot fluctuate for a system with constant axial radius of gyration; it can only grow or decay monotonically, or stay constant. It can then also be inferred that a time-varying axial radius of gyration can lead to fluctuations in the spin rate. Thus, the radius of gyration has a crucial effect on the spin rate. These comments on the spin rate are in agreement with the work of Snyder and Warner <span id="id3">[<a class="reference internal" href="refs.html#id32" title="V.W. Snyder and G.G. Warner. A re-evaluation of jet damping. J. Spacecr. Rockets, 5:364-366, 1968.">Snyder and Warner, 1968</a>]</span>, and Wang and Eke <span id="id4">[<a class="reference internal" href="refs.html#id13" title="S.-M. Wang and F.O. Eke. Rotational dynamics of axisymmetric variable mass systems. J. Appl. Mech., 62(4):970-974, 1995.">Wang and Eke, 1995</a>]</span>.</p>
 </section>
 <section id="transverse-rate">
 <h2>Transverse rate<a class="headerlink" href="#transverse-rate" title="Permalink to this heading">#</a></h2>
diff --git a/searchindex.js b/searchindex.js
index 2a10071..7ffb281 100644
--- a/searchindex.js
+++ b/searchindex.js
@@ -1 +1 @@
-Search.setIndex({"docnames": ["README", "geometry-of-motion", "intro", "refs", "translation_curtis", "vms_end_burn", "vms_end_burn_normallized_quaternion", "vms_radial_burn", "vms_uniform_burn"], "filenames": ["README.md", "geometry-of-motion.ipynb", "intro.ipynb", "refs.md", "translation_curtis.ipynb", "vms_end_burn.ipynb", "vms_end_burn_normallized_quaternion.ipynb", "vms_radial_burn.ipynb", "vms_uniform_burn.ipynb"], "titles": ["variable-mass-systems-book", "Geometry of Torque-free Motion and Nutation Angle Solution", "Equations of Motion and the Angular Velocity of an Axisymmetric Systems", "References", "Constants", "Define the global physical parameters for the variable mass system", "Define global variables", "Define the global physical parameters for the variable mass system", "Uniform Burn"], "terms": {"explicit": 1, "descript": 1, "inertia": [1, 2, 8], "paramet": [1, 2], "coupl": 1, "develop": [1, 2], "regard": 1, "angular": [1, 3], "speed": [1, 2], "from": [1, 2, 4], "previou": [1, 8], "section": [1, 2], "ar": [1, 2, 8], "suffici": [1, 2], "visual": [1, 5, 6, 7, 8], "veloc": [1, 8], "vector": [1, 2, 8], "bodi": [1, 2, 3], "fix": [1, 2, 5, 6, 8], "frame": [1, 5, 6, 7, 8], "howev": [1, 2], "also": [1, 2], "possibl": 1, "its": [1, 2, 5, 6, 7, 8], "an": 1, "inerti": [1, 2], "refer": [1, 2], "which": [1, 2], "focu": 1, "thi": [1, 2, 3], "These": [1, 2], "us": [1, 2], "studi": [1, 2], "variabl": [1, 2, 3], "mass": [1, 2, 3], "cylind": [1, 2, 3, 5, 6, 7, 8], "problem": 1, "follow": [1, 2], "exploit": 1, "fact": 1, "system": [1, 3], "s": [1, 2, 3, 4, 8], "momentum": [1, 3, 8], "space": 1, "nanjangud": [1, 2, 3], "2016": [1, 3], "ek": [1, 2, 3], "2018": [1, 3], "enabl": [1, 8], "evolut": [1, 2], "The": [1, 2], "gener": [1, 2], "can": [1, 2], "written": [1, 2], "bf": [1, 2], "h": [1, 3, 4, 5, 6, 7, 8], "int_": [1, 2], "v_0": [1, 2], "rho": [1, 2], "p": [1, 2, 3], "time": [1, 2], "v": [1, 3, 4], "_r": 1, "boldsymbol": [1, 2], "omega": [1, 2], "mathrm": [1, 2], "d": [1, 2, 3], "retain": [1, 2], "assumpt": [1, 2], "steadi": [1, 2], "axisymmetr": [1, 3, 5, 6, 7, 8], "rel": [1, 2], "intern": [1, 2, 3], "flow": [1, 2, 4], "reduc": 1, "equat": [1, 3], "34": 1, "i": [1, 2, 5, 6, 7, 8], "cdot": [1, 2, 8], "With": 1, "definit": 1, "3": [1, 2, 5, 6, 7, 8], "4": [1, 2, 3, 5, 6, 7, 8], "respect": 1, "compon": [1, 2, 8], "form": [1, 2], "omega_1": [1, 2], "b": [1, 2, 5, 6, 7, 8], "_1": [1, 2], "omega_2": [1, 2], "_2": [1, 2], "j": [1, 2, 3, 5, 6, 7, 8], "omega_3": [1, 2], "_3": [1, 2], "_": [1, 2, 5], "t": [1, 2, 3, 4, 8], "31": [1, 2], "37": 1, "give": [1, 2], "omega_0": [1, 2], "gamma": [1, 2, 4], "12": [1, 2, 5, 6, 7, 8], "38": 1, "altern": 1, "rigid": [1, 2], "constant": [1, 2], "uniqu": 1, "alwai": [1, 2], "same": 1, "plane": [1, 2, 5, 6], "evid": [1, 2], "compar": 1, "33": 1, "shown": [1, 2], "fig": [1, 2, 5, 6, 7, 8], "In": [1, 3], "figur": [1, 2, 4, 5, 6, 7, 8], "theta": [1, 5, 6, 7, 8], "between": 1, "beta": 1, "comput": [1, 8], "tan": 1, "1": [1, 2, 3, 4, 5, 6, 7, 8], "bigg": [1, 2], "frac": [1, 2], "omega_": [1, 2], "k_1": [1, 2], "2": [1, 2, 3, 4, 5, 6, 7, 8], "0": [1, 2, 4, 5, 6, 7, 8], "k_3": [1, 2], "evalu": [1, 2, 3], "along": [1, 2], "permit": 1, "understand": [1, 2], "both": [1, 2], "coordin": 1, "case": [1, 2], "spin": [1, 3, 5, 7, 8], "rate": [1, 4], "scalar": [1, 2], "consequ": [1, 2], "39": 1, "40": 1, "result": [1, 2], "format": 1, "two": 1, "cone": 1, "rotat": [1, 2, 3], "about": [1, 2, 8], "anoth": [1, 2], "n": [1, 2, 4], "_h": 1, "unit": [1, 2], "parallel": [1, 2], "appropri": [1, 2], "name": 1, "typic": [1, 2], "function": [1, 2], "mai": [1, 2], "grow": [1, 2], "decai": [1, 2], "dictat": 1, "instantan": [1, 2], "valu": [1, 2, 8], "radii": 1, "gyrat": [1, 2, 5], "It": [1, 2], "exce": 1, "certain": 1, "instant": [1, 2], "while": [1, 2], "other": [1, 2, 5], "repres": [1, 2], "where": [1, 2, 5, 6, 7, 8], "occur": 1, "when": 1, "doe": 1, "trace": 1, "thu": [1, 2], "surfac": [1, 2, 4], "document": [1, 3], "thei": 1, "limit": [1, 5, 6, 7, 8], "spacecraft": [1, 3], "dynam": [1, 2, 3], "literatur": [1, 2], "kaplan": [1, 3], "1976": [1, 3], "classic": 1, "mechan": [1, 2, 3, 8], "identifi": 1, "second": [1, 4], "euler": [1, 2], "likin": [1, 3], "1973": [1, 3], "said": 1, "growth": 1, "undesir": 1, "potenti": 1, "hazard": 1, "excess": 1, "might": 1, "poor": 1, "initi": [1, 2], "condit": [1, 2], "subsequ": [1, 2], "mission": 1, "oper": 1, "loss": [1, 2, 3], "commun": 1, "due": 1, "incorrect": 1, "point": [1, 4, 5, 6, 7, 8], "antenna": 1, "base": [1, 2], "station": 1, "One": 1, "wa": [1, 8], "observ": [1, 2], "power": 1, "star": 1, "48": 1, "srm": 1, "cochran": [1, 3], "kang": [1, 3], "1991": [1, 3], "flandro": [1, 3], "et": 1, "al": 1, "1987": [1, 3], "halsmer": [1, 3], "mingori": [1, 3], "1995": [1, 2, 3], "meyer": [1, 3], "1983": [1, 3], "yam": [1, 3], "1986": [1, 3], "Or": [1, 3], "1992": [1, 3], "To": 1, "mitig": 1, "activ": [1, 3], "damper": 1, "webster": [1, 3], "1985": [1, 3], "close": [1, 2], "express": [1, 2], "avail": 1, "themselv": 1, "invalu": 1, "onli": [1, 2, 8], "stabil": [1, 2, 3], "control": [1, 3], "strategi": 1, "formul": 1, "temper": 1, "ani": [1, 2, 5, 6, 7, 8], "instabl": [1, 3], "deriv": [1, 2, 5, 6, 7, 8], "next": 1, "special": 1, "class": 1, "util": 1, "construct": 1, "choic": 1, "cylindr": 1, "ideal": [1, 5, 6, 7, 8], "sever": [1, 2], "prior": 1, "rocket": [1, 2, 3, 4], "flight": [1, 4], "snyder": [1, 2, 3], "warner": [1, 2, 3], "1968": [1, 2, 3], "2004": [1, 3], "wang": [1, 2, 3], "some": [1, 2], "exampl": [1, 8], "investig": 1, "who": 1, "have": [1, 2], "incorpor": 1, "propel": [1, 3, 4], "grain": 1, "attitud": [1, 2, 3], "vehicl": 1, "compris": 2, "consum": 2, "fluid": [2, 3], "phase": 2, "f": [2, 3, 5, 6, 7, 8], "A": [2, 3], "massless": 2, "shell": 2, "c": [2, 3, 4], "volum": [2, 3], "area": [2, 4], "s_0": 2, "rigidli": 2, "attach": 2, "assum": [2, 4, 8], "enter": 2, "exit": [2, 4, 5, 6], "through": 2, "region": 2, "dash": 2, "circl": 2, "radiu": [2, 4, 5], "r": [2, 3, 5, 6, 7, 8], "At": 2, "everi": 2, "everyth": 2, "within": 2, "consid": 2, "interest": 2, "law": 2, "newton": 2, "lagrang": [2, 3], "method": 2, "etc": 2, "invok": 2, "reynold": 2, "transport": 2, "theorem": 2, "approach": 2, "ha": [2, 8], "been": [2, 8], "discuss": 2, "2012": [2, 3], "mathbf": 2, "m": [2, 3, 4, 5, 6, 7, 8], "alpha": 2, "dt": [2, 5, 6, 7, 8], "v_r": 2, "ds": 2, "dv": 2, "most": 2, "tractabl": 2, "examin": 2, "decoupl": 2, "translat": 2, "integr": [2, 8], "hard": 2, "sinc": 2, "field": 2, "known": 2, "reason": 2, "made": 2, "gase": 2, "further": 2, "lack": 2, "whirl": 2, "forc": [2, 3], "last": 2, "three": 2, "term": [2, 5, 6, 7, 8], "disappear": 2, "make": 2, "less": 2, "cumbersom": 2, "becom": 2, "big": 2, "remaind": 2, "paper": [2, 3], "concern": 2, "torqu": 2, "free": [2, 3], "so": 2, "note": 2, "though": 2, "non": 2, "realist": 2, "magnitud": 2, "unaffect": 2, "As": 2, "inform": 2, "back": 2, "dextral": 2, "set": [2, 5, 6, 7, 8], "affix": 2, "princip": 2, "direct": 2, "central": 2, "dyadic": 2, "triangleq": 2, "moment": [2, 8], "center": [2, 8], "necessarili": 2, "locat": 2, "symmetri": 2, "axi": [2, 5, 6, 8], "acceler": 2, "dot": 2, "first": 2, "right": 2, "hand": 2, "side": 2, "constrain": 2, "uniform": [2, 5, 6, 7], "u": [2, 3, 8], "text": 2, "yield": 2, "z_e": 2, "given": 2, "rho_": 2, "pi": [2, 5, 6, 7, 8], "densiti": 2, "exhaust": [2, 4], "substitut": 2, "6": [2, 4, 5, 6, 7, 8], "7": [2, 4, 5, 6, 7, 8], "8": [2, 5, 6, 7, 8], "9": [2, 4, 5, 6, 7, 8], "11": [2, 5, 6, 7, 8], "13": 2, "commonli": 2, "character": 2, "solv": 2, "easili": 2, "obtain": [2, 8], "30": 2, "exp": [2, 5], "phi": [2, 5, 6, 7, 8], "int_0": 2, "If": 2, "vari": [2, 5, 6, 7, 8], "axial": [2, 3], "ln": 2, "14": 2, "17": 2, "stai": 2, "fluctuat": 2, "ad": [2, 5, 6, 7, 8], "polar": 2, "posit": 2, "neg": 2, "vice": 2, "versa": 2, "analysi": [2, 3], "chapter": 2, "addition": 2, "2k_3": 2, "18": [2, 3], "assert": 2, "cannot": 2, "monoton": 2, "infer": 2, "lead": 2, "crucial": 2, "effect": [2, 4], "comment": 2, "agreement": 2, "work": 2, "varieti": 2, "differ": 2, "wai": 2, "defin": 2, "complex": 2, "omega_c": 2, "allow": 2, "combin": 2, "differenti": 2, "linear": 2, "_c": 2, "20": 2, "co": [2, 4, 5], "chi": [2, 5, 8], "23": 2, "21": 2, "19": 2, "10": [2, 3, 4, 5, 6, 7, 8], "sin": [2, 4, 5], "oscillatori": 2, "natur": 2, "li": 2, "whose": 2, "sqrt": [2, 5, 6, 7, 8], "chosen": 2, "simplifi": 2, "appear": 2, "solut": [2, 4, 8], "25": 2, "26": 2, "27": 2, "clearli": 2, "orthogon": 2, "geometri": 2, "22": 2, "influenc": 2, "show": [2, 4, 5, 6, 7, 8], "rescal": 2, "either": 2, "reveal": 2, "24": 2, "depend": 2, "advanc": [2, 8], "around": 2, "coincid": 2, "nutat": [2, 3, 5], "attempt": 2, "analyt": 2, "Its": 2, "numer": 2, "involv": 2, "reformul": 2, "here": 2, "model": [2, 3], "constraint": [2, 5, 6, 7, 8], "28": 2, "variat": 2, "passag": 2, "should": 2, "under": 2, "circumst": 2, "final": [2, 8], "now": 2, "mimic": 2, "referenc": 3, "book": 3, "ck91": 3, "jr": 3, "e": [3, 5, 6, 7, 8], "y": [3, 4, 5, 6, 7, 8], "nonlinear": 3, "motion": [3, 5, 6, 7, 8], "upper": 3, "stage": 3, "adv": 3, "astronaut": 3, "sci": 3, "75": 3, "345": 3, "364": 3, "emm04": 3, "o": 3, "mao": 3, "morri": 3, "substanti": 3, "appl": 3, "mech": 3, "71": 3, "190": 3, "194": 3, "ew95": 3, "behavior": 3, "62": 3, "935": 3, "940": 3, "fvm": 3, "87": 3, "g": 3, "w": [3, 5, 6, 7, 8], "k": [3, 5, 6, 7, 8], "van": 3, "moorhem": 3, "shorthil": 3, "chen": 3, "woolsei": 3, "technic": 3, "report": 3, "tr": 3, "86": 3, "072": 3, "air": 3, "propuls": 3, "lab": 3, "edward": 3, "afb": 3, "ca": 3, "jan": 3, "hm95": 3, "l": [3, 5, 6, 7, 8], "passiv": 3, "guid": 3, "dyn": 3, "5": [3, 4, 5, 6, 7, 8], "1197": 3, "1203": 3, "kap76": 3, "modern": 3, "wilei": 3, "new": 3, "york": 3, "lik73": 3, "element": 3, "engin": 3, "mcgraw": 3, "hill": 3, "mey83": 3, "x": [3, 4], "convect": 3, "solid": 3, "motor": 3, "54": 3, "173": 3, "aa": 3, "83": 3, "368": 3, "my86": 3, "thrust": [3, 4], "aiaa": 3, "astrodynam": 3, "confer": 3, "proceed": 3, "367": 3, "375": 3, "williamsburg": 3, "va": 3, "2271": 3, "nan16": 3, "On": 3, "phd": 3, "thesi": 3, "univers": 3, "california": 3, "davi": 3, "ne12": 3, "type": 3, "rev": 3, "aerosp": 3, "eng": 3, "256": 3, "260": 3, "ne18": 3, "partial": 3, "conserv": 3, "tech": 3, "doi": 3, "1016": 3, "ast": 3, "03": 3, "003": 3, "or92": 3, "rotor": 3, "pendulum": 3, "perige": 3, "assist": 3, "modul": 3, "anomali": 3, "15": 3, "297": 3, "303": 3, "sw68": 3, "re": 3, "jet": 3, "damp": 3, "spacecr": 3, "366": 3, "we95": 3, "970": 3, "974": 3, "web85": 3, "american": 3, "institut": 3, "aeronaut": 3, "import": [4, 5, 6, 7, 8], "numpi": [4, 5, 6, 7, 8], "np": [4, 5, 6, 7, 8], "matplotlib": [4, 5, 6, 7, 8], "pyplot": [4, 5, 6, 7, 8], "plt": [4, 5, 6, 7, 8], "scipi": [4, 5, 6, 7, 8], "odeint": 4, "m_e": 4, "974e20": 4, "6742e11": 4, "r_e": 4, "6371e3": 4, "earth": 4, "isp": 4, "300": 4, "specif": 4, "impuls": 4, "g0": 4, "81": 4, "a_": 4, "nozzl": 4, "p_a": 4, "101325": 4, "atmospher": 4, "pressur": 4, "sea": 4, "level": 4, "pa": 4, "p_e": 4, "vacuum": 4, "simplic": 4, "v0": 4, "gamma0": 4, "radian": 4, "90": 4, "path": 4, "angl": 4, "vertic": 4, "launch": 4, "x0": 4, "downrang": 4, "distanc": [4, 5, 6], "h0": 4, "m_wet": 4, "500e3": 4, "wet": 4, "kg": [4, 5, 6, 7, 8], "m_dry": 4, "50e3": 4, "dry": 4, "5e6": 4, "linspac": 4, "3600": 4, "1000": [4, 5, 6, 7, 8], "def": [4, 5, 6, 7, 8], "graviti": 4, "return": [4, 5, 6, 7, 8], "rocket_dynam": 4, "mdot_": 4, "els": [4, 5, 6, 7, 8], "more": 4, "burn": 4, "t_dynam": 4, "dvdt": 4, "print": [4, 8], "dgamma_dt": 4, "dxdt": 4, "dhdt": 4, "dmdt": 4, "consumpt": 4, "y0": [4, 5, 6, 7, 8], "arg": 4, "figsiz": [4, 5, 6, 7, 8], "1e3": 4, "convert": 4, "kilomet": 4, "readabl": 4, "titl": [4, 5, 6, 7, 8], "xlabel": [4, 5, 6, 7, 8], "km": 4, "ylabel": [4, 5, 6, 7, 8], "true": [4, 5, 6, 7, 8], "z": 4, "1000e3": 4, "g_z": 4, "tonn": 4, "solve_ivp": [5, 6, 7, 8], "anim": [5, 6, 7, 8], "funcanim": [5, 6, 7, 8], "rho_ini": [5, 6, 7, 8], "z0": [5, 6], "w0": [5, 6, 7, 8], "w10": [5, 6, 7, 8], "w20": [5, 6, 7, 8], "w30": [5, 6, 7, 8], "mf0": [5, 6, 7, 8], "i10": [5, 6, 7, 8], "i30": [5, 6, 7, 8], "theta0": [5, 6, 7], "f0": [5, 6, 7], "tb": [5, 6, 7, 8], "100": [5, 6, 7], "m_dot": [5, 6, 7, 8], "z_dot": [5, 6], "q0": [5, 6, 7, 8], "i1": [5, 6, 7, 8], "i3": [5, 6, 7, 8], "w1": [5, 6, 7, 8], "w2": [5, 6, 7, 8], "w3": [5, 6, 7, 8], "md": [5, 6, 7, 8], "i1d": [5, 6, 7, 8], "i3d": [5, 6, 7, 8], "ze": [5, 6, 7], "undergo": [5, 6, 7, 8], "w1d": [5, 6, 7, 8], "w2d": [5, 6, 7, 8], "w3d": [5, 6, 7, 8], "thetad": [5, 6, 7], "fd": [5, 6, 7, 8], "zd": [5, 6, 7, 8], "omega_quat": [5, 6, 7, 8], "arrai": [5, 6, 7, 8], "quat": [5, 6, 7, 8], "quat_dot": [5, 6, 7, 8], "quat_mult": [5, 6, 7, 8], "wd": [5, 6, 7, 8], "tolist": [5, 6], "gam": 5, "zero": [5, 6, 7, 8], "transvers": [5, 7, 8], "calcul": [5, 6, 7, 8], "h12": 5, "h3": [5, 6, 7, 8], "arctan": 5, "q": [5, 6, 7, 8], "x1": [5, 6, 7, 8], "y1": [5, 6, 7, 8], "z1": [5, 6, 7, 8], "x2": [5, 6, 7, 8], "y2": [5, 6, 7, 8], "z2": [5, 6, 7, 8], "01": [5, 6, 7, 8], "t_eval": [5, 6, 7, 8], "arang": [5, 6, 7, 8], "sol": [5, 6, 7, 8], "lambda": 5, "atol": [5, 6, 7, 8], "1e": [5, 6, 7, 8], "rtol": [5, 6, 7, 8], "statu": [5, 6, 7, 8], "fail": [5, 6, 7, 8], "messag": [5, 6, 7, 8], "omega1": [5, 6, 7, 8], "omega2": [5, 6, 7, 8], "omega3": [5, 6, 7, 8], "quat_norm": [5, 6, 7, 8], "norm_violation_indic": [5, 6, 7, 8], "ab": [5, 6, 7, 8], "len": [5, 6, 7, 8], "step": [5, 6, 7, 8], "satisfi": [5, 6, 7, 8], "throughout": [5, 6, 7, 8], "rang": [5, 6, 7, 8], "quat_to_euler_zxz": [5, 6, 7, 8], "psi": [5, 6, 7, 8], "arctan2": [5, 6, 7, 8], "arcco": [5, 6, 7, 8], "subplot": [5, 6, 7, 8], "label": [5, 6, 7, 8], "rad": [5, 6, 7, 8], "vs": [5, 6, 7, 8], "legend": [5, 6, 7, 8], "tight_layout": [5, 6, 7, 8], "linewidth": [5, 6, 7, 8], "ylim": [5, 6, 7, 8], "min": [5, 6, 7, 8], "max": [5, 6, 7, 8], "h1": [5, 6, 7, 8], "h2": [5, 6, 7, 8], "e1": [5, 6, 7, 8], "e2": [5, 6, 7, 8], "e3": [5, 6, 7, 8], "quat_to_rot_matrix": [5, 6, 7, 8], "animate_t_handle_quat": [5, 6, 7, 8], "specifi": [5, 6, 7, 8], "dimens": [5, 6, 7, 8], "lag": [5, 6, 7, 8], "cm": [5, 6, 7, 8], "lbc": [5, 6, 7, 8], "lad": [5, 6, 7, 8], "store": [5, 6, 7, 8], "orient": [5, 6, 7, 8], "kei": [5, 6, 7, 8], "xa": [5, 6, 7, 8], "ya": [5, 6, 7, 8], "za": [5, 6, 7, 8], "xb": [5, 6, 7, 8], "yb": [5, 6, 7, 8], "zb": [5, 6, 7, 8], "xc": [5, 6, 7, 8], "yc": [5, 6, 7, 8], "zc": [5, 6, 7, 8], "xd": [5, 6, 7, 8], "yd": [5, 6, 7, 8], "up": [5, 6, 7, 8], "window": [5, 6, 7, 8], "ax": [5, 6, 7, 8], "add_subplot": [5, 6, 7, 8], "111": [5, 6, 7, 8], "project": [5, 6, 7, 8], "3d": [5, 6, 7, 8], "set_xlabel": [5, 6, 7, 8], "set_ylabel": [5, 6, 7, 8], "set_zlabel": [5, 6, 7, 8], "set_xlim": [5, 6, 7, 8], "set_ylim": [5, 6, 7, 8], "set_zlim": [5, 6, 7, 8], "set_titl": [5, 6, 7, 8], "draw": [5, 6, 7, 8], "bc": [5, 6, 7, 8], "updat": [5, 6, 7, 8], "current": [5, 6, 7, 8], "set_data": [5, 6, 7, 8], "set_3d_properti": [5, 6, 7, 8], "interv": [5, 6, 7, 8], "blit": [5, 6, 7, 8], "axhlin": [5, 6, 7, 8], "color": [5, 6, 7, 8], "linestyl": [5, 6, 7, 8], "99": [5, 6, 7, 8], "adjust": [5, 6, 7, 8], "better": [5, 6, 7, 8], "normal": [6, 7], "prevent": [6, 7], "drift": [6, 7], "linalg": [6, 7], "reset": [6, 7], "becaus": [6, 7], "rsq0": 7, "radialburn": 7, "rsq": 7, "rsqd": 7, "mf": 7, "violat": [7, 8], "intel": 8, "mkl": 8, "warn": 8, "support": 8, "stream": 8, "simd": 8, "extens": 8, "sse4": 8, "processor": 8, "deprec": 8, "oneapi": 8, "math": 8, "kernel": 8, "librari": 8, "2025": 8, "requir": 8, "avx": 8, "instruct": 8, "1020": 8, "rho_exhaust": 8, "002": 8, "chi0": 8, "precess": 8, "fd0": 8, "uniburn": 8, "qw": 8, "qx": 8, "qy": 8, "qz": 8, "chid": 8, "success": 8, "check": 8, "reach": 8, "end": 8, "norm": 8, "over": 8, "handl": 8, "total": 8, "energi": 8, "matrix": 8, "usag": 8, "setup": 1}, "objects": {}, "objtypes": {}, "objnames": {}, "titleterms": {"variabl": [0, 5, 6, 7, 8], "mass": [0, 4, 5, 6, 7, 8], "system": [0, 2, 5, 7, 8], "book": 0, "geometri": 1, "torqu": 1, "free": 1, "motion": [1, 2], "nutat": 1, "angl": [1, 5, 6, 7, 8], "solut": [1, 5, 6, 7], "equat": [2, 4, 5, 6, 7, 8], "angular": [2, 5, 6, 7, 8], "veloc": [2, 4, 5, 6, 7], "an": 2, "axisymmetr": 2, "spin": 2, "rate": [2, 5, 6, 7, 8], "transvers": 2, "refer": 3, "constant": 4, "initi": [4, 5, 6, 7, 8], "condit": [4, 5, 6, 7, 8], "time": [4, 5, 6, 7, 8], "span": [4, 5, 6, 7, 8], "gravit": 4, "acceler": 4, "function": [4, 5, 6, 7, 8], "altitud": 4, "differenti": [4, 5, 6, 7, 8], "state": 4, "integr": [4, 5, 6, 7], "over": [4, 5, 6, 7], "grid": 4, "extract": [4, 5, 6, 7, 8], "result": [4, 5, 6, 7, 8], "plot": [4, 5, 6, 7, 8], "trajectori": 4, "g": 4, "vs": 4, "defin": [5, 6, 7, 8], "global": [5, 6, 7, 8], "physic": [5, 7, 8], "paramet": [5, 7, 8], "properti": [5, 6, 7, 8], "quaternion": [5, 6, 7, 8], "repres": [5, 6, 7, 8], "rotat": [5, 6, 7, 8], "combin": [5, 6, 7, 8], "endburn": [5, 6], "multipli": [5, 6, 7, 8], "two": [5, 6, 7, 8], "simul": [5, 6, 7, 8], "solv": [5, 6, 7, 8], "od": [5, 6, 7, 8], "check": [5, 6, 7], "wa": [5, 6, 7], "success": [5, 6, 7], "print": [5, 6, 7], "final": [5, 6, 7], "reach": [5, 6, 7], "end": [5, 6, 7], "burn": [5, 6, 7, 8], "comput": [5, 6], "norm": [5, 6, 7], "violat": [5, 6], "us": 5, "addit": 5, "ext": 5, "convert": [5, 6, 7, 8], "euler": [5, 6, 7, 8], "z": [5, 6, 7, 8], "x": [5, 6, 7, 8], "sequenc": [5, 6, 7, 8], "from": [5, 6, 7, 8], "moment": [5, 6, 7], "inertia": [5, 6, 7], "t": [5, 6, 7], "handl": [5, 6, 7], "s": [5, 6, 7], "total": [5, 6, 7], "mechan": [5, 6, 7], "energi": [5, 6, 7], "compon": [5, 6, 7], "momentum": [5, 6, 7], "about": [5, 6, 7], "center": [5, 6, 7], "matrix": [5, 6, 7], "exampl": [5, 6, 7], "usag": [5, 6, 7], "assum": [5, 6, 7], "qw": [5, 6, 7], "qx": [5, 6, 7], "qy": [5, 6, 7], "qz": [5, 6, 7], "ar": [5, 6, 7], "obtain": [5, 6, 7], "previou": [5, 6, 7], "ensur": 7, "continu": 7, "unwrap": 7, "uniform": 8}, "envversion": {"sphinx.domains.c": 2, "sphinx.domains.changeset": 1, "sphinx.domains.citation": 1, "sphinx.domains.cpp": 6, "sphinx.domains.index": 1, "sphinx.domains.javascript": 2, "sphinx.domains.math": 2, "sphinx.domains.python": 3, "sphinx.domains.rst": 2, "sphinx.domains.std": 2, "sphinx.ext.intersphinx": 1, "sphinxcontrib.bibtex": 9, "sphinx": 56}})
\ No newline at end of file
+Search.setIndex({"docnames": ["README", "geometry-of-motion", "intro", "refs", "translation_curtis", "vms_end_burn", "vms_end_burn_normallized_quaternion", "vms_radial_burn", "vms_uniform_burn"], "filenames": ["README.md", "geometry-of-motion.ipynb", "intro.ipynb", "refs.md", "translation_curtis.ipynb", "vms_end_burn.ipynb", "vms_end_burn_normallized_quaternion.ipynb", "vms_radial_burn.ipynb", "vms_uniform_burn.ipynb"], "titles": ["variable-mass-systems-book", "Geometry of Torque-free Motion and Nutation Angle Solution", "Equations of Motion and the Angular Velocity of an Axisymmetric Systems", "References", "Constants", "Define the global physical parameters for the variable mass system", "Define global variables", "Define the global physical parameters for the variable mass system", "Uniform Burn"], "terms": {"explicit": 1, "descript": 1, "inertia": [1, 2, 8], "paramet": [1, 2], "coupl": 1, "develop": [1, 2], "regard": 1, "angular": [1, 3], "speed": [1, 2], "from": [1, 2, 4], "previou": [1, 8], "section": [1, 2], "ar": [1, 2, 8], "suffici": [1, 2], "visual": [1, 5, 6, 7, 8], "veloc": [1, 8], "vector": [1, 2, 8], "bodi": [1, 2, 3], "fix": [1, 2, 5, 6, 8], "frame": [1, 5, 6, 7, 8], "howev": [1, 2], "also": [1, 2], "possibl": 1, "its": [1, 2, 5, 6, 7, 8], "an": 1, "inerti": [1, 2], "refer": [1, 2], "which": [1, 2], "focu": 1, "thi": [1, 2, 3], "These": [1, 2], "us": [1, 2], "studi": [1, 2], "variabl": [1, 2, 3], "mass": [1, 2, 3], "cylind": [1, 2, 3, 5, 6, 7, 8], "problem": 1, "follow": [1, 2], "exploit": 1, "fact": 1, "system": [1, 3], "s": [1, 2, 3, 4, 8], "momentum": [1, 3, 8], "space": 1, "nanjangud": [1, 2, 3], "2016": [1, 3], "ek": [1, 2, 3], "2018": [1, 3], "enabl": [1, 8], "evolut": [1, 2], "The": [1, 2], "gener": [1, 2], "can": [1, 2], "written": [1, 2], "bf": [1, 2], "h": [1, 3, 4, 5, 6, 7, 8], "int_": [1, 2], "v_0": [1, 2], "rho": [1, 2], "p": [1, 2, 3], "time": [1, 2], "v": [1, 3, 4], "_r": 1, "boldsymbol": [1, 2], "omega": [1, 2], "mathrm": [1, 2], "d": [1, 2, 3], "retain": [1, 2], "assumpt": [1, 2], "steadi": [1, 2], "axisymmetr": [1, 3, 5, 6, 7, 8], "rel": [1, 2], "intern": [1, 2, 3], "flow": [1, 2, 4], "reduc": 1, "equat": [1, 3], "34": 1, "i": [1, 2, 5, 6, 7, 8], "cdot": [1, 2, 8], "With": 1, "definit": 1, "3": [1, 2, 5, 6, 7, 8], "4": [1, 2, 3, 5, 6, 7, 8], "respect": 1, "compon": [1, 2, 8], "form": [1, 2], "omega_1": [1, 2], "b": [1, 2, 5, 6, 7, 8], "_1": [1, 2], "omega_2": [1, 2], "_2": [1, 2], "j": [1, 2, 3, 5, 6, 7, 8], "omega_3": [1, 2], "_3": [1, 2], "_": [1, 2, 5], "t": [1, 2, 3, 4, 8], "31": [1, 2], "37": 1, "give": [1, 2], "omega_0": [1, 2], "gamma": [1, 2, 4], "12": [1, 2, 5, 6, 7, 8], "38": 1, "altern": 1, "rigid": [1, 2], "constant": [1, 2], "uniqu": 1, "alwai": [1, 2], "same": 1, "plane": [1, 2, 5, 6], "evid": [1, 2], "compar": 1, "33": 1, "shown": [1, 2], "fig": [1, 2, 5, 6, 7, 8], "2": [1, 2, 3, 4, 5, 6, 7, 8], "setup": 1, "In": [1, 3], "figur": [1, 2, 4, 5, 6, 7, 8], "theta": [1, 5, 6, 7, 8], "between": 1, "beta": 1, "comput": [1, 8], "tan": 1, "1": [1, 2, 3, 4, 5, 6, 7, 8], "bigg": [1, 2], "frac": [1, 2], "omega_": [1, 2], "k_1": [1, 2], "0": [1, 2, 4, 5, 6, 7, 8], "k_3": [1, 2], "evalu": [1, 2, 3], "along": [1, 2], "permit": 1, "understand": [1, 2], "both": [1, 2], "coordin": 1, "case": [1, 2], "spin": [1, 3, 5, 7, 8], "rate": [1, 4], "scalar": [1, 2], "consequ": [1, 2], "39": 1, "40": 1, "result": [1, 2], "format": 1, "two": 1, "cone": 1, "rotat": [1, 2, 3], "about": [1, 2, 8], "anoth": [1, 2], "n": [1, 2, 4], "_h": 1, "unit": [1, 2], "parallel": [1, 2], "appropri": [1, 2], "name": 1, "typic": [1, 2], "function": [1, 2], "mai": [1, 2], "grow": [1, 2], "decai": [1, 2], "dictat": 1, "instantan": [1, 2], "valu": [1, 2, 8], "radii": 1, "gyrat": [1, 2, 5], "It": [1, 2], "exce": 1, "certain": 1, "instant": [1, 2], "while": [1, 2], "other": [1, 2, 5], "repres": [1, 2], "where": [1, 2, 5, 6, 7, 8], "occur": 1, "when": 1, "doe": 1, "trace": 1, "thu": [1, 2], "surfac": [1, 2, 4], "document": [1, 3], "thei": 1, "limit": [1, 5, 6, 7, 8], "spacecraft": [1, 3], "dynam": [1, 2, 3], "literatur": [1, 2], "kaplan": [1, 3], "1976": [1, 3], "classic": 1, "mechan": [1, 2, 3, 8], "identifi": 1, "second": [1, 4], "euler": [1, 2], "likin": [1, 3], "1973": [1, 3], "said": 1, "growth": 1, "undesir": 1, "potenti": 1, "hazard": 1, "excess": 1, "might": 1, "poor": 1, "initi": [1, 2], "condit": [1, 2], "subsequ": [1, 2], "mission": 1, "oper": 1, "loss": [1, 2, 3], "commun": 1, "due": 1, "incorrect": 1, "point": [1, 4, 5, 6, 7, 8], "antenna": 1, "base": [1, 2], "station": 1, "One": 1, "wa": [1, 8], "observ": [1, 2], "power": 1, "star": 1, "48": 1, "srm": 1, "cochran": [1, 3], "kang": [1, 3], "1991": [1, 3], "flandro": [1, 3], "et": 1, "al": 1, "1987": [1, 3], "halsmer": [1, 3], "mingori": [1, 3], "1995": [1, 2, 3], "meyer": [1, 3], "1983": [1, 3], "yam": [1, 3], "1986": [1, 3], "Or": [1, 3], "1992": [1, 3], "To": 1, "mitig": 1, "activ": [1, 3], "damper": 1, "webster": [1, 3], "1985": [1, 3], "close": [1, 2], "express": [1, 2], "avail": 1, "themselv": 1, "invalu": 1, "onli": [1, 2, 8], "stabil": [1, 2, 3], "control": [1, 3], "strategi": 1, "formul": 1, "temper": 1, "ani": [1, 2, 5, 6, 7, 8], "instabl": [1, 3], "deriv": [1, 2, 5, 6, 7, 8], "next": 1, "special": 1, "class": 1, "util": 1, "construct": 1, "choic": 1, "cylindr": 1, "ideal": [1, 5, 6, 7, 8], "sever": [1, 2], "prior": 1, "rocket": [1, 2, 3, 4], "flight": [1, 4], "snyder": [1, 2, 3], "warner": [1, 2, 3], "1968": [1, 2, 3], "2004": [1, 3], "wang": [1, 2, 3], "some": [1, 2], "exampl": [1, 8], "investig": 1, "who": 1, "have": [1, 2], "incorpor": 1, "propel": [1, 3, 4], "grain": 1, "attitud": [1, 2, 3], "vehicl": 1, "compris": 2, "consum": 2, "fluid": [2, 3], "phase": 2, "f": [2, 3, 5, 6, 7, 8], "A": [2, 3], "massless": 2, "shell": 2, "c": [2, 3, 4], "volum": [2, 3], "area": [2, 4], "s_0": 2, "rigidli": 2, "attach": 2, "assum": [2, 4, 8], "enter": 2, "exit": [2, 4, 5, 6], "through": 2, "region": 2, "dash": 2, "circl": 2, "radiu": [2, 4, 5], "r": [2, 3, 5, 6, 7, 8], "At": 2, "everi": 2, "everyth": 2, "within": 2, "consid": 2, "interest": 2, "law": 2, "newton": 2, "lagrang": [2, 3], "method": 2, "etc": 2, "invok": 2, "reynold": 2, "transport": 2, "theorem": 2, "approach": 2, "ha": [2, 8], "been": [2, 8], "discuss": 2, "2012": [2, 3], "mathbf": 2, "m": [2, 3, 4, 5, 6, 7, 8], "alpha": 2, "dt": [2, 5, 6, 7, 8], "v_r": 2, "ds": 2, "dv": 2, "most": 2, "tractabl": 2, "examin": 2, "decoupl": 2, "translat": 2, "integr": [2, 8], "hard": 2, "sinc": 2, "field": 2, "known": 2, "reason": 2, "made": 2, "gase": 2, "further": 2, "lack": 2, "whirl": 2, "forc": [2, 3], "last": 2, "three": 2, "term": [2, 5, 6, 7, 8], "disappear": 2, "make": 2, "less": 2, "cumbersom": 2, "becom": 2, "big": 2, "remaind": 2, "paper": [2, 3], "concern": 2, "torqu": 2, "free": [2, 3], "so": 2, "note": 2, "though": 2, "non": 2, "realist": 2, "magnitud": 2, "unaffect": 2, "As": 2, "inform": 2, "back": 2, "dextral": 2, "set": [2, 5, 6, 7, 8], "affix": 2, "princip": 2, "direct": 2, "central": 2, "dyadic": 2, "triangleq": 2, "moment": [2, 8], "center": [2, 8], "necessarili": 2, "locat": 2, "symmetri": 2, "axi": [2, 5, 6, 8], "acceler": 2, "dot": 2, "first": 2, "right": 2, "hand": 2, "side": 2, "constrain": 2, "uniform": [2, 5, 6, 7], "u": [2, 3, 8], "text": 2, "yield": 2, "z_e": 2, "given": 2, "rho_": 2, "pi": [2, 5, 6, 7, 8], "densiti": 2, "exhaust": [2, 4], "substitut": 2, "6": [2, 4, 5, 6, 7, 8], "7": [2, 4, 5, 6, 7, 8], "8": [2, 5, 6, 7, 8], "9": [2, 4, 5, 6, 7, 8], "11": [2, 5, 6, 7, 8], "13": 2, "commonli": 2, "character": 2, "solv": 2, "easili": 2, "obtain": [2, 8], "30": 2, "exp": [2, 5], "phi": [2, 5, 6, 7, 8], "int_0": 2, "If": 2, "vari": [2, 5, 6, 7, 8], "axial": [2, 3], "ln": 2, "14": 2, "17": 2, "stai": 2, "fluctuat": 2, "ad": [2, 5, 6, 7, 8], "polar": 2, "posit": 2, "neg": 2, "vice": 2, "versa": 2, "analysi": [2, 3], "chapter": 2, "addition": 2, "2k_3": 2, "18": [2, 3], "assert": 2, "cannot": 2, "monoton": 2, "infer": 2, "lead": 2, "crucial": 2, "effect": [2, 4], "comment": 2, "agreement": 2, "work": 2, "varieti": 2, "differ": 2, "wai": 2, "defin": 2, "complex": 2, "omega_c": 2, "allow": 2, "combin": 2, "differenti": 2, "linear": 2, "_c": 2, "20": 2, "co": [2, 4, 5], "chi": [2, 5, 8], "23": 2, "21": 2, "19": 2, "10": [2, 3, 4, 5, 6, 7, 8], "sin": [2, 4, 5], "oscillatori": 2, "natur": 2, "li": 2, "whose": 2, "sqrt": [2, 5, 6, 7, 8], "chosen": 2, "simplifi": 2, "appear": 2, "solut": [2, 4, 8], "25": 2, "26": 2, "27": 2, "clearli": 2, "orthogon": 2, "geometri": 2, "22": 2, "influenc": 2, "show": [2, 4, 5, 6, 7, 8], "rescal": 2, "either": 2, "reveal": 2, "24": 2, "depend": 2, "advanc": [2, 8], "around": 2, "coincid": 2, "nutat": [2, 3, 5], "attempt": 2, "analyt": 2, "Its": 2, "numer": 2, "involv": 2, "reformul": 2, "here": 2, "model": [2, 3], "constraint": [2, 5, 6, 7, 8], "28": 2, "variat": 2, "passag": 2, "should": 2, "under": 2, "circumst": 2, "final": [2, 8], "now": 2, "mimic": 2, "referenc": 3, "book": 3, "ck91": 3, "jr": 3, "e": [3, 5, 6, 7, 8], "y": [3, 4, 5, 6, 7, 8], "nonlinear": 3, "motion": [3, 5, 6, 7, 8], "upper": 3, "stage": 3, "adv": 3, "astronaut": 3, "sci": 3, "75": 3, "345": 3, "364": 3, "emm04": 3, "o": 3, "mao": 3, "morri": 3, "substanti": 3, "appl": 3, "mech": 3, "71": 3, "190": 3, "194": 3, "ew95": 3, "behavior": 3, "62": 3, "935": 3, "940": 3, "fvm": 3, "87": 3, "g": 3, "w": [3, 5, 6, 7, 8], "k": [3, 5, 6, 7, 8], "van": 3, "moorhem": 3, "shorthil": 3, "chen": 3, "woolsei": 3, "technic": 3, "report": 3, "tr": 3, "86": 3, "072": 3, "air": 3, "propuls": 3, "lab": 3, "edward": 3, "afb": 3, "ca": 3, "jan": 3, "hm95": 3, "l": [3, 5, 6, 7, 8], "passiv": 3, "guid": 3, "dyn": 3, "5": [3, 4, 5, 6, 7, 8], "1197": 3, "1203": 3, "kap76": 3, "modern": 3, "wilei": 3, "new": 3, "york": 3, "lik73": 3, "element": 3, "engin": 3, "mcgraw": 3, "hill": 3, "mey83": 3, "x": [3, 4], "convect": 3, "solid": 3, "motor": 3, "54": 3, "173": 3, "aa": 3, "83": 3, "368": 3, "my86": 3, "thrust": [3, 4], "aiaa": 3, "astrodynam": 3, "confer": 3, "proceed": 3, "367": 3, "375": 3, "williamsburg": 3, "va": 3, "2271": 3, "nan16": 3, "On": 3, "phd": 3, "thesi": 3, "univers": 3, "california": 3, "davi": 3, "ne12": 3, "type": 3, "rev": 3, "aerosp": 3, "eng": 3, "256": 3, "260": 3, "ne18": 3, "partial": 3, "conserv": 3, "tech": 3, "doi": 3, "1016": 3, "ast": 3, "03": 3, "003": 3, "or92": 3, "rotor": 3, "pendulum": 3, "perige": 3, "assist": 3, "modul": 3, "anomali": 3, "15": 3, "297": 3, "303": 3, "sw68": 3, "re": 3, "jet": 3, "damp": 3, "spacecr": 3, "366": 3, "we95": 3, "970": 3, "974": 3, "web85": 3, "american": 3, "institut": 3, "aeronaut": 3, "import": [4, 5, 6, 7, 8], "numpi": [4, 5, 6, 7, 8], "np": [4, 5, 6, 7, 8], "matplotlib": [4, 5, 6, 7, 8], "pyplot": [4, 5, 6, 7, 8], "plt": [4, 5, 6, 7, 8], "scipi": [4, 5, 6, 7, 8], "odeint": 4, "m_e": 4, "974e20": 4, "6742e11": 4, "r_e": 4, "6371e3": 4, "earth": 4, "isp": 4, "300": 4, "specif": 4, "impuls": 4, "g0": 4, "81": 4, "a_": 4, "nozzl": 4, "p_a": 4, "101325": 4, "atmospher": 4, "pressur": 4, "sea": 4, "level": 4, "pa": 4, "p_e": 4, "vacuum": 4, "simplic": 4, "v0": 4, "gamma0": 4, "radian": 4, "90": 4, "path": 4, "angl": 4, "vertic": 4, "launch": 4, "x0": 4, "downrang": 4, "distanc": [4, 5, 6], "h0": 4, "m_wet": 4, "500e3": 4, "wet": 4, "kg": [4, 5, 6, 7, 8], "m_dry": 4, "50e3": 4, "dry": 4, "5e6": 4, "linspac": 4, "3600": 4, "1000": [4, 5, 6, 7, 8], "def": [4, 5, 6, 7, 8], "graviti": 4, "return": [4, 5, 6, 7, 8], "rocket_dynam": 4, "mdot_": 4, "els": [4, 5, 6, 7, 8], "more": 4, "burn": 4, "t_dynam": 4, "dvdt": 4, "print": [4, 8], "dgamma_dt": 4, "dxdt": 4, "dhdt": 4, "dmdt": 4, "consumpt": 4, "y0": [4, 5, 6, 7, 8], "arg": 4, "figsiz": [4, 5, 6, 7, 8], "1e3": 4, "convert": 4, "kilomet": 4, "readabl": 4, "titl": [4, 5, 6, 7, 8], "xlabel": [4, 5, 6, 7, 8], "km": 4, "ylabel": [4, 5, 6, 7, 8], "true": [4, 5, 6, 7, 8], "z": 4, "1000e3": 4, "g_z": 4, "tonn": 4, "solve_ivp": [5, 6, 7, 8], "anim": [5, 6, 7, 8], "funcanim": [5, 6, 7, 8], "rho_ini": [5, 6, 7, 8], "z0": [5, 6], "w0": [5, 6, 7, 8], "w10": [5, 6, 7, 8], "w20": [5, 6, 7, 8], "w30": [5, 6, 7, 8], "mf0": [5, 6, 7, 8], "i10": [5, 6, 7, 8], "i30": [5, 6, 7, 8], "theta0": [5, 6, 7], "f0": [5, 6, 7], "tb": [5, 6, 7, 8], "100": [5, 6, 7], "m_dot": [5, 6, 7, 8], "z_dot": [5, 6], "q0": [5, 6, 7, 8], "i1": [5, 6, 7, 8], "i3": [5, 6, 7, 8], "w1": [5, 6, 7, 8], "w2": [5, 6, 7, 8], "w3": [5, 6, 7, 8], "md": [5, 6, 7, 8], "i1d": [5, 6, 7, 8], "i3d": [5, 6, 7, 8], "ze": [5, 6, 7], "undergo": [5, 6, 7, 8], "w1d": [5, 6, 7, 8], "w2d": [5, 6, 7, 8], "w3d": [5, 6, 7, 8], "thetad": [5, 6, 7], "fd": [5, 6, 7, 8], "zd": [5, 6, 7, 8], "omega_quat": [5, 6, 7, 8], "arrai": [5, 6, 7, 8], "quat": [5, 6, 7, 8], "quat_dot": [5, 6, 7, 8], "quat_mult": [5, 6, 7, 8], "wd": [5, 6, 7, 8], "tolist": [5, 6], "gam": 5, "zero": [5, 6, 7, 8], "transvers": [5, 7, 8], "calcul": [5, 6, 7, 8], "h12": 5, "h3": [5, 6, 7, 8], "arctan": 5, "q": [5, 6, 7, 8], "x1": [5, 6, 7, 8], "y1": [5, 6, 7, 8], "z1": [5, 6, 7, 8], "x2": [5, 6, 7, 8], "y2": [5, 6, 7, 8], "z2": [5, 6, 7, 8], "01": [5, 6, 7, 8], "t_eval": [5, 6, 7, 8], "arang": [5, 6, 7, 8], "sol": [5, 6, 7, 8], "lambda": 5, "atol": [5, 6, 7, 8], "1e": [5, 6, 7, 8], "rtol": [5, 6, 7, 8], "statu": [5, 6, 7, 8], "fail": [5, 6, 7, 8], "messag": [5, 6, 7, 8], "omega1": [5, 6, 7, 8], "omega2": [5, 6, 7, 8], "omega3": [5, 6, 7, 8], "quat_norm": [5, 6, 7, 8], "norm_violation_indic": [5, 6, 7, 8], "ab": [5, 6, 7, 8], "len": [5, 6, 7, 8], "step": [5, 6, 7, 8], "satisfi": [5, 6, 7, 8], "throughout": [5, 6, 7, 8], "rang": [5, 6, 7, 8], "quat_to_euler_zxz": [5, 6, 7, 8], "psi": [5, 6, 7, 8], "arctan2": [5, 6, 7, 8], "arcco": [5, 6, 7, 8], "subplot": [5, 6, 7, 8], "label": [5, 6, 7, 8], "rad": [5, 6, 7, 8], "vs": [5, 6, 7, 8], "legend": [5, 6, 7, 8], "tight_layout": [5, 6, 7, 8], "linewidth": [5, 6, 7, 8], "ylim": [5, 6, 7, 8], "min": [5, 6, 7, 8], "max": [5, 6, 7, 8], "h1": [5, 6, 7, 8], "h2": [5, 6, 7, 8], "e1": [5, 6, 7, 8], "e2": [5, 6, 7, 8], "e3": [5, 6, 7, 8], "quat_to_rot_matrix": [5, 6, 7, 8], "animate_t_handle_quat": [5, 6, 7, 8], "specifi": [5, 6, 7, 8], "dimens": [5, 6, 7, 8], "lag": [5, 6, 7, 8], "cm": [5, 6, 7, 8], "lbc": [5, 6, 7, 8], "lad": [5, 6, 7, 8], "store": [5, 6, 7, 8], "orient": [5, 6, 7, 8], "kei": [5, 6, 7, 8], "xa": [5, 6, 7, 8], "ya": [5, 6, 7, 8], "za": [5, 6, 7, 8], "xb": [5, 6, 7, 8], "yb": [5, 6, 7, 8], "zb": [5, 6, 7, 8], "xc": [5, 6, 7, 8], "yc": [5, 6, 7, 8], "zc": [5, 6, 7, 8], "xd": [5, 6, 7, 8], "yd": [5, 6, 7, 8], "up": [5, 6, 7, 8], "window": [5, 6, 7, 8], "ax": [5, 6, 7, 8], "add_subplot": [5, 6, 7, 8], "111": [5, 6, 7, 8], "project": [5, 6, 7, 8], "3d": [5, 6, 7, 8], "set_xlabel": [5, 6, 7, 8], "set_ylabel": [5, 6, 7, 8], "set_zlabel": [5, 6, 7, 8], "set_xlim": [5, 6, 7, 8], "set_ylim": [5, 6, 7, 8], "set_zlim": [5, 6, 7, 8], "set_titl": [5, 6, 7, 8], "draw": [5, 6, 7, 8], "bc": [5, 6, 7, 8], "updat": [5, 6, 7, 8], "current": [5, 6, 7, 8], "set_data": [5, 6, 7, 8], "set_3d_properti": [5, 6, 7, 8], "interv": [5, 6, 7, 8], "blit": [5, 6, 7, 8], "axhlin": [5, 6, 7, 8], "color": [5, 6, 7, 8], "linestyl": [5, 6, 7, 8], "99": [5, 6, 7, 8], "adjust": [5, 6, 7, 8], "better": [5, 6, 7, 8], "normal": [6, 7], "prevent": [6, 7], "drift": [6, 7], "linalg": [6, 7], "reset": [6, 7], "becaus": [6, 7], "rsq0": 7, "radialburn": 7, "rsq": 7, "rsqd": 7, "mf": 7, "violat": [7, 8], "intel": 8, "mkl": 8, "warn": 8, "support": 8, "stream": 8, "simd": 8, "extens": 8, "sse4": 8, "processor": 8, "deprec": 8, "oneapi": 8, "math": 8, "kernel": 8, "librari": 8, "2025": 8, "requir": 8, "avx": 8, "instruct": 8, "1020": 8, "rho_exhaust": 8, "002": 8, "chi0": 8, "precess": 8, "fd0": 8, "uniburn": 8, "qw": 8, "qx": 8, "qy": 8, "qz": 8, "chid": 8, "success": 8, "check": 8, "reach": 8, "end": 8, "norm": 8, "over": 8, "handl": 8, "total": 8, "energi": 8, "matrix": 8, "usag": 8}, "objects": {}, "objtypes": {}, "objnames": {}, "titleterms": {"variabl": [0, 5, 6, 7, 8], "mass": [0, 4, 5, 6, 7, 8], "system": [0, 2, 5, 7, 8], "book": 0, "geometri": 1, "torqu": 1, "free": 1, "motion": [1, 2], "nutat": 1, "angl": [1, 5, 6, 7, 8], "solut": [1, 5, 6, 7], "equat": [2, 4, 5, 6, 7, 8], "angular": [2, 5, 6, 7, 8], "veloc": [2, 4, 5, 6, 7], "an": 2, "axisymmetr": 2, "spin": 2, "rate": [2, 5, 6, 7, 8], "transvers": 2, "refer": 3, "constant": 4, "initi": [4, 5, 6, 7, 8], "condit": [4, 5, 6, 7, 8], "time": [4, 5, 6, 7, 8], "span": [4, 5, 6, 7, 8], "gravit": 4, "acceler": 4, "function": [4, 5, 6, 7, 8], "altitud": 4, "differenti": [4, 5, 6, 7, 8], "state": 4, "integr": [4, 5, 6, 7], "over": [4, 5, 6, 7], "grid": 4, "extract": [4, 5, 6, 7, 8], "result": [4, 5, 6, 7, 8], "plot": [4, 5, 6, 7, 8], "trajectori": 4, "g": 4, "vs": 4, "defin": [5, 6, 7, 8], "global": [5, 6, 7, 8], "physic": [5, 7, 8], "paramet": [5, 7, 8], "properti": [5, 6, 7, 8], "quaternion": [5, 6, 7, 8], "repres": [5, 6, 7, 8], "rotat": [5, 6, 7, 8], "combin": [5, 6, 7, 8], "endburn": [5, 6], "multipli": [5, 6, 7, 8], "two": [5, 6, 7, 8], "simul": [5, 6, 7, 8], "solv": [5, 6, 7, 8], "od": [5, 6, 7, 8], "check": [5, 6, 7], "wa": [5, 6, 7], "success": [5, 6, 7], "print": [5, 6, 7], "final": [5, 6, 7], "reach": [5, 6, 7], "end": [5, 6, 7], "burn": [5, 6, 7, 8], "comput": [5, 6], "norm": [5, 6, 7], "violat": [5, 6], "us": 5, "addit": 5, "ext": 5, "convert": [5, 6, 7, 8], "euler": [5, 6, 7, 8], "z": [5, 6, 7, 8], "x": [5, 6, 7, 8], "sequenc": [5, 6, 7, 8], "from": [5, 6, 7, 8], "moment": [5, 6, 7], "inertia": [5, 6, 7], "t": [5, 6, 7], "handl": [5, 6, 7], "s": [5, 6, 7], "total": [5, 6, 7], "mechan": [5, 6, 7], "energi": [5, 6, 7], "compon": [5, 6, 7], "momentum": [5, 6, 7], "about": [5, 6, 7], "center": [5, 6, 7], "matrix": [5, 6, 7], "exampl": [5, 6, 7], "usag": [5, 6, 7], "assum": [5, 6, 7], "qw": [5, 6, 7], "qx": [5, 6, 7], "qy": [5, 6, 7], "qz": [5, 6, 7], "ar": [5, 6, 7], "obtain": [5, 6, 7], "previou": [5, 6, 7], "ensur": 7, "continu": 7, "unwrap": 7, "uniform": 8}, "envversion": {"sphinx.domains.c": 2, "sphinx.domains.changeset": 1, "sphinx.domains.citation": 1, "sphinx.domains.cpp": 6, "sphinx.domains.index": 1, "sphinx.domains.javascript": 2, "sphinx.domains.math": 2, "sphinx.domains.python": 3, "sphinx.domains.rst": 2, "sphinx.domains.std": 2, "sphinx.ext.intersphinx": 1, "sphinxcontrib.bibtex": 9, "sphinx": 56}})
\ No newline at end of file