Fundamental Sets and General Solutions

Consider the equation y′′+p(t)y′+q(t)y=0y'' + p(t)y' + q(t)y = 0 on an interval where pp and qq are continuous. Two students each guess a solution — say y1(t)y_1(t) and y2(t)y_2(t) — and both are correct. A natural question follows: have they found all solutions, or only a sliver of them? This is the central organizing question of linear ODE theory. Answering it precisely is what makes the subject tractable, rather than an endless search for one-off solutions. The answer turns out to be sharply binary. Either y1y_1 and y2y_2 are enough. Every solution is some combination c1y1+c2y2c_1y_1 + c_2y_2. Or they are redundant, secretly proportional, and miss almost everything. The dividing line is exactly linear independence. The tool that certifies it is the Wronskian. This article proves the theorem that upgrades "a pair of independent solutions" into "a complete description of the solution space": the fundamental set theorem. The big picture is that the solution set of a linear homogeneous ODE of order nn is not a wild, unstructured collection of functions. Instead it is a vector space of dimension exactly nn. Once you know this, solving the ODE becomes a finite, algebraic task. Find nn independent solutions, and you are done. Every other solution is a guaranteed linear combination of them, no exceptions, no missed cases. This is the payoff of building the theory carefully, rather than guessing solutions one at a time. This article states and proves the fundamental set theorem for second-order linear homogeneous equations. The argument generalizes verbatim to order nn. It builds the intuition for why the Wronskian is the exact right certificate of completeness. It also works through examples where the solutions are non-obvious functions rather than clean exponentials. It closes by connecting this structure to the constant-coefficient methods and higher-order theory that build on it directly.

Key takeaways
  • A fundamental set of solutions is a linearly independent set of nn solutions to an nnth-order linear homogeneous ODE.
  • The general solution is the set of all linear combinations of a fundamental set — nothing more is needed.
  • Linear independence of solutions on an interval II is certified by the Wronskian being nonzero at even one point of II.
  • The Existence–Uniqueness Theorem for linear equations is the engine that forces every solution to be a combination of the fundamental set.
  • The solution space of an nnth-order linear homogeneous ODE is a genuine vector space of dimension nn.
  • Key idea: finding nn independent solutions is not a partial victory — it is total victory, because the theorem guarantees they span everything.

Where this lives today

The fundamental set theorem is the reason engineers can write "the general solution of a damped spring-mass system is c1er1t+c2er2tc_1e^{r_1t}+c_2e^{r_2t}" and mean it exactly, not approximately — this underlies vibration analysis, circuit theory (RLC circuits), and control systems, all covered under Homogeneous Equations with Constant Coefficients. In quantum mechanics, the time-independent Schrödinger equation for simple potentials is a second-order linear ODE, and identifying a fundamental set of solutions is precisely how bound-state and scattering solutions are classified. In numerical analysis, the theorem justifies why solving an nnth-order linear ODE numerically only ever requires tracking nn independent trajectories (or an nn-dimensional state vector), which is the basis of state-space methods and the transition from scalar ODEs to first-order systems discussed under Higher-Order Linear Equations. In algebraic and differential Galois theory, the dimension of the solution space of a linear ODE (a D\mathcal{D}-module) is a basic invariant used to classify special functions; hypergeometric and Bessel equations are studied precisely by exhibiting fundamental sets and their monodromy. Finally, in statistics and signal processing, the state-space (Kalman filter) formulation of linear time-varying systems relies on exactly this vector-space structure of the solution set to propagate uncertainty consistently through time.

:::quiz[If y1,y2y_1,y_2 solve a linear homogeneous 2nd-order ODE with continuous coefficients on II and W[y1,y2](t1)=0W[y_1,y_2](t_1)=0 for one t1∈It_1\in I, what can you conclude about WW at other points of II?]

  • It must be nonzero everywhere else in II
  • It is identically zero on all of II ✓
  • It could be zero at some points and nonzero at others
  • Nothing can be concluded without more information

:::quiz[The Cauchy–Euler equation 2t2y′′−ty′+y=02t^2y''-ty'+y=0 has solutions y1=ty_1=t, y2=t1/2y_2=t^{1/2} on (0,∞)(0,\infty). What is W[y1,y2](4)W[y_1,y_2](4)?]

  • 00
  • −1-1 ✓
  • 22
  • 4=2\sqrt{4}=2, so undefined sign

:::quiz[For the repaired fundamental set {e2t,te2t}\{e^{2t}, te^{2t}\} of y′′−4y′+4y=0y''-4y'+4y=0, what is W[e2t,te2t](0)W[e^{2t},te^{2t}](0)?]

  • 00
  • 11 ✓
  • 44
  • e4⋅0=1e^{4\cdot0}=1 but with a minus sign, so −1-1

Ask the tutor: Why does a nonzero Wronskian at just one point guarantee it's nonzero everywhere on the interval, and how does that fact depend on the equation being linear and homogeneous?

DefinitionFundamental set of solutions and general solution

Let p,qp, q be continuous real-valued functions on an open interval II, and consider the second-order linear homogeneous equation

y′′+p(t)y′+q(t)y=0,t∈I.(∗)y'' + p(t)y' + q(t)y = 0, \qquad t \in I. \tag{$\ast$}

A pair of solutions y1,y2y_1, y_2 of (∗)(\ast) on II is called a fundamental set of solutions if y1y_1 and y2y_2 are linearly independent on II, equivalently if their Wronskian

W[y1,y2](t)=y1(t)y2′(t)−y1′(t)y2(t)W[y_1, y_2](t) = y_1(t)y_2'(t) - y_1'(t)y_2(t)

satisfies W[y1,y2](t0)≠0W[y_1,y_2](t_0) \neq 0 for some (equivalently, by The Wronskian and Linear Independence, every) t0∈It_0 \in I. Given a fundamental set {y1,y2}\{y_1, y_2\}, the general solution of (∗)(\ast) on II is the family

y(t)=c1y1(t)+c2y2(t),c1,c2∈R arbitrary constants,y(t) = c_1 y_1(t) + c_2 y_2(t), \qquad c_1, c_2 \in \mathbb{R} \text{ arbitrary constants,}

and the claim to be proved is that this family is exactly the solution set of (∗)(\ast) on II: every solution of (∗)(\ast) arises this way for a unique choice of (c1,c2)(c_1, c_2). The definition extends verbatim to order nn: an nnth-order linear homogeneous equation y(n)+pn−1(t)y(n−1)+⋯+p0(t)y=0y^{(n)} + p_{n-1}(t)y^{(n-1)} + \cdots + p_0(t)y = 0 has a fundamental set consisting of nn linearly independent solutions y1,…,yny_1,\dots,y_n, and its general solution is y=c1y1+⋯+cnyny = c_1y_1 + \cdots + c_ny_n.

Intuition

Think of the solution set of (∗)(\ast) as a landscape you are trying to map. Each individual solution is a single point discovered by trial, substitution, or luck. The fundamental set theorem says this landscape is flat and two-dimensional (for a second-order equation): once you have two independent "directions" — two solutions that are not scalar multiples of each other — you can reach every point in the landscape by walking a straight combination of those two directions. There is no hidden third direction, no exotic solution lurking outside the span. This is a strong structural claim, and it is exactly what makes linear ODEs solvable in a systematic way rather than an art of clever guesses. The Wronskian is the instrument that detects whether two solutions are truly two independent directions or secretly the same direction in disguise. Geometrically, (y1(t),y1′(t))(y_1(t), y_1'(t)) and (y2(t),y2′(t))(y_2(t), y_2'(t)) are two points in the phase plane R2\mathbb{R}^2 at each time tt; the Wronskian is (twice) the signed area of the parallelogram they span. If that area is ever nonzero, the two vectors are not parallel, and — as proved under The Wronskian and Linear Independence — this non-parallelism persists for all t∈It \in I whenever y1,y2y_1, y_2 solve (∗)(\ast). The fundamental set theorem then says: two non-parallel initial vectors are enough to hit any target initial vector (y0,y0′)(y_0, y_0') by scaling and adding, and each such target initial vector corresponds to exactly one solution curve by uniqueness. Combining these two facts is the whole proof.

Figure
The Wronskian as a non-degeneracy certificate: two solution vectors in the phase plane at t0t_0, and the target vector they must reach
y y′ (y₁(t₀), y₁′(t₀)) (y₂(t₀), y₂′(t₀)) c₁·v₁ + c₂·v₂ = (y₀, y₀′) W[y₁,y₂](t₀) ≠ 0 ⇒ v₁, v₂ span the plane
TheoremFundamental Set Theorem for Second-Order Linear Homogeneous Equations

Let p,q:I→Rp, q : I \to \mathbb{R} be continuous on an open interval II, and let y1,y2y_1, y_2 be solutions of

y′′+p(t)y′+q(t)y=0y'' + p(t)y' + q(t)y = 0

on II such that W[y1,y2](t0)≠0W[y_1,y_2](t_0) \neq 0 for some point t0∈It_0 \in I. Then for every solution ϕ\phi of the equation on II, there exist unique constants c1,c2∈Rc_1, c_2 \in \mathbb{R} such that

ϕ(t)=c1y1(t)+c2y2(t)for all t∈I.\phi(t) = c_1 y_1(t) + c_2 y_2(t) \quad \text{for all } t \in I.

Consequently the solution set of the equation on II is exactly {c1y1+c2y2:c1,c2∈R}\{c_1y_1 + c_2y_2 : c_1, c_2 \in \mathbb{R}\}, a two-dimensional real vector space with basis {y1,y2}\{y_1, y_2\}.

Proof

Proof technique: existence–uniqueness for the IVP, reduced to a linear algebra step via evaluation at a single point.

Step 1 (Setup at t0t_0). Let ϕ\phi be an arbitrary solution of the ODE on II. Since p,qp, q are continuous on II and t0∈It_0 \in I, the hypotheses of the Existence and Uniqueness Theorem for linear second-order equations are satisfied on II: for any prescribed values α,β\alpha, \beta, the initial value problem

y′′+p(t)y′+q(t)y=0,y(t0)=α,y′(t0)=βy'' + p(t)y' + q(t)y = 0, \qquad y(t_0) = \alpha,\quad y'(t_0) = \beta

has a solution on all of II, and that solution is unique. Set

α=ϕ(t0),β=ϕ′(t0).\alpha = \phi(t_0), \qquad \beta = \phi'(t_0).

By construction ϕ\phi itself is a solution of this IVP.

Step 2 (Solve a 2×22\times2 linear system). We seek constants c1,c2c_1, c_2 so that ψ(t):=c1y1(t)+c2y2(t)\psi(t) := c_1y_1(t) + c_2y_2(t) satisfies ψ(t0)=α\psi(t_0) = \alpha and ψ′(t0)=β\psi'(t_0) = \beta. This is the linear system

c1y1(t0)+c2y2(t0)=α,c1y1′(t0)+c2y2′(t0)=β,\begin{aligned} c_1 y_1(t_0) + c_2 y_2(t_0) &= \alpha, \\ c_1 y_1'(t_0) + c_2 y_2'(t_0) &= \beta, \end{aligned}

whose coefficient matrix has determinant y1(t0)y2′(t0)−y1′(t0)y2(t0)=W[y1,y2](t0)y_1(t_0)y_2'(t_0) - y_1'(t_0)y_2(t_0) = W[y_1,y_2](t_0). By hypothesis this determinant is nonzero, so by Cramer's rule the system has a unique solution:

c1=αy2′(t0)−βy2(t0)W[y1,y2](t0),c2=βy1(t0)−αy1′(t0)W[y1,y2](t0).c_1 = \frac{\alpha y_2'(t_0) - \beta y_2(t_0)}{W[y_1,y_2](t_0)}, \qquad c_2 = \frac{\beta y_1(t_0) - \alpha y_1'(t_0)}{W[y_1,y_2](t_0)}.

These c1,c2c_1, c_2 are real numbers, uniquely determined by α,β\alpha, \beta and by y1,y2y_1, y_2 at t0t_0.

Step 3 (ψ\psi is a solution of the ODE). Since y1y_1 and y2y_2 are each solutions of the linear homogeneous equation and the equation is linear, any constant-coefficient combination ψ=c1y1+c2y2\psi = c_1y_1 + c_2y_2 is also a solution: substituting,

ψ′′+pψ′+qψ=c1(y1′′+py1′+qy1)+c2(y2′′+py2′+qy2)=c1⋅0+c2⋅0=0.\psi'' + p\psi' + q\psi = c_1(y_1'' + py_1' + qy_1) + c_2(y_2'' + py_2' + qy_2) = c_1 \cdot 0 + c_2 \cdot 0 = 0.

(This uses only that differentiation is linear and that y1,y2y_1, y_2 each satisfy the homogeneous equation — no circularity, since this is a direct computation, not an appeal to the theorem being proved.)

Step 4 (ψ\psi and ϕ\phi solve the same IVP). By Step 2, ψ(t0)=α=ϕ(t0)\psi(t_0) = \alpha = \phi(t_0) and ψ′(t0)=β=ϕ′(t0)\psi'(t_0) = \beta = \phi'(t_0). By Step 3, ψ\psi solves the ODE on II. So both ϕ\phi and ψ\psi are solutions of the same initial value problem

y′′+p(t)y′+q(t)y=0,y(t0)=α,y′(t0)=βy'' + p(t)y' + q(t)y = 0, \qquad y(t_0) = \alpha, \quad y'(t_0) = \beta

on the interval II. By the uniqueness clause of the Existence–Uniqueness Theorem invoked in Step 1, this IVP has at most one solution on II. Therefore

ϕ(t)=ψ(t)=c1y1(t)+c2y2(t)for all t∈I.\phi(t) = \psi(t) = c_1y_1(t) + c_2y_2(t) \quad \text{for all } t \in I.

This proves existence of the representation.

Step 5 (Uniqueness of c1,c2c_1, c_2). Suppose also ϕ=c1′y1+c2′y2\phi = c_1'y_1 + c_2'y_2 on II for some other constants. Evaluating at t0t_0 and differentiating and evaluating at t0t_0 gives the same linear system as in Step 2 with the same right-hand side (α,β)(\alpha,\beta); since that system's coefficient matrix is invertible (determinant W[y1,y2](t0)≠0W[y_1,y_2](t_0)\neq0), the solution (c1,c2)(c_1,c_2) is unique. Hence c1′=c1c_1' = c_1 and c2′=c2c_2' = c_2.

Step 6 (Conclusion). Every solution ϕ\phi of the ODE on II equals c1y1+c2y2c_1y_1+c_2y_2 for unique c1,c2∈Rc_1,c_2\in\mathbb R, and conversely every such combination is a solution by Step 3. Hence the solution set equals {c1y1+c2y2:c1,c2∈R}\{c_1y_1+c_2y_2 : c_1,c_2\in\mathbb R\} exactly. Because this set is closed under addition and scalar multiplication (Step 3's computation) and {y1,y2}\{y_1,y_2\} is linearly independent by hypothesis (via The Wronskian and Linear Independence) and spans the set (Steps 4–5), it is a two-dimensional vector space with basis {y1,y2}\{y_1,y_2\}. ■\blacksquare

∎
Corollarynn-th order generalization, stated

Let p0,…,pn−1p_0,\dots,p_{n-1} be continuous on II and let y1,…,yny_1,\dots,y_n solve y(n)+pn−1y(n−1)+⋯+p0y=0y^{(n)} + p_{n-1}y^{(n-1)} + \cdots + p_0y = 0 on II with Wronskian W[y1,…,yn](t0)≠0W[y_1,\dots,y_n](t_0)\neq 0 for some t0∈It_0\in I. Then every solution ϕ\phi on II is a unique combination ϕ=c1y1+⋯+cnyn\phi = c_1y_1+\cdots+c_ny_n, and the solution space is an nn-dimensional vector space. The proof is identical in structure: replace the 2×22\times2 linear system in Step 2 by the n×nn\times n system obtained from matching ϕ(k)(t0)\phi^{(k)}(t_0) for k=0,…,n−1k=0,\dots,n-1, whose coefficient matrix is exactly the Wronskian matrix at t0t_0, invertible by hypothesis; Steps 3–6 go through verbatim with the nn-th order Existence–Uniqueness Theorem in place of the second-order one.

Worked examples

Example 1 (Cauchy–Euler equation with an asymmetric interval and non-integer exponents). Consider on I=(0,∞)I = (0,\infty) the equation

2t2y′′−ty′+y=0.2t^2 y'' - t y' + y = 0.

Try y=try = t^r: then y′=rtr−1y' = rt^{r-1}, y′′=r(r−1)tr−2y'' = r(r-1)t^{r-2}, and substitution gives

2t2⋅r(r−1)tr−2−t⋅rtr−1+tr=[2r(r−1)−r+1]tr=(2r2−3r+1)tr.2t^2 \cdot r(r-1)t^{r-2} - t\cdot rt^{r-1} + t^r = \big[2r(r-1) - r + 1\big] t^r = \big(2r^2 - 3r + 1\big)t^r.

Factor: 2r2−3r+1=(2r−1)(r−1)2r^2 - 3r + 1 = (2r-1)(r-1), with roots r=1r = 1 and r=12r = \tfrac12. So

y1(t)=t,y2(t)=t1/2y_1(t) = t, \qquad y_2(t) = t^{1/2}

are both solutions on (0,∞)(0,\infty). Compute the Wronskian directly:

W[y1,y2](t)=y1y2′−y1′y2=t⋅12t−1/2−1⋅t1/2=12t1/2−t1/2=−12t1/2.W[y_1,y_2](t) = y_1 y_2' - y_1' y_2 = t\cdot\tfrac12 t^{-1/2} - 1\cdot t^{1/2} = \tfrac12 t^{1/2} - t^{1/2} = -\tfrac12 t^{1/2}.

For any t0>0t_0 > 0, W[y1,y2](t0)=−12t0≠0W[y_1,y_2](t_0) = -\tfrac12\sqrt{t_0} \neq 0. First put the equation in standard form to confirm continuity hypotheses: dividing by 2t22t^2,

y′′−12ty′+12t2y=0,y'' - \frac{1}{2t}y' + \frac{1}{2t^2}y = 0,

so p(t)=−12tp(t) = -\tfrac{1}{2t} and q(t)=12t2q(t) = \tfrac{1}{2t^2} are continuous on I=(0,∞)I=(0,\infty) (the singularity at t=0t=0 is excluded from II, which is exactly why we work on (0,∞)(0,\infty) rather than all of R\mathbb{R}). All hypotheses of the theorem hold, so {y1,y2}\{y_1,y_2\} is a fundamental set on (0,∞)(0,\infty) and the general solution is

y(t)=c1t+c2t1/2,t>0.y(t) = c_1 t + c_2 t^{1/2}, \qquad t > 0.

By the theorem, any solution of the original equation on (0,∞)(0,\infty) — however it was found, by series methods, numerically, or otherwise — must be expressible this way. For instance, one can now immediately solve the IVP y(4)=10y(4) = 10, y′(4)=3y'(4) = 3: substituting, 4c1+2c2=104c_1 + 2c_2 = 10 and c1+14c2=3c_1 + \tfrac14 c_2 = 3, giving (solving) c2=2c_2 = 2, c1=5/2c_1 = 5/2, so y(t)=52t+2ty(t) = \tfrac52 t + 2\sqrt t is the unique solution — uniqueness guaranteed by the theorem, not just a hopeful guess.

Example 2 (verifying two solutions are NOT a fundamental set, and repairing it). Consider y′′−4y′+4y=0y'' - 4y' + 4y = 0 on I=RI = \mathbb{R} (standard form already, p=−4p=-4, q=4q=4 continuous everywhere). Two students propose

y1(t)=e2t,y2(t)=5e2t−3e2t=2e2t.y_1(t) = e^{2t}, \qquad y_2(t) = 5e^{2t} - 3e^{2t} = 2e^{2t}.

Both are genuine solutions (check: y1′′−4y1′+4y1=4e2t−8e2t+4e2t=0y_1'' - 4y_1' + 4y_1 = 4e^{2t}-8e^{2t}+4e^{2t}=0, and y2y_2 is a constant multiple so it also satisfies the linear homogeneous equation). Compute the Wronskian:

W[y1,y2](t)=e2t⋅4e2t−2e2t⋅2e2t=4e4t−4e4t=0for all t.W[y_1,y_2](t) = e^{2t}\cdot 4e^{2t} - 2e^{2t}\cdot 2e^{2t} = 4e^{4t} - 4e^{4t} = 0 \quad\text{for all } t.

The Wronskian is identically zero, so the hypothesis of the theorem fails: {y1,y2}\{y_1, y_2\} is not a fundamental set (indeed y2=2y1y_2 = 2y_1, manifestly dependent). The theorem gives no conclusion here — and correctly so, since c1y1+c2y2=(c1+2c2)e2tc_1y_1 + c_2y_2 = (c_1+2c_2)e^{2t} only ever produces multiples of e2te^{2t}, never reaching, e.g., the solution te2tt e^{2t} (readily checked: (te2t)′′−4(te2t)′+4(te2t)=0(te^{2t})'' - 4(te^{2t})' + 4(te^{2t}) = 0). To repair this, use reduction of order or the characteristic-root method to find the second independent solution y2true(t)=te2ty_2^{\text{true}}(t) = te^{2t}. Now

W[y1,y2true](t)=e2t(e2t+2te2t)−2e2t(te2t)=e4t+2te4t−2te4t=e4t≠0W[y_1, y_2^{\text{true}}](t) = e^{2t}(e^{2t}+2te^{2t}) - 2e^{2t}(te^{2t}) = e^{4t} + 2te^{4t} - 2te^{4t} = e^{4t} \neq 0

for every t∈Rt \in \mathbb{R}, so {e2t,te2t}\{e^{2t}, te^{2t}\} is a fundamental set, and the general solution is y(t)=c1e2t+c2te2ty(t) = c_1e^{2t} + c_2te^{2t}, matching the standard repeated-root theory previewed under Homogeneous Equations with Constant Coefficients.

Counterexampledropping the nonvanishing-Wronskian hypothesis

The theorem requires W[y1,y2](t0)≠0W[y_1,y_2](t_0)\neq 0 at some point. Drop this and completeness can fail even though y1,y2y_1, y_2 both genuinely solve the equation and are formally "two different functions." Take again y′′−4y′+4y=0y'' - 4y' + 4y = 0 on R\mathbb{R}, with y1(t)=e2ty_1(t) = e^{2t} and y2(t)=2e2ty_2(t) = 2e^{2t} as in Example 2. Both solve the ODE; they are two distinct functions; but W[y1,y2]≡0W[y_1,y_2] \equiv 0. The set of combinations

S={c1e2t+c2⋅2e2t:c1,c2∈R}={ke2t:k∈R}S = \{c_1 e^{2t} + c_2\cdot 2e^{2t} : c_1, c_2 \in \mathbb{R}\} = \{k e^{2t} : k \in \mathbb{R}\}

is only a one-dimensional family, even though we wrote it with two arbitrary constants and two solutions. The genuine solution te2tte^{2t} satisfies the ODE (direct substitution: 2e2t+4te2t−4(e2t+2te2t)+4te2t=(2−4)e2t+(4−8+4)te2t=02e^{2t}+4te^{2t} - 4(e^{2t}+2te^{2t}) + 4te^{2t} = (2-4)e^{2t} + (4-8+4)te^{2t} = 0) but te2t∉Ste^{2t} \notin S: if te2t=ke2tte^{2t} = ke^{2t} for all tt, dividing by e2t≠0e^{2t}\neq0 forces t=kt = k for all tt, impossible for a constant kk. So SS is a strict, proper subset of the true solution set, despite having "two constants." This shows the nonvanishing-Wronskian hypothesis is not a technicality: without it, the representation c1y1+c2y2c_1y_1+c_2y_2 can silently omit most of the actual solution space, and nothing in the naive computation warns you that this has happened.

Scope, edge cases & common traps

SituationWhat to checkCommon error
Equation not in standard form y′′+py′+qy=0y''+py'+q y=0Divide by leading coefficient first; note where that coefficient vanishesApplying the theorem where p,qp,q are actually discontinuous (e.g. across t=0t=0 in Example 1)
Two solutions found by different methodsCompute W[y1,y2](t0)W[y_1,y_2](t_0) at one convenient pointAssuming "two different-looking formulas" automatically means independent (Example 2 shows this fails)
Repeated characteristic rootsSecond solution is tertte^{rt}, not another multiple of erte^{rt}Writing c1ert+c2ertc_1e^{rt}+c_2e^{rt} as if it were a 2-parameter family
Interval of validityState II explicitly; theorem is local to an interval where p,qp,q are continuousClaiming the general solution formula holds across a singularity of pp or qq
Higher order n>2n>2Need exactly nn independent solutions, Wronskian is an n×nn\times n determinantStopping after finding 2 solutions for a 3rd- or 4th-order equation

The mathematician behind it

Origins · 1 minJózef Hoëné-Wroński and how this was discovered
JH
Józef Hoëné-Wroński
1776–1853 · Polish mathematician and philosopher working in post-Napoleonic France

Wroński introduced the determinant that now bears his name while developing an ambitious, often obscure philosophical program aiming to unify mathematics, physics, and metaphysics through what he called "absolute" laws of creation. Much of his broader system was dismissed by contemporaries as impenetrable or unfounded, but the determinant construction he isolated for testing functional dependence survived and was absorbed into mainstream analysis, eventually becoming indispensable to the rigorous theory of linear differential equations developed over the following century.

History
- 1776: Józef Hoëné-Wroński born in Wolsztyn, Poland.
- 1812: Wroński introduces the determinant later named the Wronskian, in work on functional dependence.
- 1830s–1840s: The Wronskian is adopted into the developing theory of linear ordinary differential equations.
- 1890s–1900s: Rigorous existence–uniqueness theorems (Picard, Peano, Lindelöf) supply the missing engine that turns "nonvanishing Wronskian" into "spans all solutions."
- 1853: Wroński dies in Neuilly-sur-Seine, France.
What exactly is the difference between "a solution" and "the general solution"?

"A solution" is any single function satisfying the ODE (and possibly initial conditions); there can be infinitely many. "The general solution" is a formula with arbitrary constants, built from a fundamental set, that is guaranteed by the Fundamental Set Theorem to produce every solution on the interval as those constants range over all real numbers. The word "the" is doing real work: the theorem proves this formula is not just a large family but the entire solution set, with no solutions left outside it. Without the theorem, "general solution" would be an unjustified name for what might be an incomplete list.

Do I need to verify the Wronskian is nonzero at every point, or just one?

Just one point, and that is one of the more subtle payoffs of the theory. If y1,y2y_1,y_2 both solve the same linear homogeneous equation with continuous coefficients on II, Abel's formula (proved under The Wronskian and Linear Independence) shows W[y1,y2](t)=W[y1,y2](t0)exp⁡(−∫t0tp(s) ds)W[y_1,y_2](t) = W[y_1,y_2](t_0)\exp\left(-\int_{t_0}^t p(s)\,ds\right), and since the exponential factor is never zero, WW is either identically zero on II or never zero on II. So checking a single convenient point (often t0=0t_0=0) settles the matter for all t∈It\in I simultaneously.

Can a second-order linear homogeneous ODE have three linearly independent solutions?

No — this is a direct corollary of the Fundamental Set Theorem. If y1,y2y_1,y_2 form a fundamental set, every solution equals c1y1+c2y2c_1y_1+c_2y_2 for unique constants, so the solution space has dimension exactly 22. Any third solution y3y_3 must itself be expressible as c1y1+c2y2c_1y_1+c_2y_2 for some constants, which means y3,y1,y2y_3, y_1, y_2 are linearly dependent. This is why "the solution space is nn-dimensional" is not just descriptive language but a hard numerical ceiling enforced by the theorem.

What if the coefficients p(t)p(t) or q(t)q(t) are discontinuous somewhere in my interval of interest?

Then the hypotheses of the theorem fail at that point, and you must split the domain. Restrict to a subinterval on which p,qp,q are both continuous, apply the theorem there, and treat each side of a discontinuity or singularity separately (as in Example 1, where (0,∞)(0,\infty) was used specifically to avoid the singularity of p,qp,q at t=0t=0). Matching solutions across such a point, if required, needs additional continuity/smoothness conditions on the solution itself and is a separate analysis, not covered by this theorem directly.

Is the fundamental set theorem specific to homogeneous equations, or does it apply to y′′+py′+qy=g(t)y''+py'+qy=g(t) too?

As stated, it is specifically about the homogeneous equation (g≡0g\equiv0), because the proof's Step 3 crucially uses that a linear combination of solutions is again a solution — true only when g=0g=0. For the nonhomogeneous equation, the solution set is not a vector space (it doesn't contain 00 unless g≡0g\equiv 0), but it is still fully described using this theorem: the general solution of the nonhomogeneous equation is y=yp+c1y1+c2y2y = y_p + c_1y_1+c_2y_2, where ypy_p is any single particular solution and c1y1+c2y2c_1y_1+c_2y_2 is the general solution of the associated homogeneous equation supplied by this theorem.

How does this generalize to systems of first-order linear ODEs?

Any nnth-order linear ODE converts to a first-order linear system x′=A(t)x\mathbf{x}' = A(t)\mathbf{x} in Rn\mathbb{R}^n, and the exact same structure holds: a set of nn vector-valued solutions is a fundamental set exactly when the matrix Wronskian (determinant of the matrix whose columns are the solutions) is nonzero at one point, and then it is nonzero everywhere by an analogue of Abel's formula. The general solution is x(t)=c1x1(t)+⋯+cnxn(t)\mathbf{x}(t) = c_1\mathbf{x}_1(t)+\cdots+c_n\mathbf{x}_n(t), and the proof is structurally identical, replacing scalar Cramer's rule with the invertibility of the fundamental matrix. This is the natural bridge to state-space and matrix-exponential methods.

Why does linear independence of solutions have to be checked with the Wronskian rather than "by inspection"?

Because inspection is unreliable: two solutions can look completely different algebraically yet be proportional (as y2=2e2ty_2=2e^{2t} looked different from y1=e2ty_1=e^{2t} only by a constant, which is easy to spot, but subtler disguises — trigonometric identities, hidden algebraic relations — are not always obvious). The Wronskian converts the question "are these functions proportional?" into a concrete computation and, thanks to Abel's formula, a single-point check suffices for solutions of the same linear homogeneous ODE. This removes any guesswork from the independence check that the Fundamental Set Theorem requires.

Before this

Where this leads