D'Alembert's formula (source code)

= D'Alembert's formula
{c}
{title2=$u(x,t)=\tfrac12[f(x-ct)+f(x+ct)]+\tfrac1{2c}\int_{x-ct}^{x+ct}g(s)\,ds$}
{wiki}

= D'Alembert formula
{c}
{synonym}

For constant <wave speed> $c>0$, the whole-line <Cauchy problem>
$$
u_{tt}=c^2u_{xx},\qquad u(x,0)=f(x),\qquad u_t(x,0)=g(x)
$$
has the unique <classical solution>, for $f\in C^2(\mathbb R)$ and $g\in C^1(\mathbb R)$,
$$
u(x,t)=\frac{f(x-ct)+f(x+ct)}2+\frac1{2c}\int_{x-ct}^{x+ct}g(s)\,ds.
$$
The first term evolves displacement released from rest; the second evolves initial <velocity>. The factor $1/(2c)$ converts an integral of <velocity> over distance into displacement. On a finite string the formula applies before the backward <characteristic curves> reach an endpoint, or globally after an extension that enforces the <boundary conditions>. <Math 220A: Partial Differential Equations of Applied Mathematics/5 Wave Equation in R>{ref}

For a completely smooth example, choose a <bump function> with <compact support>:
$$
b(z)=\begin{cases}\exp\!\left(1-\frac1{1-z^2}\right),&|z|<1,\\0,&|z|\geq1,\end{cases}
\qquad f(x)=A b(x/a),\qquad g(x)=0,
$$
where $a>0$ sets the half-width and $A$ sets the initial peak. Then
$$
u(x,t)=\frac A2 b\!\left(\frac{x-ct}{a}\right)+\frac A2 b\!\left(\frac{x+ct}{a}\right).
$$
At $t=0$ the two copies coincide, recovering peak $A$. At $t=a/(2c)$, $u(0,t)=A e^{-1/3}\simeq0.716531A$. At $t=3a/(2c)$ the pulses are separated, centred at $x=\pm3a/2$, each of width $2a$ and peak $A/2$. Zero initial <velocity> forces equal left and right contributions; a purely right-moving initial profile instead requires $g=-cf'$. Thus “released from rest” and “launched to the right” describe different <Cauchy data>.

\Image[/wave-pulse-propagation.png]
{title=A smooth pulse released from rest splits in two, while backward characteristics select a causal interval}

The <domain of dependence> of $(x_*,t_*)$, for $t_*\geq0$, is contained in $[x_*-ct_*,x_*+ct_*]$ on the initial line. More precisely, displacement contributes only through the two endpoints, while initial <velocity> contributes throughout the interval. With $c=1$, $x_*=0.5$ and $t_*=1.5$, that interval is $[-1,2]$. Changing either initial datum outside it cannot change $u(0.5,1.5)$. If both initial data vanish outside $[-a,a]$, the solution vanishes outside $[-a-ct,a+ct]$: this is <finite propagation speed>. The expanding interval bounds possible influence; it need not be filled by nonzero displacement, as the separated bumps demonstrate.

Historically, d'Alembert studied the vibrating-string problem in “Recherches sur la courbe que forme une corde tendue mise en vibration”, in the Berlin Academy's volume for the year 1747, pages 214–219. The memoir introduces the differential equation and integrates it through two arbitrary one-variable functions. The original geometric problem makes the formula's enduring appeal clear: arbitrary initial shapes can be handled without first decomposing them into <normal modes>. https://webusers.imj-prg.fr/~david.aubin/cours/Textes/Dalembert-HAB-1747-cordes-vibrantes.pdf[D'Alembert, vibrating-string memoir, Academy year 1747]{ref}

For a brief derivation, use <characteristic coordinates> $\xi=x-ct$ and $\eta=x+ct$. The <chain rule> gives $u_{tt}-c^2u_{xx}=-4c^2u_{\xi\eta}$, so $u=F(\xi)+G(\eta)$. The initial conditions yield $F+G=f$ and $-cF'+cG'=g$. Integrating these relations gives the displayed formula; the arbitrary integration constants cancel in their sum.