mathkaveli

In our previous posts, we’ve explored elliptical motion, energy, 3D orbits, rotating frames, and even complex numbers. Now we arrive at a thrilling conclusion: what happens when we disturb the system just a little? Today we explore the edge of predictability — welcome to the world of chaos. --- What Is Chaos in Physics? Chaos doesn’t mean randomness. In mathematics and physics, chaos means: The system is deterministic (it follows fixed rules) But it is **extremely sensitive** to initial conditions This is sometimes called the **butterfly effect**: a tiny change in one part of a system leads to large differences later. --- From Ellipses to Chaos Let’s start with elliptical motion again: r(t) = Now imagine that something disturbs the motion — a nearby object’s gravity, a small push, or even a tiny variation in mass. We’ll simulate this by introducing a **perturbation**: a small oscillation added to one of the axes. x(t) = a * cos(ωt) y(t) = b * sin(ωt) + ε * sin(μ...

So far in this series, we’ve described elliptical motion using position vectors, trigonometry, and rotation matrices.  But what if I told you there’s a more elegant way to represent all of it — using complex numbers?  Today, we’ll explore how complex numbers and Euler’s formula allow us to model elliptical motion, circular motion, and rotation with surprising simplicity. Along the way, we’ll translate our real-valued vector math into the complex plane — and simulate it all in Octave.  Complex Numbers as Vectors  A complex number is written as: z = x + i·y, where: x is the real part, y is the imaginary part, and i is the square root of -1.  You can think of this as a 2D vector: z = (x, y) Euler’s Formula: exp(iθ) = cos(θ) + i·sin(θ) This means we can represent circular motion using the exponential expression: z(t) = R · e^(i·ω·t) This traces a circle of radius R, rotating counterclockwise at angular speed ω. What About Ellipses?  To model an ellipse, we scal...

In our last post, we lifted ellipses into 3D space to explore planetary orbits. But here’s a twist — what if the ellipse itself spins while the particle moves along it? In other words, what happens if the entire orbital plane rotates over time? This is not just a fun thought experiment — it reflects real-world effects, like the precession of orbits, gyroscopes, and the Coriolis effect on Earth. Today, we’ll rotate our coordinate system and simulate the result in Octave.  Combining Two Motions Let’s start with a basic elliptical path: r(t) = < a*cos(ωt), b*sin(ωt) > Now imagine the entire coordinate system spins slowly — like placing the ellipse on a turntable. We apply a rotation angle φ(t) to the whole position vector. So the new position becomes: x_rot(t) = x(t)*cos(φ(t)) - y(t)*sin(φ(t)) y_rot(t) = x(t)*sin(φ(t)) + y(t)*cos(φ(t)) This is a classic 2D rotation matrix applied to a moving object. If φ(t) increases over time, we get compound motion: the object moves ...

We've seen how a particle can move in an elliptical path and how energy varies throughout the motion. But here's something exciting — planets also move in ellipses. This is where mathematics meets astronomy. In this post, we’ll explore Kepler’s Laws, simulate a simple planetary orbit in 3D using Octave, and discover how a flat ellipse turns into a space-filling orbit. --- Kepler’s First Law: The Ellipse of Orbits In the early 1600s, Johannes Kepler studied the planets and found something surprising: Kepler’s First Law: The orbit of a planet is an ellipse, with the sun at one focus. This means that planetary motion is not circular — but elliptical, just like the math we’ve been studying! --- Position in 3D: Tilting the Ellipse Let’s say the 2D elliptical orbit is: r(t) = < a*cos(θ), b*sin(θ) > To lift this into 3D, we rotate the ellipse out of the xy-plane. That’s done by applying a tilt angle (called inclination) and optional rotation around the z-axis. W...

In our previous posts, we explored how a particle traces out an ellipse and what happens when its speed changes. But motion isn’t just about where and how fast something goes — it’s also about energy. Today, we’ll explore the kinetic and potential energy of a particle moving along an elliptical path. We'll use physics to understand the motion and simulate the energies using Octave. --- The Position and Velocity Revisited Let’s return to our position vector: r(t) = < a*cos(ωt), b*sin(ωt) > From this, the velocity is: v(t) = < -a*ω*sin(ωt), b*ω*cos(ωt) > This gives us the magnitude of velocity (speed): |v(t)| = ω \* sqrt( a²*sin²(ωt) + b²*cos²(ωt) ) --- Kinetic Energy (KE) Kinetic energy is given by the formula: KE = (1/2) \* m \* |v(t)|² So for our particle of mass m , we get: KE(t) = (1/2) \* m \* ω² \* ( a²*sin²(ωt) + b²*cos²(ωt) ) This energy varies over time as the particle speeds up and slows down along the ellipse. --- Potential Energy (PE) ...

In our last post, we explored how a particle moves along an elliptical path, using a parametric position vector: r(t) = \ That model assumed a constant angular speed, ω (omega). But real life isn’t always so smooth. What happens when ω is not constant? In this post, we’ll explore how changing speed affects the motion of a particle. We’ll go from elegant ellipses to beautiful spirals, and we’ll simulate them using **Octave**, a free programming tool for doing math and science with code. --- Motion with Changing Angular Speed In our original setup, we had: x(t) = a*cos(ωt) y(t) = b*sin(ωt) If ω is constant, the motion traces out an ellipse over and over again. But what if we let the angular speed change over time — for example: ω(t) = ω₀ + k\*t This means that as time increases, the angle grows **faster** (if k > 0) or **slower** (if k r(t) = \ where θ(t) = ∫₀ᵗ ω(s) ds = ω₀\*t + (1/2)*k*t² Now we’re dealing with **non-uniform angular motion**, and the particle no ...

Have you ever watched a satellite orbit a planet or marveled at the way planets orbit the sun? Their paths often trace out more than just perfect circles — they move in ellipses. But how do we describe such motion mathematically? And how can we measure how far an object has moved if it's not traveling in a straight line? Let’s take a look at a mathematical model for this type of motion, break it down, and even compute how far the particle travels using some Octave code (a free alternative to MATLAB). --- The Position Vector: A Mathematical Story In physics, the position of a moving particle in 2D can be described with a position vector. Here's a special one that traces out an ellipse: r(t) = < a*cos(ωt), b*sin(ωt) > This vector has: a : the horizontal radius of the ellipse b : the vertical radius ω (omega): the angular speed of the particle t : time --- What Shape Is This? Let’s prove it's an ellipse. We define: x(t) = a*cos(ωt) y(t...