Mathematical induction Created 2026-09-24 Updated 2026-10-14
A mathematical induction proof establishes a property for every integer by proving its induction base case and its inductive step. If a counterexample existed, the well-ordering principle would give a least one; its predecessor and the inductive implication exclude it. The starting case excludes a counterexample at .
The useful idea is to cover infinitely many cases with two finite pieces of reasoning. The induction base case starts the chain, and the inductive step shows that success at an arbitrary stage propagates to the next. It is especially natural when a construction adds one item at a time: another row of a diagram, another input processed by an algorithm, or another term of a sequence. The loop invariant method is the same propagation idea applied to successive stages of a computation.
For the triangular number , the proposed formula is . At , both sides are one. Assuming the formula for a fixed positive integer , add the next row:
This is the desired formula at , proving the inductive step. For example, and ; these numerical cases illustrate the general step rather than replace it.
Figure 1.
Adding the next row gives the induction step; rearranging two triangular arrays gives the same formula geometrically
. Details.
The second picture explains where the formula comes from: two copies of the triangular arrangement tile an by rectangle. Thus . The mathematical induction argument certifies how the formula continues, while the rearrangement makes its factors visible. Neither picture is restricted to the five rows drawn: both constructions are specified for every positive integer .
In a Boolean matrix, find an axis-aligned rectangle containing only true entries and having maximum area. One approach first records each row's consecutive true run ending at each column. For each possible top-right corner, extend downwards, maintaining the minimum run width and maximizing height times width. This uses time complexity on an by matrix after preprocessing.
Midpoint quadrature rule Created 2026-10-03 Updated 2026-10-14
The midpoint quadrature rule replaces the integrand on , , by its value halfway across the interval:
Geometrically, is the signed area of a rectangle with width and height . Its single quadrature node is and its quadrature weight is . Sampling in the interior is useful when endpoint evaluation is inconvenient, although an endpoint singularity can invalidate the quadrature error estimates below.
The midpoint quadrature rule is exact for every affine function : the average height over an interval equals its midpoint height. For example, on , gives and . Positive and negative area discrepancies cancel for a sloping straight line. For the polynomial on , however, and , giving approximation error . Thus its degree of precision of a quadrature rule is exactly one.[ref]
Figure 1.
Midpoint rectangles and trapezoids approximate the area under a quadratic
. Details.
For the polynomial on , four panels give composite midpoint quadrature rule value and composite trapezoidal rule value , bracketing the exact integral . Individual midpoint rectangles can cross the curve; it is their total signed area that underestimates the integral.
For , meaning that its first two derivatives exist and are continuous functions on the closed interval, the signed approximation error has the form
Hence a convex function with is underestimated, and a concave function with is overestimated. A short justification uses the Taylor formula with integral remainder: writing , it gives , where for and for . This nonnegative Peano kernel has integral , so the weighted mean value theorem for integrals gives the stated sign and coefficient.
Minkowski addition Created 2026-10-06 Updated 2026-10-14
The Minkowski sum of subsets of a vector space collects the result of vector addition for every possible pair:
Think of fixing a point of , making a translation of the whole of by , and then taking the union of these copies over every . The points of specify displacement vectors; they need not be the centre of a shape. Minkowski addition applies to nonconvex sets too. It is different from set union: for on the real line, but .
For axis-aligned rectangles in the Euclidean plane, the coordinate intervals add separately:
Every coordinate between and can be split into one coordinate in and another in , and likewise for the vertical coordinate. Thus a rectangle plus a rectangle gives a square of area , although the two input areas total only . The Minkowski sum adds possible displacements, so its area need not be the sum of input areas.
For closed Euclidean disks centred at the origin, the triangle inequality gives . Conversely, if , split any as . Therefore , including the case . Disks of radii and give a disk of radius , with areas . For disks with centres , the resulting centre is . These two examples lead to the root-of-area comparison in the Brunn-Minkowski theorem.
Figure 1.
Minkowski addition of two rectangles and of two closed disks: exact toy shapes, sample displacement arrows, and the square-root-of-area comparison
. Details.
Minkowski addition is commutative and associative because vector addition is. Its identity is the singleton , and . For displacement vectors and , translations and uniform dilations satisfy
where uses scalar multiplication. A Minkowski sum of convex sets is convex: a line segment between and is obtained by adding the corresponding points of the line segments in and . Closure requires care; see sum of two closed sets need not be closed and sum of a compact set and a closed set.
Suppose a robot occupies a closed set relative to its chosen reference point, and a closed obstacle occupies . Keep the robot's orientation fixed and allow only translations. At reference position , its occupied set is . If touching counts as collision, the forbidden positions in this configuration space are exactly
The minus sign reflects the robot through its reference origin: solving gives . Thus collision checking for a moving shape becomes checking whether one point lies in a fixed Minkowski sum. Steven M. LaValle derives this translational configuration space construction in Planning Algorithms, section 4.3.2, equation 4.37 and Example 4.13.
For a concrete toy example, take the robot rectangle with its reference point at the lower-left corner, and obstacle rectangle . Their collision region is the rectangle
For example, is forbidden: meets in , of area . At there is just corner contact, which our convention also forbids. Using instead would wrongly miss both reference positions. Only when , as for an origin-centred disk, does the reflection leave the robot set unchanged.
Figure 1.
Toy translational collision geometry: a corner-referenced robot requires , of area , and its positive Minkowski sum puts forbidden positions on the wrong side
. Details.
The same calculation handles bounded additive uncertainty. If the actual reference position is for any , the union of possible occupied sets is . There exists an allowed error causing contact exactly when
Consequently avoiding this set guarantees avoidance for every in this model. For , the example's forbidden region becomes , with area . Here because the error bounds are symmetric; a one-sided error set must still be negated. This is a worst-case geometric bound, not a probability of collision. Rotating the robot requires an orientation-dependent Minkowski sum for each orientation.
Multiplication Created 2026-09-24 Updated 2026-10-13
Multiplication combines two factors into their product of numbers. For a natural number , is the sum of copies of ; with no copies the result is zero. This repeated-addition interpretation directly explains integer grouping, while multiplication of arbitrary real numbers also describes scaling by noninteger factors.
A rectangle of side lengths and has area . Split the second side into lengths and . The two areas are and , giving the visible identity
Rotating the same rectangle interchanges the sides without changing its area, explaining . These pictures illustrate distributivity and commutativity for ordinary positive numbers; matrix multiplication and other operations require their own rules.
Figure 1.
Partitioning a rectangle explains distributivity; partitioning a unit square gives an exact product of fractions
. Details.
For a fractional example, take of a quantity and then of that result: the combined scale is . The two fractions are multiplicative inverses. In the unit square on the right, the overlap of a strip occupying of the width with one occupying of the height contains six of twelve equal cells. Its area is . The product stays exact even though the decimal expansion of repeats. [ref]
Non-isolated point 2026-10-25
A non-isolated point of a topological space is one for which is not open, equivalently every neighbourhood of contains another point of . A homeomorphism preserves this property because it and its inverse preserve open sets, in particular open singletons. In a product space , a point is isolated exactly when both coordinates are isolated: use basic open rectangles and fix the other coordinate. Thus the non-isolated points of are , where the primes denote the non-isolated points in the spaces themselves.
Numerical integration Created 2026-10-03 Updated 2026-10-18
Numerical integration approximates a definite integral using finitely many sampled function values. For an integrand on a finite interval , , a quadrature rule has the form
The quadrature nodes tell us where to evaluate ; the quadrature weights turn those values into an approximation to the accumulated signed area. The exact definite integral depends on the integrand throughout the interval, whereas the quadrature rule sees only these samples. Integrating a chosen interpolant exactly still generally approximates the original integrand's integral.[ref]
The fundamental theorem of calculus evaluates as when an antiderivative is available. Numerical integration is useful when no convenient antiderivative is known, or each evaluation of the integrand requires an expensive computation. The task then becomes choosing informative quadrature nodes and controlling the quadrature error with as few evaluations as possible. Fixed data also permit a quadrature rule, but refinement cannot supply new information between observations.[ref]
For an exact comparison, take the polynomial on . Its antiderivative gives . On a partition of an interval into equal panels, the composite midpoint quadrature rule uses midpoint samples, each with quadrature weight ; the composite trapezoidal rule uses endpoint samples, with outer quadrature weights and interior quadrature weights . The existing derivations under these two quadrature rules give
In particular, and already bracket ; two panels give and . The convex function explains the opposite signs. A midpoint rectangle can cross the curve within its panel, so its total signed area, rather than pointwise ordering, determines the underestimate.
Figure 1.
One and two panels for numerical integration of : midpoint and endpoint samples with their quadrature weights
. Details.
Define the signed approximation error as . In this example the absolute errors are and : doubling reduces each by exactly four. A plot with both axes on a logarithmic scale therefore has slope , expressing order two with respect to panel refinement. This is an exact calculation for this polynomial, including , rather than an observed rate extrapolated from a few computed answers.
Figure 2.
Exact absolute errors for midpoint and trapezoidal numerical integration of : slope and a factor-four reduction
. Details.
For a general integrand, the composite midpoint quadrature rule and composite trapezoidal rule give second-order quadrature error bounds when has two continuous derivatives on and a known bound on . These hypotheses describe behavior between samples. A formal error formula containing an unknown derivative is not a numerical accuracy certificate or a justified sample count. Finite samples do not certify quadrature accuracy makes the missing information explicit. Errors in the evaluations and roundoff error must also be considered; noise sensitivity of a quadrature rule treats the former.
More elaborate choices exploit other information: Gaussian quadrature obtains high polynomial exactness from carefully chosen quadrature nodes, and Simpson's rule uses a quadratic interpolant. Adaptive numerical integration compares approximations and concentrates samples where they disagree, often saving evaluations. Such a difference is an error estimate, not a certified bound without further assumptions: both approximations can miss the same feature.[ref]
Objective function 2026-10-13
An objective function is a real-valued function assigning a score to each decision under consideration. An optimization problem seeks a feasible solution maximizing or minimizing this score. The optimizing decision and its score are different objects: in the budget constraint example, the decision is and the optimal value is .
Choosing an objective function expresses what matters. With and , maximizing makes every feasible point equally good; maximizing selects ; maximizing rewards balance. If are side lengths, is a rectangle's area, and the budget constraint fixes its perimeter at . The same mathematics then says that a square encloses the largest area among these rectangles.
For the product objective function, the positive level sets are hyperbolas. Each collects decisions with identical scores. The highest attainable level set touches the feasible set at ; higher ones miss it entirely.
Figure 1.
Product contours meet the feasible budget segment, with a unique maximum of 25 at (5,5)
. Details.
The product is neither a convex function nor a concave function on the positive quadrant: along it is , while along it is . Yet its restriction to this feasible set is a strictly concave function. Thus the geometry of the allowed decisions can make an otherwise nonconcave maximization problem tractable. Replacing by always converts a maximization problem into an equivalent minimization problem.
On two copies of the real line with their usual topologies, let and . The coordinate projection maps give
The first two open sets are unbounded strips; their set intersection is an open rectangle. The point passes both coordinate tests, whereas fails the first because the endpoint is excluded.
Figure 1.
An open rectangle is the intersection of two coordinate-projection strips
. Dashed boundaries and hollow endpoints are excluded. Arrows at the ends of the strips indicate that the strips continue indefinitely. The final panel projects the interior point onto the two tested intervals. Details.
Every point of an open set in the Euclidean plane has such an open rectangle around it inside that open set. Conversely, every open rectangle contains a small Euclidean ball around each of its points. Thus the product topology here is the usual metric topology of the Euclidean plane. With finitely many factors all coordinates may be restricted; on an infinite Cartesian product a basic cylinder set must leave all but finitely many factors free.
The Young diagrams of partitions in a rectangle with at most rows and at most cells in each row correspond to the lattice paths in the lattice-path interpretation of Gaussian binomial coefficients. Read the shaded row lengths from top to bottom; the result is a nonincreasing list of nonnegative integers. This bijection preserves the number of cells.
In a two-by-two rectangle, the possibilities are the empty integer partition, , , , , and . Their sizes are , giving . In particular, there are two allowed diagrams of size two, explaining the middle coefficient.[ref]
Take the lower coordinate in each direction from the smaller lower bound, and the upper coordinate from the larger upper bound. This gives an axis-aligned bounding box, rather than the generally nonrectangular set union of the two regions. Use defensive copying because the paper's Point has mutable fields.
final class Rect {
    private final Point lowerLeft, upperRight;

    public Rect(Point ll, Point ur) {
        if (ll == null || ur == null)
            throw new IllegalArgumentException("null corner");
        double xl = ll.getX(), yl = ll.getY();
        double xh = ur.getX(), yh = ur.getY();
        if (!Double.isFinite(xl) || !Double.isFinite(yl)
                || !Double.isFinite(xh) || !Double.isFinite(yh)
                || xl > xh || yl > yh)
            throw new IllegalArgumentException("invalid rectangle");
        lowerLeft = new Point(xl, yl);
        upperRight = new Point(xh, yh);
    }

    public Rect Union(Rect r) {
        return new Rect(
            new Point(Math.min(lowerLeft.getX(), r.lowerLeft.getX()),
                      Math.min(lowerLeft.getY(), r.lowerLeft.getY())),
            new Point(Math.max(upperRight.getX(), r.upperRight.getX()),
                      Math.max(upperRight.getY(), r.upperRight.getY())));
    }

    public Point lowerLeft() { return lowerLeft.clone(); }
    public Point upperRight() { return upperRight.clone(); }
}
Any containing rectangle must reach at least as far as each input rectangle in all four coordinate directions. The constructed rectangle reaches exactly those four limits, which proves both containment and minimality:
The corner accessors return copies, so changing a constructor argument or returned Point cannot alter an existing Rect.
For each sample, include both its marker and its bar. A bar's lower vertical limit is , including negative bars below the horizontal axis. Its horizontal half-width is , rather than . Its marker extends by in both directions. Therefore that sample's axis-aligned bounding box is
Taking successive bounding unions includes every sample. Finally unite with to include the stipulated axes. The printed upper-corner formula has one extra closing parenthesis; both the lower-corner formula and the picture confirm the coordinatewise interpretation used here.[ref]
The following commented Java object method may be placed in a graph-view utility software class. It uses Rect from part (a) and one iterator snapshot from part (e).
static Rect barBounds(DataList data, double C, double M) {
    if (!Double.isFinite(C) || !Double.isFinite(M) || C <= 0 || M <= 0)
        throw new IllegalArgumentException("positive finite margins required");
    List<Point> p = new ArrayList<Point>();
    for (Point point : data) p.add(point);
    int n = p.size();
    if (n < 2) throw new IllegalArgumentException("at least two points required");

    // This central square supplies the required minimum axis extents.
    Rect bounds = new Rect(new Point(-M, -M), new Point(M, M));
    for (int i = 0; i < n; ++i) {
        double x = p.get(i).getX(), y = p.get(i).getY();
        double width;
        if (i == 0) {
            width = p.get(1).getX() - x;
        } else if (i == n - 1) {
            width = x - p.get(i - 1).getX();
        } else {
            width = Math.min(x - p.get(i - 1).getX(),
                             p.get(i + 1).getX() - x);
        }
        // Include whichever is wider: marker or half of the bar.
        double half = Math.max(C, width / 2.0);
        // The bar reaches the baseline; the marker reaches y +/- C.
        Rect item = new Rect(new Point(x - half, Math.min(0.0, y - C)),
                             new Point(x + half, Math.max(0.0, y + C)));
        bounds = bounds.Union(item);
    }
    return bounds;
}
This uses the usual assumption that coordinate differences and expanded bounds are representable as finite double values. Rect rejects overflowed bounds; numerical range extension is not part of the requested geometry. As with all floating-point arithmetic, the calculation approximates the corresponding exact real number formulas.
Every side produced by the coordinatewise extrema is required by some bar, marker enclosure or axis extent. The result is therefore the smallest guaranteed rectangle for the specified marker enclosures and axes. If denote the graph-only bounds, the final answer is
Choose the holomorphic square root in the complex upper half-plane whose boundary value is positive on . Its primitive has derivative . The integral from the boundary basepoint zero is well-defined by a limiting path: the local singularity is only .
The three finite prevertices have derivative exponent , hence image angle . At infinity , so the local coordinate also gives image angle . These are the four corners of a Schwarz-Christoffel mapping to a rectangle. To establish that it is a square rather than a general rectangle, put
The length along is also by . The length along is
by , and the fourth side has the same length by symmetry. With the chosen branch, the boundary derivatives on the intervals have phases , respectively. Therefore
The extended real boundary maps continuously and once around this square. Each side is traversed monotonically, and the four corner limits are finite. The argument principle, applied after small indentations around the boundary singularities, counts one inverse image of each interior value and zero of each exterior value. Since never vanishes in the complex upper half-plane, those interior inverse images are simple. Thus
Take a real modulus and complementary modulus . The Jacobi elliptic sine is constructed by inverting a Schwarz-Christoffel mapping. The main reason for the Schwarz-Christoffel derivative is local angle behavior: if a prevertex corresponds to a polygon angle , straightening that corner gives . Hence has exponent . Along each straight boundary side its argument is constant. Dividing by removes those angle changes; Schwarz reflection extends the quotient across the real boundary and the corner singularities. With all vertices finite, its behavior at infinity is also regular, because the sum of the derivative exponents is . It is therefore a nonzero constant on the sphere. This yields
This argument explains the exponents, the branch choices and the role of boundary reflection, rather than just writing down the formula.
For a rectangle, four angles are , so choose prevertices and normalize the inverse map as
Use the branch positive on . Set
The two middle corner values are . The right side length is
as the substitution shows. The selected branch makes the increment along that side . Thus the other corner values are and . The two boundary intervals through infinity have combined horizontal length ; the substitution gives length from to infinity. In particular . As in the square mapping, the boundary goes once around a rectangle, and the argument principle gives a conformal bijection
Its inverse maps that rectangle onto the complex upper half-plane, and near zero has , . Differentiating the inverse relation gives
These identities continue meromorphically.
All four rectangle edges have real images, so Schwarz reflection extends across them. Reflection in the bottom edge gives . Reflection in the right edge, combined with the first reflection, gives . The even inverse-integral derivative gives , and hence
Reflection in the top edge and then the bottom gives . The resulting period lattice of the Jacobi elliptic sine is
To see the poles, use near the boundary point infinity. Then , so the inverse has a simple pole at with residue . The real half-period shift gives another pole at with opposite residue. Reflection tiles the plane with copies of the rectangle, so these are exactly two simple poles per displayed fundamental cell. They also show there is no additional period: a period must permute these two pole classes; exchanging them would shift by modulo the displayed lattice, which reverses the function's sign and residues rather than preserving it.
For any finite , integrate around a fundamental parallelogram avoiding its zeros and poles. The opposite edges cancel by periodicity, so the number of zeros of equals the number of poles, namely two. Therefore every value has two inverse images modulo periods, counted with multiplicity. The value infinity likewise has the two simple poles as preimages. Usually the finite-value preimages are distinct; the involution interchanges them. At , a critical point supplies one double inverse image instead. For example , , . Thus the printed “exactly two” requires the standard multiplicity convention; it would be false if interpreted as two distinct solutions at these branch values of a holomorphic map.
Now choose and put . Since , is even, elliptic for , and has degree two. The Weierstrass elliptic function for the same lattice is even and has one double pole per cell, hence also degree two. Its generic fibers are precisely . Consequently an even is constant on those fibers and descends to a meromorphic function of . At the four fixed classes of , the even local Laurent expansion makes this descent meromorphic in the squared local coordinate. A meromorphic function on the Riemann sphere is rational, so . Generic map degrees multiply:
Thus has degree one and the degree-two even elliptic function conclusion is
The determinant condition follows because is nonconstant. This proves the existence requested, with the explicit allowable shift .
Let be the difference of two solutions. It is a harmonic function in the open rectangle, continuous on its closure, with
Interpret these Neumann boundary conditions classically on the open edges, with the corresponding first derivatives extending there. We use the strong maximum principle for harmonic functions together with even reflection across a homogeneous Neumann boundary. This avoids imposing differentiability at the two bottom corners, where the original data fail a corner compatibility condition for mixed boundary data.
By continuity on the compact rectangle, attains a maximum . Suppose . It cannot occur on the bottom edge, where . If it occurs in the interior, the strong maximum principle for harmonic functions makes constant, contradicting its bottom value. If it occurs on an open vertical or top edge, the zero normal derivative allows an even reflection across that edge. The reflected harmonic function then has an interior maximum, giving the same contradiction.
The only remaining possibilities are the two top corners. Reflect evenly across both meeting edges. This gives a bounded harmonic function on a punctured neighbourhood of the corner. The removable singularity for a bounded harmonic function fills the centre with its continuous value. The corner is now an interior maximum, again impossible. Therefore . Applying the identical argument to gives .
Thus , and the solution is unique among the continuous classical solutions with the stated edge conditions.
For separation of variables, set , with . The side Neumann boundary conditions give . As in the standard Sturm-Liouville problem, integration by parts shows . The zero eigenvalue gives constant , and together with makes constant too.
For positive eigenvalues, the side conditions give , , . The top condition chooses the hyperbolic cosine factor . Normalize it to one at , obtaining
The orthogonality of integer Fourier modes and the permitted interchange of summation and integration give
For the remaining Fourier coefficients, integrate by parts twice. The sine boundary term vanishes, and with ,
The last integral vanishes because and is a positive integer. Thus : only even Fourier modes survive, as expected from reflection symmetry about . Writing gives
The factor involving hyperbolic cosine lies between zero and one on the closed rectangle. Comparison with proves uniform convergence, so the solution is continuous on the entire boundary and takes the prescribed bottom values. On every strip , the Fourier coefficients have additional exponential decay, allowing differentiation and verifying Laplace's equation and all three Neumann boundary conditions.
To justify the endpoint calculation in b, differentiate at , in the interior-edge sense:
Subtracting from the hyperbolic tangent leaves an absolutely convergent series. For the remaining part, the power series gives, by taking real parts at and then ,
The same formula follows as the limit of the differentiated interior solution, since its leading multiplier is . This establishes the boundary normal derivative used in b. It also explains why the failure of the corner compatibility condition for mixed boundary data affects first derivatives at the bottom corners without spoiling continuity of itself.
The inequalities become the rectangle , in parabolic polar coordinates. The change of variables formula and the Jacobian determinant from a give the area
The middle panel in b shows the curved region in Cartesian coordinates; its straight edges in parabolic polar coordinates explain the simple integration limits.
Use the change of variables , . Its Jacobian determinant in the first quadrant is , and the inverse is uniquely determined by
Both are positive for , so the transformation is one-to-one on the region. It becomes the rectangle , . The integrand is , and . Thus
Parametrize the function graph by . Its two tangent vectors are and . A small parameter rectangle maps, to first order, to a tangent parallelogram with area . Since
the surface area of a graph is
For example, ensures a regular surface chart; more general domains can be handled by the corresponding convergent improper integral.