 Derivation of Wien's displacement law

At fixed $T>0$, substitute $x=hc/(\lambda k_{B}T)>0$ into the Planck function:Here $h$ is the Planck constant, $k_B$ the Boltzmann constant, and $c$ the speed of light. Since $dx/d\lambda =-x/\lambda \ne 0$, an interior maximum in wavelength is also a stationary point in $x$. A logarithmic derivative givesThe required positive root is $x_{\lambda }=4.9651142317\dots$. The rearranged equation also has $x=0$, but that endpoint is not a peak. For $g(x)=5(1-e^{-x})-x$, $g^{\prime }(x)=5e^{-x}-1$ changes sign once, at $x=\log 5$. Thus $g$ first increases and then decreases to negative values, crossing zero exactly once for $x>0$. Since $x^5/(e^x-1)$ tends to zero at both endpoints, this stationary point is the unique global maximum. Solving the definition of $x$ for wavelength givesThe inverse-temperature dependence follows before any numerical root-finding. More generally $B_{\lambda }(T)=T^{5}F(\lambda T)$ for one fixed function $F$, so spectra at all temperatures have the same rescaled shape. This is stronger than a relation between two peak positions. Replacing the denominator by $e^x$ in the Wien approximation would instead give $x=5$ and a peak wavelength about $0.70\%$ too short.