The Direction Problem: On the Causal Origin of Oriented Motion in Complex Systems
Author: Fabio Ghioni Ph.D. Affiliation: Ordinative Sciences Foundation Research Labs — ordinativescience.foundation Date: May 2026 (v1.0) · September 2026 (v1.1) Version: Preprint v1.1 Framework: github.com/anckhalion/ordinative_sciences_framework Companion to: The Collapse Equation (DOI: 10.5281/zenodo.19932312)
Abstract
While constructing a predictive model for phase transitions in complex systems, we encountered a dilemma: the empirically measured acceleration constant (g_j) increases as the system approaches the transition point. In a conventional causal framework — where the past drives the present — a system consuming its internal resources should decelerate, not accelerate. The measured 67% increase in g_j is incompatible with a push-only dynamics.
Starting from a simple observation — that every functional system exhibits a vector oriented toward the future — we follow a chain of logical eliminations. The orientation cannot originate in the past (which does not contain the novelty toward which the system moves). It cannot be random (randomness has no direction). The only remaining option is that the source of the orientation is in the destination itself.
We formalize this as the Principle of Causal Inversion: every determined result occupies a coordinate in a space where time is treated as a navigable dimension, and from that coordinate it emits a signal — a pull — that intensifies with proximity. Observable events are not caused by the past pushing the present forward; they are the local expressions of a future attractor pulling the system toward it.
The principle resolves five open problems: the origin of genuine novelty, the arrow of time, the source of spontaneous impulses in neuroscience, the anomalous pre-transition acceleration in complex systems, and the multi-scale coherence of simultaneous events without lateral coordination. We present empirical evidence from a reaction-diffusion model of civilizational dynamics, where two predicted events were confirmed in timing and structural type. The principle is falsifiable: we specify the conditions under which it would fail.
Terminological note: We use coherent and decoherent not in the quantum-mechanical sense but structurally: coherent denotes the domain of determined potential (structured but not yet manifested); decoherent denotes the domain of observable expression. The parallel with QM usage is intentional — both describe the transition from superposed potential to definite state — but the principle operates at all scales, not only the quantum scale. We use terminal to denote any bounded subsystem (individual, institution, organism, nation) that receives a signal and responds according to its internal state. Collapse denotes the moment when potential becomes expression, not destruction.
Keywords: causal inversion, attractor dynamics, direction problem, teleological causation, pre-transition acceleration, complex systems, arrow of time, principle of least action
Status: v1.1 — notation aligned with the Ordinative Sciences Symbol Canon v1.2; under observation and continuous refinement.
1. The Dilemma
1.1 Context
In the course of constructing a mathematical model for phase transitions in civilizational systems [1], we derived an empirical constant, g_j, that measures the rate at which a complex system accelerates toward its next transition point. The model treats the system as traversing a phase space following a quadratic trajectory:
IC(t) = v₀·t + ½·g_j·t²
The preliminary calibration yielded g_j ≈ 0.045 IC/month². Two subsequent events — predicted by the model in both timing and structural type — provided direct calibration:
- μ₁: 7 April 2026 (predicted: ~7 April). Structural type confirmed.
- μ₂: 29 April 2026, effective 1 May (predicted: ~1 May). Structural type confirmed.
Recalibration yielded g_j = 0.075 IC/month² — 67% higher than the preliminary estimate.
1.2 The Problem
In the original model, g_j was interpreted as a consumption rate: the system degrades an internal capacity as it traverses the phase. Under this reading — the push model — the system accelerates because it is losing internal resources.
But a system that is losing internal resources should not accelerate. It should decelerate. A machine running out of fuel slows down; it does not speed up. The 67% increase in g_j is incompatible with a push-only dynamics.
The energy for the acceleration must come from somewhere. Not from the past (the system has less capacity than before). Not from the present (nothing in the current state explains the increase). The question becomes:
Where does the additional acceleration come from?
2. The Observation
2.1 Every Functional System Has a Direction
Every functional system — every entity that does something — exhibits a vector oriented toward the future: a cell divides toward the next cell, an organism grows toward maturity, a river flows toward the sea, a relationship moves toward stabilization or dissolution, an enterprise moves toward its next cycle, a civilization moves toward its next form.
In each case, the movement has a direction — it points toward something that does not yet exist in the present.
2.2 The Property Is Scale-Invariant
This orientation appears at every scale of observation: subatomic, molecular, cellular, organismic, social, civilizational, cosmological. If every part of a system has this property, then the system as a whole has it. And the system is itself part of a larger system, which has it at its own scale. The property is fractal.
2.3 Direction Implies Destination
This is a logical necessity. A direction is a vector with orientation. Orientation means: pointing toward something. If there is nothing to point toward, there is no orientation. If there is no orientation, there is no direction. If there is no direction, there is no systematic motion — only diffusion.
But the systems we observe are not diffusing. They are converging. The acceleration constant g_j — which is not decreasing — measures the rate of this convergence.
Therefore: the destination must exist.
3. The Elimination
If the direction exists and the destination exists, where is the destination, and what is its causal relationship to the system?
Before examining the options formally, consider the argument in its simplest form. A directed vector requires a point of arrival. If the push comes from the past, then the past must contain the point of arrival. But if the past contains it, the point of arrival is already given — it is not in the future but in the past, and the system is not moving toward something new but unfolding something old. If the point of arrival is genuinely new, the past does not contain it. Therefore it is in the future. But what moves the system toward a point in the future? Not the past (which does not contain it). Not the present (which is where the system is, not where it is going). The movement must be generated by the future itself — by the destination.
We now examine this reasoning formally through three options.
3.1 Option 1: The Direction Comes from the Past
This is the standard causal framework. The system's state at time t is determined by its state at t−δt plus laws.
Problem — novelty. If the direction is fully determined by the past, then the destination is already contained in the initial conditions. But the destination is typically new. A caterpillar's body does not contain a butterfly. If the past contains the destination, it is not new. If the past does not contain it, it cannot cause it. "Emergence" names this phenomenon but does not explain it.
Problem — the energy of acceleration. In the push model, g_j should remain constant or decrease. The measured increase is incompatible.
3.2 Option 2: The Direction Is Random
Problem — coherence. Random motion is diffusion. Diffusion produces dispersion, not convergence. g_j is not a diffusion coefficient; it is an acceleration constant.
Problem — multi-scale synchrony. If direction were random at each scale independently, simultaneous convergent events at multiple scales with the same structural signature would be negligibly improbable. Yet this is precisely what is observed [1].
3.3 Option 3: The Direction Comes from the Destination
The destination exists at a future coordinate and exerts a pull on the system. This option has no logical contradiction. It explains the increasing g_j (pull intensifies with proximity). It explains novelty (the new form is not in the past; it is in the destination). It explains multi-scale synchrony (the destination pulls all scales simultaneously). It explains the direction itself.
The only objection is intuitive: it requires the destination to "exist" before it is "reached." This objection assumes that time is a one-directional flow in which the future does not exist until it becomes the present. We examine this assumption next.
4. The Principle
4.1 Time as a Navigable Dimension
In classical mechanics, space and time are treated differently: space is navigable, time is not. In relativity, this distinction weakens: space and time merge into spacetime. In quantum mechanics, it weakens further: the Schrödinger equation is time-symmetric; the arrow of time is not a property of fundamental laws but of boundary conditions.
We propose taking this weakening to its logical conclusion: time, like space, is a dimension in which coordinates can be occupied by real structures. Just as a point in space can contain an object that exerts gravitational pull, a point in time can contain a determined configuration — an attractor — that exerts a pull on systems at earlier coordinates.
4.2 The Principle of Least Action — A Precedent
The most foundational principle in classical and quantum mechanics is the Principle of Least Action: a system evolves along the trajectory that extremizes the action integral over the entire path from initial to final state:
S = ∫[t_i → t_f] L(q, dq/dt, t) dt
This principle is structurally teleological: it requires knowledge of the final state t_f to determine the trajectory at every intermediate point. Physicists have used it for over 250 years without calling it teleological, because the formalism can be reformulated locally (Euler-Lagrange equations). But the global formulation — the more fundamental one — is explicitly future-referencing.
What we propose is making explicit what the Principle of Least Action already implies: the endpoint participates causally in the trajectory.
4.3 Formal Statement
Definition 1 — Ordinative Space 𝕋: (from Latin ordinare, to arrange, to give order — the attractor imposes structural order on the temporal dimension; the term derives from a broader research programme, the Ordinative Sciences, whose full framework is available at the repository listed on the title page, but this paper is self-contained) A space in which time is treated as a navigable dimension alongside spatial dimensions. Every point has coordinates (x, τ). A point in 𝕋 can be occupied by a real structure.
Definition 2 — Terminal: A terminal is any bounded subsystem within a larger system — an individual, an institution, an organism, a nation — that receives the attractor signal and responds according to its internal state. We denote terminals as T_i with internal state s_i(t).
Definition 3 — Attractor: A determined configuration 𝔸 occupying a coordinate (x_𝔸, τ_𝔸) in 𝕋, where τ_𝔸 > τ_now. The attractor emits a signal σ_𝔸 that propagates through 𝕋.
Definition 4 — Signal Intensity:
σ_𝔸(S, 𝔸) = κ / d(S, 𝔸)^α
where d is distance in 𝕋, κ is a coupling constant, and α > 0. As S approaches 𝔸, d → 0 and σ_𝔸 → ∞.
4.4 Theorem — Causal Inversion
Let S be a system in 𝕋 with non-zero direction vector d(t) ≠ 0 for t ∈ [t₀, t_now]. Then:
(i) There exists 𝔸 ∈ 𝕋 at τ_𝔸 > τ_now such that d(t) · (𝔸 − S(t)) > 0. Direction implies destination.
(ii) The acceleration is proportional to signal intensity: g_j(t) ∝ σ_𝔸(S(t), 𝔸). Therefore g_j increases as S approaches 𝔸. Acceleration increases with proximity. Empirically: g_j = 0.045 → 0.075 with decreasing distance. ✓
(iii) For subsystems {T_i} ⊂ S, each with state s_i(t), the response is r_i(t) = φ_i(σ_𝔸(t), s_i(t)). Signal σ_𝔸 is shared; responses r_i differ by state. One signal, many expressions. This produces the observed multi-scale harmonic pattern. ✓
(iv) If s_i(t) is incompatible with σ_𝔸(t) — i.e., r_i(t) · d(t) ≤ 0 — then T_i is resolved (terminated) at time and in form coherent with 𝔸, independent of the terminal's age, scale, or longevity. Incompatibility is resolved by the attractor, not by the system.
4.5 The Resistance of Form
The attractor signal does not operate in a void. Every system through which the signal propagates has a form — a coherent configuration that maintains itself. The properties that make a form coherent (stability, self-reproduction, predictability) are the same properties that resist change. Coherence is resistance, viewed from the perspective of the transition.
Definition — Form Resistance ϱ: The form resistance ϱ_i of a terminal T_i is the structural capacity of T_i's current configuration to maintain itself against the attractor signal. ϱ is not a separate force added to the form — it IS the form, viewed from the perspective of the signal that would dissolve it.
This resistance manifests identically at every scale: a crystal resists melting, a cell resists apoptosis, an institution resists reform, a person resists identity change, a civilization resists phase transition. A form with zero resistance would have zero coherence and would not exist as a form.
Between successive transition events, ϱ temporarily dominates the local dynamics. This manifests as nostalgia for the previous state, attempts at restoration, declarations that "things will return to normal," defense of institutions that have lost function. This is not irrationality — it is the survival instinct of the form, operating through the terminals that compose it.
4.6 Free Fall: Vacuum and Medium
In gravitational physics, a body falling in vacuum accelerates at g indefinitely. A body falling through a medium (air, water) experiences drag that increases with velocity. At terminal velocity, drag equals gravitational pull and acceleration stops.
Vacuum (pure ordinative space): If a terminal had zero form resistance (ϱ = 0), it would accelerate toward the attractor at the rate determined solely by σ_𝔸. All terminals would fall identically.
Medium (decoherent realm): Every terminal has ϱ > 0 because every terminal has form. The form resistance acts as drag. Terminals with low ϱ (e.g., financial systems) express the signal rapidly. Terminals with high ϱ (e.g., agricultural systems) express it slowly.
Proposition — Terminal Velocity: A terminal reaches terminal velocity when ϱ(v) = σ_𝔸(t). Different terminals reach different terminal velocities because ϱ differs. The signal σ_𝔸 is the same for all.
The Δt between senza vista (Italian, literally "without sight": an event that has occurred in the coherent domain but is not yet observable) and con vista ("with sight": the same event as it becomes observable) is therefore a property of the terminal's form resistance, not of the signal. Financial terminals express in hours (low ϱ). Agricultural terminals express in months (high ϱ). The signal arrived at the same moment for both.
Critical observation: During decomposition, the form loses coherence progressively. As ϱ decreases, the medium thins. The system continues to accelerate because the air is disappearing. The apparent increase in g_j may reflect not an intensifying signal but a thinning medium.
4.7 G_j and Apparent g_j
Definition — Fundamental Ordinative Constant G_j: G_j is the fundamental constant of the ordinative field — the intrinsic pull of the attractor in 𝕋. Constant. Same everywhere, every scale, every terminal. Property of the space itself. Isomorphic to G (gravitational constant) in physics. G ≠ g; G is universal, g is local.
Definition — Apparent Acceleration g_j: g_j is the acceleration measured by a decoherent observer in clock-time. It is G_j filtered through:
(a) Temporal compression: Near the attractor, temporal bubbles compress (less clock-time per bubble). Same G_j in bubble-space appears as higher g_j in clock-space.
(b) Medium thinning: As the form decomposes, ϱ decreases and drag drops. The terminal accelerates through a rarefying medium.
Formally:
g_j(t) = G_j × (dτ_bubble / dτ_clock)⁻² × 1/(1 + ϱ(t)/σ_𝔸(t))
The value g_j = 0.075 IC/month² is g_j apparent, not G_j. The 67% increase from 0.045 may reflect temporal compression, medium thinning, or both. G_j itself may be constant.
4.8 The Ordinative Equivalence Principle
In 1907, Einstein realized that a person in free fall cannot distinguish between floating in empty space and falling in a gravitational field. All bodies fall at the same rate regardless of composition.
Proposition — Ordinative Equivalence: In the ordinative field of an attractor 𝔸, all terminals at the same distance experience the same fundamental acceleration G_j, regardless of their internal structure, scale, complexity, or nature.
Observable differences arise entirely from: (i) ϱ_i (form resistance / drag), (ii) s_i (internal state at collapse), (iii) φ_i (response function).
In vacuum (ϱ = 0): all terminals fall identically. Galileo's feather and cannonball.
In the decoherent medium (ϱ > 0): each terminal falls through its own "atmosphere," producing different apparent velocities. But G_j is the same for all.
The acceleration toward the attractor is not a property of the system. It is a property of the space. Every terminal is subject to it equally. The only variable is how much each terminal's form resists — and even that resistance is finite and decreasing as the attractor approaches.
5. Consequences: Five Open Problems Resolved
5.1 The Origin of Genuine Novelty
Problem: If every state is determined by the previous state, nothing genuinely new can appear.
Resolution: Novelty comes from the attractor, not from the past. The past provides the material; the attractor provides the direction and the form.
5.2 The Arrow of Time
Problem: Fundamental equations are time-symmetric. Why does the universe have a preferred direction?
Resolution: The arrow of time is the pull of the attractor. Systems move toward the future because the attractor is at τ_𝔸 > τ_now. The arrow emerges from the structure; it does not need to be postulated separately.
5.3 The Source of Spontaneous Impulses
Problem: Neuroscience describes the mechanism of spontaneous impulses but not their origin. Why this impulse? Why now? Why toward this object? For reactive impulses, the cause is in the environment. For spontaneous ones — creative desire, unexplained attraction, intuition that precedes data — the cause is not in the environment.
Resolution: A spontaneous impulse is the attractor signal received by an individual terminal. The neural circuitry is the channel, not the source. The form depends on the person's state at reception — different people express the same signal differently, as different instruments play the same score with different timbres.
5.4 Anomalous Pre-Transition Acceleration
Problem: Systems approaching phase transitions accelerate. The push explanation is circular (unstable because approaching transition, approaching because unstable). The energy is unaccounted for.
Resolution: The energy comes from the attractor, not from the system's past. Structurally identical to gravitational free fall: a body falling toward Earth accelerates because the attractor pulls, not because the body consumes internal energy. g_j = 0.045 → 0.075 is the direct measurement.
5.5 Multi-Scale Coherence Without Coordination
Problem: Structurally analogous events occur simultaneously at different scales without lateral communication.
Resolution: The events are harmonic: each scale receives the same attractor signal and responds with its own expression. Like piano strings resonating from a tuning fork — not because they communicate, but because the source is the same.
6. Empirical Evidence
6.1 The Collapse Equation
Two predicted events confirmed in both timing and structural type:
| Event | Predicted | Observed | Δt | Type confirmed |
|---|---|---|---|---|
| μ₁ | ~7 Apr 2026 | 7 Apr 2026 | 0 days | Yes |
| μ₂ | ~1 May 2026 | 29 Apr/1 May 2026 | 0 days | Yes |
Recalibrated g_j = 0.075 (67% increase). Retroactive t₀ = 5–6 February 2026, coinciding with original prediction.
6.2 Multi-Scale Harmonic Signature
Within 48 hours of μ₂, structurally analogous events (fragmentation of relational fields) manifested at four scales without lateral coordination: geopolitical (OPEC exit), financial (non-dollar petroleum settlement), national (regional refusal of federal protocols), social (non-unionized logistics strikes).
7. Falsifiability
F1. g_j must increase with proximity. If constant or decreasing → falsified.
F2. Multi-scale events must share structural signatures. If unrelated types → falsified.
F3. Spontaneous impulses must have retrospective coherence. If consistently incoherent → falsified.
F4. Six remaining micro-junctions predicted (May–July 2026). If neither timing nor type confirmed for two consecutive → fundamental revision required.
F5. No unfalsifiable auxiliary hypotheses permitted. We commit to publishing disconfirmations.
8. Relation to Existing Physics
8.1 The Principle of Least Action
The Lagrangian/Hamiltonian formulation is structurally teleological — it requires the endpoint to compute the trajectory. We make the causal implication explicit.
8.2 Retrocausality in Quantum Mechanics
Cramer's Transactional Interpretation (1986), Aharonov's Two-State Vector Formalism (1964), Wheeler's Delayed-Choice experiments (1978, realized by Jacques et al. 2007) — all include the future participating causally in the present. These are mainstream physics, not fringe proposals.
8.3 Attractors in Dynamical Systems
Standard nonlinear dynamics uses attractors as future states that pull trajectories. We extend the concept from abstract phase space to 𝕋 where time is a navigable dimension.
8.4 Teleology in Biology
Deacon's Incomplete Nature (2012) and Kauffman's Investigations (2000) both argue that biological systems exhibit goal-directed behavior irreducible to mechanical causation.
9. Discussion
9.1 What This Principle Is Not
Not predestination. The attractor determines direction, not path. The path is free.
Not mysticism. Formalized, falsifiable, with specific conditions for failure.
Not denial of push causation. Past conditions still matter. Push is a valid first-order approximation far from the attractor. Pull becomes dominant near the attractor — precisely where push fails and anomalous acceleration appears.
9.2 The Gravitational Isomorphism
| Property | Gravity | Attractor pull |
|---|---|---|
| Fundamental constant | G (universal) | G_j (universal) |
| Local acceleration | g (depends on mass, distance) | g_j apparent (depends on bubbles, ϱ) |
| Source | Mass at spatial coordinate | Determined configuration at temporal coordinate |
| Signal | Gravitational field | Attractor signal σ_𝔸 |
| Medium | Atmosphere (drag) | Form coherence ϱ (drag) |
| Terminal velocity | v_term when drag = mg | v_term when ϱ = σ_𝔸 |
| Equivalence | All bodies fall at same g | All terminals fall at same G_j |
| Vacuum | Same trajectory for all | Same trajectory for all |
The isomorphism suggests gravity and attractor pull may be the same structural principle in different dimensions of 𝕋.
9.3 Open Questions
Q1. The value of G_j. The measured g_j = 0.075 is apparent, filtered through temporal compression and form resistance. What is the underlying G_j?
Q2. The decay exponent α (gravity: α = 2; attractor: to be determined).
Q3. Universality of G_j (truly universal like G, or attractor-mass-dependent?).
Q4. Mathematical relationship between harmonics at different scales.
Q5. The temporal bubble compression factor: how does dτ_bubble/dτ_clock change as function of distance from attractor?
Q6. Empirical measurement of ϱ: candidate proxies include rate of institutional restoration attempts, frequency of "return to normal" declarations, criminalization rate of early adopters.
10. Conclusion
We began with a dilemma: an acceleration constant that should not increase was increasing. We followed the logic: every functional system has a direction; direction implies destination; the destination is not in the past and is not random. The only remaining option is that the destination exists at a future coordinate and pulls.
The formalization yields a layered structure. At the deepest level: G_j, the fundamental ordinative constant — the intrinsic pull of the attractor in 𝕋, the same for every terminal at every scale. At the observable level: g_j apparent — G_j filtered through the compression of temporal bubbles near the attractor and the drag of form resistance ϱ. The Ordinative Equivalence Principle guarantees that all terminals fall with the same G_j; the differences we observe arise from the medium (form resistance) and the measurement frame (clock-time vs bubble-time), not from the field itself.
The form resistance ϱ — the survival instinct of every coherent structure — is not a force opposed to the attractor. It is the form itself, viewed from the perspective of the transition. As the system decomposes and the form loses coherence, ϱ decreases, the medium thins, and the apparent acceleration increases. What we observe as "everything accelerating" may be not the signal getting louder but the resistance getting quieter.
This is not a radical departure. The Principle of Least Action has been structurally teleological for 250 years. Retrocausal quantum mechanics is mainstream. The Equivalence Principle has been the foundation of general relativity for over a century. What we have done is extend these structures from spatial dimensions into the temporal dimension of 𝕋.
What people experience as desire — the pull toward something that does not yet exist — may be the individual-scale expression of the same field that G_j governs at every scale. The neural circuitry is the channel. The form resistance is the drag. The source is the attractor. The note is different for each person, but the fundamental frequency is the same — and the field is one.
This work was developed within the Ordinative Sciences research programme. It is offered as a formal tool for anyone — in any discipline, at any scale of inquiry — who encounters the direction problem and finds the push-only framework insufficient.
Version Note — v1.1 (September 2026)
This revision aligns the paper's notation with the Ordinative Sciences Symbol Canon v1.2 (ratified 19 August 2026). Four symbols used in v1.0 collided with glyphs reserved elsewhere in the programme register and have been replaced throughout:
| v1.0 | v1.1 | Reason |
|---|---|---|
| 𝒜 (attractor) | 𝔸 | 𝒜 is reserved programme-wide for the Author |
| 𝒯 (ordinative space) | 𝕋 | 𝒯 is reserved for the Expressive Terminal |
| σ (signal intensity) | σ_𝔸 | bare σ is reserved for the singularity (OST) |
| ρ (form resistance) | ϱ | ρ is reserved for resonance (TE/PA/SA) |
No definition, claim, equation structure, or empirical statement has been altered: the change is notational. The companion reference has been updated to The Collapse Equation v1.3. Empirical status note: the pre-registered predictions of §7 (F4) have since been evaluated — the 2026 phase closed with eight micro-junctions, six confirmed in both timing and structural type, and the macro-junction confirmed as an extended window centred on 10–11 July 2026 (see the Collapse Equation v1.3, §12). The v1.0 text of §6–§7 is otherwise preserved as the pre-registration record.
References
- F. Ghioni, "The Collapse Equation: Predicting Phase Transitions," Ordinative Sciences Foundation, July 2026, v1.3. DOI: 10.5281/zenodo.19932312
- F. Ghioni, "Ordinative Set Theory: A Concise Operational Guide for AI, v2.1," Ordinative Sciences Foundation, March 2026. DOI: 10.5281/zenodo.18944713
- J. G. Cramer, "The transactional interpretation of quantum mechanics," Rev. Mod. Phys., 58(3), 647–687, 1986.
- Y. Aharonov, P. G. Bergmann, J. L. Lebowitz, "Time symmetry in the quantum process of measurement," Phys. Rev., 134, B1410–B1416, 1964.
- J. A. Wheeler, "The 'past' and the 'delayed-choice' double-slit experiment," in Mathematical Foundations of Quantum Theory, 1978.
- V. Jacques et al., "Experimental realization of Wheeler's delayed-choice gedanken experiment," Science, 315(5814), 966–968, 2007.
- T. W. Deacon, Incomplete Nature: How Mind Emerged from Matter, W. W. Norton, 2012.
- S. A. Kauffman, Investigations, Oxford University Press, 2000.
- M. J. Feigenbaum, "Quantitative universality for a class of nonlinear transformations," J. Stat. Phys., 19(1), 25–52, 1978.
- S. H. Strogatz, Nonlinear Dynamics and Chaos, 2nd ed., Westview Press, 2015.
- H. Goldstein, C. Poole, J. Safko, Classical Mechanics, 3rd ed., Addison-Wesley, 2002.