First prove additivity of variation: for aβ€x1β<x2ββ€b, Vax2ββ(f)=Vax1ββ(f)+Vx1βx2ββ(f).
The inequality β€ goes as follows: given any partition P of [a,x2β], if x1β is absent we may insert it, since the sum can only grow (β£f(v)βf(u)β£β€β£f(x1β)βf(u)β£+β£f(v)βf(x1β)β£).
A partition containing x1β has its sum split into a sum over [a,x1β] plus a sum over [x1β,x2β], each bounded by the variation on that interval; taking the supremum over P gives β€.
The reverse β₯: concatenating a partition P1β of [a,x1β] with a partition P2β of [x1β,x2β] gives a partition P1ββͺP2β of [a,x2β] whose sum is the sum of the two sums.
Taking suprema on each interval gives β₯.
All suprema are finite since f is BV.
Now g(x)=Vaxβ(f) is increasing: for x1β<x2β, additivity gives g(x2β)=g(x1β)+Vx1βx2ββ(f)β₯g(x1β) (variation is nonnegative).
And h(x2β)βh(x1β)=(g(x2β)βg(x1β))β(f(x2β)βf(x1β))=Vx1βx2ββ(f)β(f(x2β)βf(x1β)); but Vx1βx2ββ(f)β₯β£f(x2β)βf(x1β)β£β₯f(x2β)βf(x1β) from the single partition {x1β,x2β}, so h(x2β)βh(x1β)β₯0 β h is increasing too.
The identity f=gβh is the definition of h.