Take the partition PNβ={0,vNβ,uNβ,β¦,v1β,u1β,1}.
The contributions from [0,vNβ] and [u1β,1] are nonnegative, so dropping them only lowers the sum; on each pair (vkβ,ukβ) we get β£f(ukβ)βf(vkβ)β£=β£ukββ(βvkβ)β£=ukβ+vkβ.
Hence VPNββ(f)β₯βk=1Nβ(ukβ+vkβ).
Since the partition variations diverge, supPβVPβ(f)=β: f is not of bounded variation on [0,1]. Numerical check: u1β=2.5Ο1ββ0.12732, v1β=3.5Ο1ββ0.09095, βk=1200β(ukβ+vkβ)β1.6849>1 β the scene's V200ββ2.1388 already exceeds this lower bound. (Contrast: x2sin(1/x) IS of bounded variation, with Vβ€3, since β«01ββ£fβ²β£β€2.)