Wednesday, June 17, 2026
The Tribe of Math Mentors: 11 Questions Every Educator Should Answer
Friday, June 5, 2026
Beyond Taylor Series: The Magic and History of Padé Approximations
If you have ever taken a calculus class, you probably remember the Taylor series. It is the mathematical magic trick that lets you turn complicated functions—like sines, cosines, and exponentials—into simple, infinitely long polynomials. For centuries, it has been a cornerstone of numerical mathematics.
But the Taylor series has a dark secret: it frequently breaks. If a function has a vertical asymptote (a pole) or if you move too far from your starting point, the Taylor series spirals out of control into infinity. It is strictly bounded by what mathematicians call a "radius of convergence."
Enter the Padé approximation. Instead of using a single polynomial to estimate a function, a Padé approximant uses a fraction (a ratio of two polynomials). This simple structural change unlocks a profound level of mathematical power, allowing us to see past the limits of Taylor series and model complex, chaotic systems in modern physics and engineering.
A Brief History: From France to Quantum Physics
While brilliant minds like Carl Gustav Jacob Jacobi and Ferdinand Georg Frobenius tinkered with rational approximations, the rigorous study of this method belongs to late 19th-century France. The prominent mathematician Charles Hermite utilized these rational fractions in 1873 to brilliantly prove that e (Euler's number) is a transcendental number.
Hermite passed this fascination down to his student, Henri Eugène Padé (1863–1953), at the École Normale Supérieure in Paris. In 1892, Padé published his doctoral dissertation, Sur la représentation approchée d'une fonction par des fractions rationnelles. Padé didn't just invent the approximation; he organized it. He created the "Padé table," a meticulous grid where columns represent the degree of the numerator and rows represent the degree of the denominator, proving deep theorems about how sequences of these fractions converge.
Sadly, despite his brilliant thesis, Padé spent the rest of his career mostly in administration as a university chancellor. Polynomials were simply easier to calculate by hand, so for decades, his rational fractions were ignored.
The Modern Revival
Everything changed in the mid-20th century with the invention of computers. Physicists realized that traditional power series were failing them in complex fields like quantum mechanics and fluid dynamics. They needed a way to model phase transitions, scattering resonances, and chaotic particle interactions—phenomena that mathematically look like "division by zero."
Computers made solving the heavy linear algebra required for Padé approximations instantaneous. Today, Padé approximants are essential in theoretical physics, control theory, and electrical engineering. They are actively used to map out critical temperatures in thermodynamics, design stable electronic filters, and reduce incredibly complex computer graphics algorithms into lightweight, fast-running equations.
The Underlying Mechanism: Why It Works
The brilliance of the Padé approximation is how it turns a messy, non-linear division problem into a highly solvable system of linear equations.
Suppose we want to approximate a function f(x). We know its standard Taylor series looks like this:
f(x) = c0 + c1*x + c2*x^2 + c3*x^3 + ...
Instead of stopping there, we define a rational function R(x) as the ratio of a numerator polynomial P(x) of degree m, and a denominator polynomial Q(x) of degree n. By convention, we set the first term of the denominator to 1:
R(x) = P(x) / Q(x) P(x) = a0 + a1*x + a2*x^2 + ... + am*x^m Q(x) = 1 + b1*x + b2*x^2 + ... + bn*x^n
We want R(x) to match the Taylor series of f(x) as closely as possible, up to the term x^(m+n).
The Linearization Trick
Trying to divide P(x) by Q(x) and match it to f(x) directly is an algebraic nightmare. Padé's underlying mechanism bypasses this by simply multiplying both sides of the equation by the denominator, Q(x). This is called linearization:
f(x) * Q(x) - P(x) = 0
This equation must hold true for all terms up to x^(m+n). When you expand f(x) * Q(x), you group the terms by their powers of x. Since the result must be zero, the coefficient for every power of x must equal exactly zero. This creates a clear, two-part system of equations.
- Step 1: Finding the Denominator (The 'b' coefficients): We look at the higher powers of x (from x^(m+1) up to x^(m+n)). Because the numerator P(x) stops at degree m, its 'a' coefficients completely disappear from these equations. We are left with a neat system of linear equations (specifically, a Toeplitz matrix) to solve for the 'b' coefficients.
- Step 2: Finding the Numerator (The 'a' coefficients): Once we have the 'b' values for the denominator, we look at the lower powers of x (from x^0 to x^m). Finding the 'a' coefficients is now just basic multiplication and substitution.
Why Rationality Defeats Polynomials
Why go through this trouble? Because fractions can do things that straight polynomials simply cannot:
- Modeling Singularities: A polynomial can never divide by zero. If your function shoots off to infinity (a vertical asymptote), a Taylor series fails entirely. But a Padé approximant has a denominator. When Q(x) = 0, the function naturally produces a vertical asymptote, perfectly modeling the singularity.
- Handling Infinite Horizons: Polynomials always race toward positive or negative infinity as x gets very large. If you are modeling a physical system that levels off (like terminal velocity or a capacitor charging), Taylor series are useless. A Padé approximant where the degree of P(x) is equal to or less than Q(x) has a built-in horizontal asymptote, making it stable at infinity.
- Breaking the Radius of Convergence: In complex analysis, Taylor series get blocked by singularities, limiting their "radius of convergence." Padé approximations mathematically "absorb" these singularities into their denominator. This allows the function to converge far beyond the original limits—a powerful process related to analytic continuation.
The next time you see a seemingly impossible power series or a diverging equation, remember Henri Padé. Sometimes, the solution isn't to add more terms to your polynomial; the solution is to divide it.
Here's a GeoGebra file I made with this idea in mind. Note the degree of the numerator and denominator are switched from the description above; otherwise, it should make sense.
Monday, June 1, 2026
The Möbius Assassin: Sifting Primes over Morning Coffee
When you spend enough time looking at the structural foundations of mathematics, you start to appreciate the tools that do the heavy lifting behind the scenes. Today, I want to talk about one of those tools—a beautiful proof showing how the Möbius function acts as an elegant filter for prime numbers.
If you are a multiplicative function, the Möbius function is the ultimate editor. Here is how it works, broken down piece by piece.
The Cast of Characters
Before we get to the proof, let's define the three main concepts at play.
1. Multiplicative Functions A function f(n) is "multiplicative" if it respects numbers that have no common factors. If "a" and "b" share no prime factors, then: f(a * b) = f(a) * f(b)
Every integer is just a unique fingerprint of prime numbers multiplied together. Because our function is multiplicative, we don't need to evaluate it for every massive number to infinity; we only need to understand how it behaves for prime powers (like 2^3 or 5^2), and we can build the rest from there.
2. The Möbius Function, μ(n) Think of μ(n) as a strict bouncer at the prime factorization club. It only outputs three possible numbers: 1, -1, or 0.
Rule 1: μ(1) = 1
Rule 2: If a number is made of distinct, single primes (like 6 = 2 * 3), then μ(n) = (-1)^k, where k is the number of primes. So, μ(p) = -1 for any single prime p.
Rule 3: If a number has a squared prime in its factorization (like 4 = 2^2, or 18 = 2 * 3^2), the bouncer kicks it out: μ(n) = 0.
3. Dirichlet Convolution When you see a sum that looks like this: Σ (for d|n) [ μ(d) * f(d) ] ...it just means "take every divisor 'd' of 'n', plug it into this expression, and add them all up." A massive theorem tells us that if you combine two multiplicative functions using this exact "sum over divisors" method, the resulting master function is also perfectly multiplicative.
The Setup
We want to prove that if f(1) = 1, then: Σ (for d|n) [ μ(d) * f(d) ] = (1 - f(p1)) * (1 - f(p2)) * ... * (1 - f(pk)) where n = p1^a1 * p2^a2 * ... * pk^ak is the prime-power factorization of n.
The Proof
Step 1: Define our master function Let's define a new function, g(n), to represent our big sum: g(n) = Σ (for d|n) [ μ(d) * f(d) ]
We know f(n) and μ(n) are multiplicative. Multiplying them together creates another multiplicative function. Combining them via our divisor sum preserves this property. Therefore, g(n) is a multiplicative function.
Step 2: Break it into prime building blocks Because g(n) is multiplicative, we can split our function into a product of its prime pieces based on the factorization of n: g(n) = g(p1^a1) * g(p2^a2) * ... * g(pk^ak)
Step 3: The Möbius Assassin goes to work Let's zoom in on just one of those prime blocks, p^a. The divisors of a prime power like p^a are simply 1, p, p^2, ..., p^a. Let's write out the sum for g(p^a):
g(p^a) = μ(1)f(1) + μ(p)f(p) + μ(p^2)f(p^2) + ... + μ(p^a)f(p^a)
Now, apply the rules of the Möbius bouncer. We know μ(1) = 1 and μ(p) = -1. But any term with a squared prime or higher gets a zero! So μ(p^2), μ(p^3), all the way to μ(p^a) are instantly erased.
Our long, intimidating sum collapses beautifully: g(p^a) = (1)f(1) + (-1)f(p) + 0 + 0 + ... + 0
Since we are given that f(1) = 1, this simplifies down to exactly: g(p^a) = 1 - f(p)
Step 4: Putting it all together In Step 2, we established that g(n) is just all of its prime-power pieces multiplied together. In Step 3, we proved that every single prime-power piece simplifies to [ 1 - f(p) ].
Substitute our simplified prime pieces back into the product: g(n) = (1 - f(p1)) * (1 - f(p2)) * ... * (1 - f(pk))
Since g(n) was just our placeholder for the original sum, we arrive right at our destination.
Q.E.D.
Friday, May 29, 2026
From Peasant Arithmetic to Modern Cryptography: The Hidden Magic of Halving and Doubling
As educators, one of our greatest challenges is humanizing mathematics. We often present mathematical operations as sterile, modern inventions delivered from on high, stripping away the messy, brilliant, human history behind them.
Take multiplication, for example. Most of us learned the standard algorithm—stack the numbers, multiply the digits, carry the tens, add it all up. But what if there was a completely different way to multiply? A way that requires no memorization of times tables beyond multiplying and dividing by two?
Enter the Russian Peasant Algorithm. It is a brilliant piece of historical arithmetic that seems like a parlor trick at first glance. But pull back the curtain, and you will find it is the exact same mathematical engine driving the encryption of your modern web browser.
Let's break down how this centuries-old method builds a bridge straight into abstract algebra and computer science.
The Algorithm: Halving and Doubling
The rules of Russian Peasant Multiplication are remarkably simple. Let's say we want to multiply 57 by 623.
Create two columns. Write 57 on the left and 623 on the right.
In the left column, continually halve the number (ignoring any remainders/decimals) until you reach 1.
In the right column, continually double the number to match the rows on the left.
Finally, look at the left column. If a number is EVEN, cross out that entire row. If a number is ODD, keep it.
Add up all the "kept" numbers in the right column.
Let's run the numbers:
| Left (Halve) | Right (Double) | Action | Keep |
| 57 | 623 | 57 is odd -> Keep | 623 |
| 28 | 1246 | 28 is even -> Ignore | - |
| 14 | 2492 | 14 is even -> Ignore | - |
| 7 | 4984 | 7 is odd -> Keep | 4984 |
| 3 | 9968 | 3 is odd -> Keep | 9968 |
| 1 | 19936 | 1 is odd -> Keep | 19936 |
Now, sum the kept values on the right:
623 + 4984 + 9968 + 19936 = 35511.
Check it with a calculator. 57 x 623 is exactly 35511.
Why It Works: Binary in Disguise
Why does this seemingly random process of crossing out numbers work? Because this algorithm is a mechanical, analog way of converting a number into binary (base-2) and applying the distributive property.
When you repeatedly halve the number on the left and note whether it is odd or even, you are actually finding the powers of 2 that make up that number.
Odd means there is a remainder of 1 (a "1" in binary).
Even means there is a remainder of 0 (a "0" in binary).
Looking at our left column for 57:
57 is odd -> includes 2^0 (which is 1)
28 is even -> skips 2^1 (which is 2)
14 is even -> skips 2^2 (which is 4)
7 is odd -> includes 2^3 (which is 8)
3 is odd -> includes 2^4 (which is 16)
1 is odd -> includes 2^5 (which is 32)
So, 57 is really just (1 + 8 + 16 + 32).
Because of the distributive property, we can rewrite our original problem:
57 x 623 = (1 + 8 + 16 + 32) x 623
Which expands to:
57 x 623 = (1 x 623) + (8 x 623) + (16 x 623) + (32 x 623)
Notice how the numbers 2 and 4 are completely missing. Because 57 does not contain a 2 or a 4 in its binary expansion, we crossed out the partial products for 623 x 2 and 623 x 4! The Russian peasants were doing binary multiplication centuries before the first computer was plugged in.
Taking It Further: The Engine of Abstract Algebra
Taking this algorithm out of basic arithmetic and dropping it into abstract algebra is where it transforms into a cornerstone of modern computational mathematics. We stop looking at it as multiplying two integers, and instead view the first column as an integer scalar, and the second column as an element in an arbitrary algebraic structure.
1. Group Theory: The "Double-and-Add" Algorithm
In an additive group, if we want to compute the scalar multiplication n * g, the naive approach is g + g + g... done n times.
The Russian Peasant algorithm reduces this from linear time to logarithmic time. "Halving" breaks down the scalar n. "Doubling" means computing g + g = 2g, then 2g + 2g = 4g. "Adding" means accumulating the elements where n is odd.
This exact logic is the engine behind Elliptic Curve Cryptography (ECC). In ECC, 'g' is a point on an elliptic curve. The Double-and-Add algorithm is the only computationally feasible way to multiply a curve point by a massive 256-bit cryptographic private key.
2. Rings and Fields: "Square-and-Multiply"
If we shift from an additive structure to a multiplicative one—like computing x^n in a finite field—the algorithm flips perfectly from addition to multiplication.
"Halving" remains the binary decomposition of the exponent n.
"Doubling" the second column becomes squaring: x -> x^2 -> x^4 -> x^8.
"Adding" the kept values becomes multiplying them together.
To compute x^57, the algorithm gives us:
x^57 = x^32 * x^16 * x^8 * x^1
This is essential for RSA encryption and primality testing, where computing a^(p-1) mod p must be done with gigantic primes without overflowing memory.
3. Polynomial Rings over GF(2)
Consider the polynomial ring GF(2)[x], where coefficients are strictly 0 or 1, and addition is equivalent to the logical XOR operation.
If we multiply two polynomials A(x) and B(x):
Halving A(x) is a bitwise right-shift.
Checking odd/even is checking if it has a +1 constant term.
Doubling B(x) is a bitwise left-shift.
Adding the kept shifted versions of B(x) is done using XOR.
This process is known as Carry-less Multiplication. It is natively built into modern processors (via the CLMUL instruction set) because it is the exact math required to compute AES encryption.
Conclusion
The beauty of the Russian Peasant algorithm isn't just that it works; it's that it reveals the underlying architecture of mathematics. It shows our students that math isn't a series of disconnected rules, but a continuous thread stretching from fields in ancient Eurasia straight into the silicon of their smartphones. By exploring these connections, we don't just teach calculation—we humanize the subject and invite students to see the wonder within the algorithms.
Friday, May 22, 2026
Unlocking the Matrix: A Guide to Solving Systems of 3 Equations
If you have ever felt like algebra is just moving letters around a page until a number falls out, solving systems of equations in three variables might seem like the ultimate boss battle. Instead of just dealing with x and y on a flat piece of paper, a third variable (like z, or in our case r, s, and t) brings the math into three-dimensional space. Every equation represents a flat plane, and solving the system means finding the exact 3D coordinate where all three planes intersect.
But before we dive into a full 3x3 system, let's talk about what happens when the number of equations doesn't perfectly match the number of variables.
The Goldilocks Problem: Too Many or Not Enough Equations?
For a system to have a single, beautifully neat solution, you generally need exactly as many equations as you have variables. But what if the balance is off?
Not Enough Equations (Underdetermined Systems): Let's say you have 3 variables but only 2 equations. This is like having two planes intersecting in 3D space. Instead of meeting at a single point, they meet along an entire line. Because that line stretches on forever, you have an infinite number of solutions. You simply don't have enough restrictions to pin down one specific answer.
Too Many Equations (Overdetermined Systems): What if you have 3 variables but 4 equations? This is like throwing a fourth plane into the mix. For a single solution to exist, that fourth plane must magically intersect the exact same point as the other three. Usually, this doesn't happen, resulting in no solution (an inconsistent system) unless one of the equations is just a redundant copy or combination of the others.
In Practice: Solving a 3x3 System
Let's look at a practical example. We can walk through a demonstration I did on YouTube
In the video, we are given the following system of equations:
-4r - s + 3t = 9 - s
3r + 2s - t = 3 - s
r + 3s - 5t = 29
Here is the step-by-step breakdown of how to tackle this using the elimination method:
1. Pair Up and Eliminate
The goal is to shrink this 3x3 system into a more manageable 2x2 system by eliminating one variable entirely. The instructor decides to group the top two equations together, and the bottom two equations together. He targets the variable r for elimination.
2. Scale the Equations
To eliminate r from the first two equations, they need opposite coefficients. By multiplying the top equation by 3 and the middle equation by 4, the r terms become -12r and 12r. Adding these new equations together eliminates r and, after combining the lingering s variables, leaves a new two-variable equation: 12s + 5t = 39.
He repeats this process for the second and third equations. By multiplying the third equation by -3 and adding it to the middle equation, r is eliminated again, yielding a second two-variable equation: -6s + 14t = -84.
3. Solve the New 2x2 System
Now, we just have 12s + 5t = 39 and -6s + 14t = -84. By multiplying the second equation by 2, the s terms become -12s, which perfectly cancels out the positive 12s when added together. This leaves an equation entirely in terms of t, which solves to t = -43/11.
4. The Domino Effect (Substitution)
Once you have the first domino (t), the rest fall easily. You plug t back into one of your 2-variable equations to find s. Finally, take both s and t and plug them back into any of the original three equations to solve for your last unknown, r.
Tips and Tricks for Saving Time
If you want to solve these quickly without tearing your hair out, keep these tips in mind:
Simplify First: In the video example, there were s variables hanging out on both the left and right sides of the equals sign. While the instructor waited to combine them later, a massive time-saver is to combine all like terms and move all variables to one side before you do any multiplication or elimination.
Target the Easy Variable: Always scan your equations for a variable that has a coefficient of 1 or -1 (like the r in the third equation or the s in the first). It is infinitely easier to multiply a 1 to match another number than it is to find a common multiple for numbers like 7 and 5.
Keep Your Columns Neat: Write your r's over your r's, and your t's over your t's. A single dropped negative sign or misaligned column will ruin the entire calculation.
The Mental Benefits of the Math Grind
Why do we force students to solve these? Beyond passing a test, solving 3x3 systems provides an incredible workout for your brain's executive functioning:
Working Memory Expansion: You are forced to hold multiple moving parts in your head at once—remembering which equation you are substituting into while tracking negative signs.
Rigor and Attention to Detail: These problems are unforgiving. A tiny mistake in step one compounds by step five. Practicing this builds a habit of double-checking work and maintaining intense, hyper-focused attention.
Cognitive Stamina: In an era of 15-second short-form videos, sitting down to execute a 10-minute, multi-step logical process trains patience and resilience. You learn not to panic when the answer isn't immediately obvious, which is a wildly valuable skill in both academics and life!
Friday, March 27, 2026
The Blueprint of Modern Analysis: Borel Sets, sigma-algebras, and Integration
Modern probability and advanced calculus rest on foundations
much deeper than the simple integration rules taught in calculus. The
transition from the classical Riemann integral to the robust Lebesgue
integral requires redefining how we perceive "size" and how
we classify subsets of a space.
This post will trace the hierarchy of concepts that build
modern analysis, moving from fundamental topology to the rigorous definition of
a probability space.
1. The Foundation: Topology and Allowable Shapes
Before we can calculate the "size" of something,
we must define what we are allowed to work with.
- Topology is
the study of openness and proximity. On the real number line, the topology
defines open intervals, like the set of all numbers x such
that 0 < x < 1. Topology gives the space its
"structure" regarding limits and continuity.
- In
integration, we need a mathematical structure that is broader than just
open intervals. We need a structure that allows us to combine, subtract,
and complement sets while remaining "measurable." This structure
is called a sigma-algebra.
A sigma-algebra is a collection of subsets
of a larger space X that is closed under complements, countable
unions, and countable intersections. If A is in the sigma-algebra,
then "not A" must also be.
- The Borel sigma-algebra (which
contains all Borel Sets) is the smallest possible sigma-algebra
generated by a topology. It includes all open sets, all closed sets, and
anything you can create by applying a countable number of unions and
intersections. Almost every set a practitioner works with (like points,
intervals, and open/closed shapes) is a Borel set.
2. Measure Theory and Probability
Once we have a universe (the set X) and a collection of
allowable subsets (the sigma-algebra, F), we can introduce a
"measuring stick."
A Measure (mu) is a function that assigns a
non-negative numerical "size" to any set in the sigma-algebra.
If we have a sequence of disjoint sets (sets that do not overlap), the measure
of their total union is exactly the sum of their individual measures. This
property is known as countable additivity.
- Probability
Theory is exactly measure theory, but with one critical
normalization constraint. We define the total measure of the entire
universe as exactly 1.
The foundational "Probability Space" is a
triple (Omega,F, P):
- Omega is
the Sample Space (the universe of all possible outcomes).
- F is
the sigma-algebra of Events (all subsets of outcomes
we can assign probabilities to).
- P is
the Probability Measure (the function assigning values from 0 to 1 to the
events).
3. The Power of Lebesgue Measure
Classical length is defined by the interval: the length
of [a, b] is b - a. However, this definition fails us when sets
are exceptionally complex or "pathological."
The Lebesgue Measure (lambda) is the
extension of this natural concept of "length" (or "volume"
in higher dimensions) to all Borel sets.
The critical insight of measure theory is that we also
include all subsets of "null sets" (sets with measure 0). The
collection of all Borel sets plus these null sets forms
the Lebesgue sigma-algebra, which is strictly larger than the
Borel sigma-algebra.
4. Lebesgue Integration
The Riemann integral learned in calculus
partitions the x-axis (the domain). If we have a very complex function (like
one that is 1 on all rational numbers and 0 on irrational numbers),
partitioning the x-axis fails completely.
The Lebesgue integral flips this process.
Instead of partitioning the domain, it partitions the y-axis (the range).
Instead of "height times width," the Lebesgue
integral sums "height times the measure of the set where
the function has that height."
This is written in full notation as the integral of f over
the space X, with respect to the measure mu:
\int_X f \, d\mu
This definition handles discontinuities gracefully. Because
the measure of all rational numbers on the real line is 0, the Lebesgue
integral of the "1 on rationals, 0 on irrationals" function is simply
0. Riemann's definition cannot reach this consistent conclusion.
This robust definition of integration is necessary for
establishing powerful results like the Monotone Convergence and Dominated
Convergence theorems, which allow us to interchange limits and integral
signs—the core operation in virtually all serious analytic proofs and advanced
probability.
Friday, February 20, 2026
Comparison: Variation of Parameters vs. Undetermined Coefficients
Building the Math Mentor Framework Together: What Are We Missing?
In my last post, we explored what happens when we apply Tim Ferriss’s Tribe of Mentors interview framework to the mathematics classroom. Ad...
-
What does it really mean for students to be mathematically fluent? If you’ve been in any math PD over the past few years, you’ve likely hea...
-
It is easy to get caught up in the vague, existential questions of education: How do I become a better teacher? How do I make math engaging?...
-
If you have ever taken a calculus class, you probably remember the Taylor series. It is the mathematical magic trick that lets you turn comp...