Existence and Uniqueness (Picard)

April 14, 2026

Problem

Show that the IVP y′ = √y with y(0) = 0 has multiple solutions: y ≡ 0 and y = t²/4. Explain why Picard's uniqueness theorem does not apply here.

Explanation

The promise and its fine print

For a first-order ODE y=f(t,y)y' = f(t, y), we usually expect exactly one solution through a given point (t0,y0)(t_0, y_0). This is the content of the Picard–Lindelöf theorem (sometimes called the existence and uniqueness theorem, or the Picard theorem):

If ff and f/y\partial f / \partial y are continuous in a rectangle around (t0,y0)(t_0, y_0), then the IVP y=f(t,y), y(t0)=y0y' = f(t, y),\ y(t_0) = y_0 has a unique solution on some interval containing t0t_0.

Two hypotheses:

  1. f(t,y)f(t, y) is continuous.
  2. f/y\partial f / \partial y is continuous (equivalently, ff is Lipschitz in yy).

Today's problem shows what happens when hypothesis 2 fails.

The given IVP

y=y,y(0)=0.y' = \sqrt{y}, \qquad y(0) = 0.

So f(t,y)=yf(t, y) = \sqrt{y}. Continuous on y0y \ge 0, so hypothesis 1 holds. But fy=12y,\frac{\partial f}{\partial y} = \frac{1}{2\sqrt{y}}, which blows up as y0+y \to 0^+. At the initial point y=0y = 0, f/y\partial f / \partial y is not even defined, let alone continuous. Hypothesis 2 fails — Picard says nothing, and indeed uniqueness fails.

Step-by-step — find multiple solutions

Solution A: y0y \equiv 0.

Check: y=0y' = 0 and y=0=0\sqrt{y} = \sqrt{0} = 0. ✓ The zero function is a perfectly valid solution.

Solution B: y=t2/4y = t^2 / 4 for t0t \ge 0.

Check: y=t/2y' = t/2 and y=t2/4=t/2\sqrt{y} = \sqrt{t^2/4} = t/2 (for t0t \ge 0). ✓

And y(0)=0y(0) = 0. ✓

So we've found two solutions to the same IVP.

Solution C (and infinitely many more): the "peel-off" family.

For any t0t^\star \ge 0, define yt(t)={0if 0tt(tt)2/4if t>ty_{t^\star}(t) = \begin{cases} 0 & \text{if } 0 \le t \le t^\star \\ (t - t^\star)^2 / 4 & \text{if } t > t^\star \end{cases}

This stays at zero until time tt^\star, then "peels off" and grows as a shifted parabola. Check that it is continuously differentiable at the peel-off point:

  • At t=tt = t^\star: y(t)=0y(t^\star) = 0 from both pieces ✓.
  • Derivative from the left: 00. From the right: y=(tt)/20y' = (t - t^\star)/2 \to 0 as tt+t \to {t^\star}^+. ✓
  • And y=yy' = \sqrt{y} on each piece.

So yty_{t^\star} is a valid solution for every t0t^\star \ge 0. Uncountably many solutions share the initial condition.

Geometric interpretation — why the peel-off is possible

The slope field of y=yy' = \sqrt{y} is very flat near the tt-axis (yy small ⇒ slope y\sqrt{y} small). A solution crawling along y=0y = 0 can sit there indefinitely — the "force pulling it off" is proportional to y\sqrt{y}, which is zero when y=0y = 0. Unlike a Lipschitz-bounded vector field, where a slight kick away from an equilibrium would grow exponentially and force uniqueness, here the slow \sqrt{\cdot} growth is too weak to pin the solution.

By contrast, y=yy' = y also has the equilibrium y=0y = 0. But Picard applies (yy=1\partial_y y = 1 is continuous), and the only solution through (0,0)(0, 0) is y0y \equiv 0. You can't peel off. Linear pull makes all the difference.

Where uniqueness failure matters in the real world

Physics models usually obey Picard, but when they don't, it's often a red flag:

  • Control systems: if the dynamics near a fixed point have y\sqrt{y}-like slow restoring force, the system's response to a perturbation is non-unique — a tell that your model is missing something (e.g. a higher-order term).
  • Hydrodynamics / dam-break problems: y=yy' = \sqrt{y}-type equations describe a self-similar flow; the non-uniqueness reflects genuine physical indeterminacy in idealised models (the actual system requires extra information to pick the "correct" solution — often a viscosity term).
  • Renormalisation of chaotic systems: some RG flows have non-unique trajectories that must be resolved by physical input.

So uniqueness failure isn't a math bug — it's a warning that the model is incomplete.

The Osgood criterion — a refinement

Picard's hypothesis is sufficient but not necessary. A weaker condition for uniqueness is Osgood's criterion: there exists a function ω\omega with 0+du/ω(u)=\int_{0^+} du / \omega(u) = \infty such that f(t,y1)f(t,y2)ω(y1y2).|f(t, y_1) - f(t, y_2)| \le \omega(|y_1 - y_2|).

For f=yf = \sqrt{y}, the modulus is y1y2ω(y1y2)|\sqrt{y_1} - \sqrt{y_2}| \le \omega(|y_1 - y_2|) with ω(u)=u\omega(u) = \sqrt{u}. But 0+du/u<\int_{0^+} du/\sqrt{u} < \infty, so Osgood doesn't save us either — and indeed uniqueness fails.

A theorem about the peel-off family

The most "eager" solution (starts growing at t=0t = 0) is y=t2/4y = t^2/4. The most "lazy" is y0y \equiv 0. Every solution is sandwiched: 0y(t)t2/40 \le y(t) \le t^2/4 for all t0t \ge 0. This is the Peano funnel — the maximum and minimum solutions form an envelope; every IVP solution lies inside.

Common mistakes

  • Thinking uniqueness always holds. It holds under mild hypotheses, and those hypotheses are satisfied for most textbook ODEs. But the counterexamples exist and they're important.
  • Assuming initial conditions pin the solution uniquely. For well-behaved ODEs yes, but not in general. Always check Picard.
  • Confusing existence with uniqueness. Peano's theorem guarantees existence from mere continuity of ff, with no condition on f/y\partial f / \partial y. But existence alone doesn't imply uniqueness.
  • Dropping the absolute value under \sqrt{\cdot}. y\sqrt{y} is only real for y0y \ge 0. For y<0y < 0 the original ODE isn't defined at all.

Why this matters pedagogically

It's the simplest, cleanest cautionary tale about the limits of ODE theory. Students often take existence and uniqueness for granted after seeing it invoked once; showing that y=yy' = \sqrt{y} has continuum many solutions through the origin makes the lesson stick.

Try it in the visualization

Plot both y0y \equiv 0 (flat along the t-axis) and y=t2/4y = t^2 / 4 (the rising parabola) through the origin. Add a slider for the "peel-off time" tt^\star from 0 to 3 and watch the hybrid solution yty_{t^\star} — flat until tt^\star, then parabolic. Overlay a visualization of f/y=1/(2y)\partial f/\partial y = 1/(2\sqrt{y}) blowing up to \infty along the tt-axis, which is why Picard fails.

Interactive Visualization

Parameters

1.00
3.00
3.00
5.00
8.00
Your turn

Got your own math or physics problem?

Turn any problem into an interactive visualization like this one — powered by AI, generated in seconds. Free to try, no credit card required.

Sign Up Free to Try It30 free visualizations every day
Existence and Uniqueness (Picard) | MathSpin