Picture the area under a curve between two fence posts. Now let something change: the curve may rise or fall, and the fence posts may move. The Leibniz integration formula tells us how fast the enclosed area changes, and it says there are exactly two sources of change.
The first source is the curve itself. If the curve moves up a little everywhere between the posts, the area grows by the sum of all those small rises.
The second source is the posts. If the right post moves to the right, a thin strip is added at the right edge. Its width is how far the post moved, and its height is the height of the curve at the post. The left post works the other way: moving it to the right cuts a strip off the region. That is why the left boundary term comes with a minus sign.
Adding the interior contribution and the two edge contributions gives the full rate of change. This is the reason the formula has one integral and two boundary terms.
Concrete
This tool gives a nice formula to the derivative of a one dimensional integral with dependencies to the varied quantity every where, with f,a and b having the right conditions and a=a(0), b=b(0),
In one dimension the formula holds, for example, when f and ∂f/∂ε are continuous on a rectangle [ε1,ε2]×[x1,x2] that contains all the points (ε,x) with x between a(ε) and b(ε), and a, b are continuously differentiable.
The general form concerns a k-form ωt that depends on a parameter t, integrated over a k-dimensional compact oriented submanifold Ωt=ϕt(Ω0) (with boundary) of a manifold M. Here ϕt is the flow of a possibly time-dependent vector field Xt, so that dtdϕt(p)=Xt(ϕt(p)). Pulling the integral back to the fixed domain Ω0 and using dtdϕt∗ωt=ϕt∗(∂tωt+LXtωt), we get
dtd∫Ωtωt=∫Ωt(∂tωt+LXtωt).
With Cartan's formula LXω=ιXdω+dιXω and Stokes' theorem this becomes
The one-dimensional formula is the case M=R, ωt=f(t,x)dx, with t in the role of ε and Xt any vector field whose flow moves the endpoints along a(t) and b(t). A 1-form on R is a top form, so dωt=0. The boundary of [a,b] consists of the point b with positive and the point a with negative orientation, and ιX(fdx)=fX with X(b)=b˙, X(a)=a˙, which reproduces f(t,b)b˙−f(t,a)a˙.
For ωt=fdV an n-form on an n-dimensional region of Rn we have LX(fdV)=∇⋅(fX)dV, and the general formula becomes the Reynolds transport theorem of continuum mechanics,
dtd∫ΩtfdV=∫Ωt∂t∂fdV+∮∂ΩtfX⋅dA.
Why is it interesting?
This formula is really nice in one dimensional Variational Calculus and with a(0)=a, b(0)=b, and for the particle mechanics Noether's Theorem. For integrals over arbitrary manifolds, it generalizes to the Lie derivative of a volume element.