The Torus Becomes the Canvas A Jacobi Theta Function lives naturally on a complex torus. Instead of drawing it on a flat plane, we wrapped the mathematics onto the torus. The surface is driven by θ₁(z|τ), with its moving zeros, phase winding, and logarithmic derivative shaping the colour, seams, and raised divisor points across the geometry. What you are see is a periodic quantum-like field painted onto the space where it actually belongs. #JacobiTheta #ComplexTorus #MathematicalArt #ComplexAnalysis #RiemannSurfaces #MathAnimation

Understanding functions and their inverses visually 📈 If f(x) maps x-->y, then the inverse maps y-->x. Their graphs are reflections across the line y=x. Mathematics becomes clearer when you see it. #Math #Mathematics #InverseFunction #Functions #MathVisualization #MathAnimation #Maths #STEM #Engineering #Calculus #Algebra #MathTeacher #MathStudent #MathEducation #MathCommunity #LearnMath #MathLovers #MathIsBeautiful #DataScience #MachineLearning #AI #STEMEducation #ViralMath

vector field visualized in motion. Watch equations turn into patterns. 🌌 #Math #VectorFields #Calculus #MathVisualization #MathArt #Maths #STEM #Science #Physics #Engineering #DynamicalSystems #Chaos #MathTwitter #MathCommunity #Visualization #MathVideo #EducationalContent #STEMLearning #MathLovers #ScienceVideo #MathAnimation #MathReels

Visualizing the impossible ✨ Rotation arms with different speeds creating stunning mathematical symmetry Math Motion = Pure Beauty ❤️ #Mathematics #MathArt #Visualization #Geometry #MathAnimation #STEM #Science #CreativeCoding #Fractals #MathLovers

Visualizing Non-Linear Transformations in action See how coordinate spaces deform, rotate, and map into new shapes through exponential maps and complex transformations. Mathematics becomes art when geometry moves! #Mathematics #NonLinearTransformation #MathVisualization #Geometry #LinearAlgebra #ComplexAnalysis #MathArt #STEM #MathAnimation #CoordinateTransformation #ExponentialMap #MathIsBeautiful #EngineeringMath #Science #Education #MathLovers #VisualMath #LearnMath #MathVideo #MathConcepts #HigherMathematics #Maths #STEMEducation #DataVisualization

Which Curve Wins the Race Under Gravity? Calculus of Variations is what happens when you stop optimizing values and start optimizing geometries. The unknown isn’t a single x…it’s the whole curve y(x). And in higher dimensions, it’s the whole surface u(x,y). And the thing you’re minimizing usually isn’t a formula you can eyeball…it’s an integral that judges the entire shape. For our first lecture, we look at the famous Brachistochrone problem. Fix two points, switch on gravity, and ask a question that sounds almost too simple…which track gets a bead from A to B in the least time? Your intuition will betray you on the first try. It’s not the straight line. It’s not the drop hard then cruise sketch either. The winner is a cycloid…the curve traced by a point on a rolling circle. In the animation, the track is the moving character. We start with an imperfect curve, run gradient descent in curve-space, and watch the geometry reshape frame by frame as T[y] collapses until it locks into the brachistochrone. Here is the breakdown of the math. Set coordinates so downhill is visually obvious without negatives. Put the start above the x-axis and the finish on it: A = (0, H), B = (L, 0), with H > 0. A track is a graph y = y(x) on x ∈ [0, L] with y(0) = H and y(L) = 0. Physics from rest: By energy conservation, height drop (H − y(x)) turns into speed, v(x) = √(2g(H − y(x))). Geometry gives the arc-length element, ds = √(1 (y′(x))²) dx. So the time element is just distance-over-speed: dt = ds / v = √(1 (y′)²) / √(2g(H − y)) · dx. That makes the travel time a single object: T[y] = ∫₀ᴸ √((1 (y′(x))²) / (2g(H − y(x)))) dx. That’s the target. Not shortest path, not steepest drop. Minimize T[y]. Then we do an honest discrete-to-variational move: parameterize curves that hit the endpoints automatically, evaluate T[y] by quadrature on a fine grid, and run gradient descent in the coefficients to shrink T. #CalculusOfVariations #Brachistochrone #Cycloid #ClassicalMechanics #Optimization #MathAnimation

What if a coffee cup and a donut are secretly the same shape? In this visual journey through Topology, we explore how shapes can transform through continuous deformation stretching and bending ,without tearing or gluing. 🔵 Circle ≈ Square ☕ Coffee Cup ≈ 🍩 Torus 🌀 Möbius Strip → One-sided surface 🔁 Genus & Euler Characteristic 🚫 Why tearing changes everything From homeomorphisms to genus classification, this video breaks down deep mathematical ideas into clean, intuitive visuals. Topology isn’t about measurements — it’s about structure. It’s the mathematics of transformation. If you’ve ever wondered how geometry becomes philosophy… this is it. #Mathematics #Topology #Homeomorphism #Genus #EulerCharacteristic #MobiusStrip #Torus #MathVisualization #PureMath #STEM #MathAnimation #LearnMath #MathIsBeautiful #AbstractMath #Geometry

How Can a Ball Bounce Infinitely Many Times in Finite Time T? 🤔 I still can’t fully wrap my head around this one. Physics trains you to believe something simple: if something happens infinitely many times, it must take infinite time. That feels obvious. A bouncing ball quietly breaks that instinct. Drop a ball from height h₀ with coefficient of restitution e, where 0 < e < 1. After each impact the speed is scaled by e, so the peak height scales like e². The bounce heights satisfy hₙ = h₀ e²ⁿ Time between bounces scales like √h. So the flight time for bounce n is Δtₙ = 2 √(2hₙ/g) = 2 √(2h₀/g) · eⁿ Now add them up: T_total = Σₙ₌₀^∞ Δtₙ = 2 √(2h₀/g) · Σₙ₌₀^∞ eⁿ = 2 √(2h₀/g) · 1/(1−e) Finite. Infinitely many impacts, but the impact times accumulate toward a limit instead of stretching out forever. Infinity shows up as an accumulation point in time, not as an endless duration. In the animation you feel it. The peaks shrink, but the rhythm tightens too. Each impact comes sooner than the last. Pulse rings on the floor go from clearly separated to a rapid blur. On a floating timeline, the impact spikes pack closer and closer until they crowd against a final marker labeled T_total. Infinite ticks. Finite interval. It feels like the universe is cheating. #ClassicalMechanics #CoefficientOfRestitution #GeometricSeries #ZenoParadox #Dynamics #MathAnimation

The Curve That Forced Mathematicians To Invent A New Calculus There is a strange curve called Brownian path. It traces a continuous path without ever jumping, but it never settles into a well-defined direction at any point, so it’s continuous everywhere but differentiable nowhere. We’re used to the idea that if something is continuous, then if you zoom in far enough it should start to look like a straight line. A Brownian path refuses. Zoom in and you don’t reveal a hidden tangent. You reveal fresh roughness. This animation is a microscope on that fact. No matter how much we magnify, the path won’t straighten out. It keeps wrinkling at every scale. #BrownianMotion #StochasticProcesses #RandomWalks #ProbabilityTheory #Fractals #MathAnimation

If Arthur Mattuck’s first-order ODE lecture taught us to see dy/dx = f(x,y) as a slope field, with solution curves that stay tangent to it, what’s the matching picture for a second-order equation? Start with a general second-order ODE, with t as the independent variable: y″ = F(t, y, y′) Now do the one substitution that unlocks the geometry: Let v = y′. One second-order equation becomes a first-order system: y′ = v v′ = F(t, y, v) So the “slope field in the (x,y) plane” idea upgrades into a vector field in the (y,v) plane, the phase plane. At each point (y,v), the system tells you the velocity (y′, v′). A solution isn’t just a graph y(t) anymore. It’s a trajectory (y(t), v(t)) that flows through the phase plane and stays tangent to that vector field at every point. If the equation is autonomous, meaning there’s no explicit t: y″ = F(y, y′), then the phase portrait in (y,v) is the natural home of the dynamics. If F does depend on t, you can still make the system autonomous by turning time into a state: y′ = v v′ = F(t, y, v) t′ = 1 Now the solution lives as a curve in 3D (t, y, v), tangent to a 3D vector field. The animation uses two examples so the idea doesn’t get fuzzy. We render trajectories as 3D curves (t, y(t), v(t)) moving forward in time. On a fixed “phase wall” at the left, you see the phase-plane vector field in (y,v), plus each trajectory projected onto that wall. A short moving tangent arrow is glued to the head of the curve so the rule stays visible: at every instant the trajectory points along (1, y′, v′). Scene A: y″ δ y′ − y y³ = γ cos(ωt) (driven Duffing, double-well) Orbits get captured by one well, then kicked across the barrier. A once-per-period Poincaré slice leaves a clean signature of the long-run motion on the wall. Scene B: y″ 2β y′ ω² y = 0 (damped oscillator) Loops collapse into spirals. The vector field tilts inward and everything drains toward the equilibrium at (0,0). #DifferentialEquations #SecondOrderODEs #ArthurMattuck #PhasePortraits #DynamicalSystems #MathAnimation

Have you ever wondered what lives between a function and its derivative? Fourier space gives a sharp answer Since dᵃ/dxᵃ turns into multiplication by (jω)ᵃ, the order a can slide continuously, so x² can morph through a whole family of in-between shapes on its way to 2x and then 2. The same dial exists in optics...a lens maps a rect aperture to a sinc pattern (a Fourier transform), and the fractional Fourier transform is the continuum of light fields between the lens plane (order 0) and the focal plane (order 1). #FractionalCalculus #FractionalDerivatives #FourierTransform #FractionalFourierTransform #SignalProcessing #Optics #MathAnimation

If Arthur Mattuck’s first-order ODE lecture taught us to see dy/dx = f(x,y) as a slope field, with solution curves that move tangent to the field, what’s the matching geometric picture when the equation is second-order?🤔 Start with a general 2nd-order ODE (think time t if you like): y″ = F(t, y, y′) Now let us do the one move that changes everything Let v = y′. Then we’ve turned one second-order equation into a first-order system y′ = v v′ = F(t, y, v) So, instead of a slope field in the (x,y) plane, we now have a vector field in the (y,v) plane (phase plane). At each point (y,v) we are told the velocity vector (y′, v′). A solution is a curve (y(t), v(t)) that stays tangent to that vector field everywhere. If the equation is autonomous (no explicit t), y″ = F(y, y′), then the phase portrait in (y,v) is the natural home of the dynamics. If it does depend on t, you can still make it autonomous by adding time as a new state y′ = v v′ = F(t, y, v) t′ = 1 Now solutions are curves in 3D (t, y, v) tangent to a 3D vector field. This animation uses two examples to keep the message sharp: We render the trajectories as 3D curves (t, y(t), v(t)) moving forward in time. On a fixed phase wall at the left, you also see the phase-plane vector field in (y,v), plus each trajectory’s projection onto that wall. A short moving tangent arrow is glued to the head of a curve so the geometric rule is visible: at every instant, the trajectory points in the direction (1, y′, v′). Scene A: y″ δ y′ − y y³ = γ cos(ωt) (Duffing, driven double-well) Orbits get captured by one well, kicked across the barrier, and the once-per-period Poincaré slice leaves a crisp signature of the long-run motion on the wall. Scene B: y″ 2β y′ ω² y = 0 (damped oscillator) Loops collapse into spirals: the vector field tilts inward and everything drains toward the equilibrium at (0,0). #DifferentialEquations #SecondOrderODEs #ArthurMattuck #PhasePortraits #DynamicalSystems #MathAnimation

We just watched Professor Arthur Mattuck kick off MIT’s ODE course with the one interpretation that most differential equations classes somehow postpone: An ordinary differential equation isn’t primarily a method hunt. It’s a geometric rule. You write dy/dx = f(x,y), and that right-hand side is literally telling you the slope your solution curve must have at each point (x,y). So I built this animation as a visual companion to that first lecture. It draws the direction field (little line elements whose slope is f(x,y)) and then shows integral curves sliding through it...curves that are tangent to the field everywhere they go. Two quick examples from the animation: For dy/dx = −x/y, the slope field steers you onto circles x² y² = R². You also see a subtle point that gets missed when everything is taught as y(x)...even when the curve exists smoothly, the graph y(x) may only exist on a limited interval (|x|<R for the upper semicircle). For dy/dx = 1 x − y, the isoclines (curves where the slope is constant) make the global behavior obvious...trajectories get funneled into a corridor and become asymptotic to the special solution y=x. You learn qualitative behavior without solving it the traditional way. #DifferentialEquations #ODEs #MITOCW #VectorFields #MathAnimation #Mathematics

When I first took ordinary differential equations, the pattern was always the same. Week 1 turns into a speedrun of methods: separation of variables, integrating factors, variation of parameters, Bernoulli, exact equations… and by Week 2 or 3 the course has quietly degenerated into hammer-picking. Spot the form, apply the recipe, move on. Mechanical! Fuuuuck!😫😫😫😫 The problem is what you don’t walk away with. You leave with a toolkit, but without a feel for what a differential equation even is, especially geometrically. And that’s a big deal, because in real modeling the equations you meet are rarely nice enough to reward memorized recipes. So you end up trained to solve toy forms, while the actual subject...the behavior, the flow, the shape of solutions stays blurry. This is why I’m biased toward the old-timers. Their old-school way of doing things always surprises me:...they’ll spend time on one idea until it sticks, instead of sprinting through a syllabus checklist. One lecture from them and you start noticing a contrast. A lot of modern teaching feels like "finish the content,". You get marched through techniques, but you’re not left with a single thought that keeps bothering you later...the kind of thought that actually pushes you toward research-level curiosity. MIT OpenCourseWare’s Professor Arthur Mattuck did that to me in his very first ODE lecture. One lecture, and your whole relationship with dy/dx = f(x,y) changes. In this segment, Prof. Mattuck is basically saying: A first-order ODE is a slope field, and a solution is a curve that moves everywhere tangent to that field. The math breakdown Write the ODE as dy/dx = f(x,y). At each point (x,y) you attach a tiny line segment with slope f(x,y). A function y = y₁(x) is a solution exactly when its graph follows those slopes:. At every x, the slope of the curve equals the slope prescribed by the field at the point on the curve. That’s the single line that unifies both viewpoints: y₁′(x) = f(x, y₁(x)). So solving the ODE and drawing an integral curve are the same statement in two languages!👌🏻 Once you see that, you can stop obsessing over whether you can write y(x) in closed form. You can start asking the questions that matter: where do solutions flow, where do they get trapped, where do they blow up, and where does existence/uniqueness fail just because the field isn’t even defined? That’s the perspective shift I wish every ODE course forces early and it’s exactly why I keep pairing math with animation. #DifferentialEquations #ODEs #VectorFields #MathAnimation #Mathematics

Spherical Harmonics 🤝 Fourier Transform I’m modeling the cell surface as a time-varying field u(θ,φ). Budding is energy shifting into higher angular frequencies. Decompose u into spherical harmonics (the Fourier basis on a sphere), evolve the coefficients aₗₘ(t) where growth/smoothing is clean, then inverse-transform back into a breathing, budding 3D surface. #SphericalHarmonics #FourierTransform #HarmonicAnalysis #SpectralMethods #SignalProcessing #AppliedMath #MathAnimation #ScientificVisualization #Biophysics

Lecture 4 on Calculus of Variations You might wonder...If I’m optimizing a shape...a curve, a surface, a whole path, what does "take the derivative and set it to zero" even mean? Do I take the damn derivative with respect to a curve/surface? 🤔 In normal calculus the variable is a number x, so the reflex is clean...f′(x)=0. In calculus of variations the variable is a whole function...the geometry itself, like a curve y(x) (or a surface z(x,y)). So the derivative can’t be a single slope. It has to be a pointwise sensitivity, i.e. how the objective reacts to tiny local deformations. You’re holding a whole shape, like a curve y(x). Your objective isn’t f(x) anymore, it’s a functional J[y], and as we've seen with our first there examples, usually an integral that depends on the entire curve (often through y and y’). To talk about a “derivative”, you do the only thing that makes sense: you nudge the entire curve by a tiny amount and see how J changes. Pick a wiggle shape η(x). It’s not random...it’s any admissible deformation direction. Admissible just means it obeys the constraints. If the endpoints are fixed, you force η(0)=η(1)=0 so the wiggle doesn’t move the endpoints. Then scale that wiggle by a small number ε and define the perturbed curve yε(x)=y(x) εη(x). Now treat ε like the usual scalar in a Taylor expansion. As ε→0, J[y εη] expands as J[y εη] = J[y] ε · (first-order term depending linearly on η) o(ε). So the difference is J[y εη] - J[y] = ε · (linear functional of η) o(ε). For the standard integral of a Lagrangian problems, that linear functional can be written as an inner product with some function of x: J[y εη] - J[y] = ε ∫ (δJ/δy)(x) η(x) dx o(ε). That’s the definition-level meaning of δJ/δy: it’s the unique pointwise sensitivity function that makes this identity true for every admissible η. If δJ/δy is positive at some x, then choosing η negative there decreases J; if δJ/δy is negative there, pushing y upward locally decreases J. It’s literally a map along the curve saying push this way to go downhill. Now translate “set the derivative to zero.” At a minimizer y*, the first-order change must vanish for every admissible wiggle: J[y* εη] − J[y*] = o(ε) for all η. Plug in the expansion and the ε-term must be zero: ∫ (δJ/δy)(x) η(x) dx = 0 for all admissible η. Here’s the crucial logic step: the only way an integral against every test function η can be zero is if the integrand itself is zero (in the usual sense used in analysis). So you get δJ/δy = 0. For the common case J[y]=∫ L(x, y, y’) dx, you can compute δJ/δy explicitly and it becomes the Euler–Lagrange expression δJ/δy = ∂L/∂y − d/dx(∂L/∂y’). So if you name the Euler–Lagrange residual as “left-hand side” R(x) = ∂L/∂y − d/dx(∂L/∂y’), then “set the derivative to zero” is exactly R(x)=0. That’s why animation works so well. You don’t have to solve R=0 in one shot. You can evolve the curve in an artificial time τ by pushing it in the downhill direction: ∂y/∂τ = −R(y). Where the residual is large, the curve moves a lot; as the residual drains toward zero, the motion dies out and the curve settles into an extremal. In our animations, we start from an intentionally ugly curve/surface. Frame by frame the functional drops, the residual drains away, and the geometry relaxes into an extremal. #CalculusOfVariations #EulerLagrange #FunctionalDerivative #GradientFlow #Optimization #MathAnimation

In Lecture 3 of our Calculus of Variations series, we hand the steering wheel to light. We don’t tell it how to bend. We only define what time costs. We render a 2D refractive-index field n(x,y). A fan of rays launches from a source and curves through the gradient, shown as volumetric glow. Then we pin a start point and an end point and do the honest thing...run gradient descent in path-space on optical time, watching the straight path deform into the Fermat minimizer. T (the optical travel time of a path) drops frame by frame, and n(y) sin θ along the converged ray stays almost flat. Snell’s law isn’t an extra rule. It’s the constant of motion you can see. 👌🏾 See the math breakdown below #CalculusOfVariations #FermatsPrinciple #SnellsLaw #GeometricOptics #Optimization #MathAnimation

Lecture 1 was our first taste of what Optimizing a Geometry really means. We didn’t tune a number, we let an entire curve y(x) reshape until the travel-time functional collapsed into the brachistochrone. Lecture 2 is the same move, just one dimension up. Now the unknown isn’t a curve...it’s a whole surface in 3D. Fix two rigid rings in space, dip them in soap, and the film that forms is nature solving a variational problem: Minimize surface area subject to those boundary circles. Your intuition will betray you again. The obvious bridge is a cylinder, but a cylinder wastes area. When you let the surface relax, it tightens, a neck forms, and the geometry settles toward the catenoid...a shape that looks engineered, but it’s just what minimizing the area functional forces. In the animation we start from a fat surface and run gradient descent in surface-space. You watch the area of the surface drop frame by frame until the film locks into its minimal shape under the boundary constraint. Pls see the comment below for the math breakdown. #CalculusOfVariations #MinimalSurfaces #PlateauProblem #Catenoid #SoapFilm #MathAnimation

Here is the breakdown of the math. Set coordinates so downhill is visually obvious without negatives. Put the start above the x-axis and the finish on it: A = (0, H), B = (L, 0), with H > 0. A track is a graph y = y(x) on x ∈ [0, L] with y(0) = H and y(L) = 0. Physics from rest: By energy conservation, height drop (H − y(x)) turns into speed, v(x) = √(2g(H − y(x))). Geometry gives the arc-length element, ds = √(1 (y′(x))²) dx. So the time element is just distance-over-speed: dt = ds / v = √(1 (y′)²) / √(2g(H − y)) · dx. That makes the travel time a single object: T[y] = ∫₀ᴸ √((1 (y′(x))²) / (2g(H − y(x)))) dx. That’s the target. Not shortest path, not steepest drop. Minimize T[y]. Then we do an honest discrete-to-variational move: parameterize curves that hit the endpoints automatically, evaluate T[y] by quadrature on a fine grid, and run gradient descent in the coefficients to shrink T. The cycloid isn’t magic...it’s where the time functional keeps dragging you. #CalculusOfVariations #Brachistochrone #Cycloid #ClassicalMechanics #Optimization #MathAnimation

manim by @manim_community Explore animated math videos & animations programmatically! Manim is an animation engine for creators. Learn more about its features & usage: github.com/manimcommunity/ma… #MathAnimation #Python #Education #OpenSource