Crossposted from [MSE.](https://math.stackexchange.com/questions/4721350/changing-variables-in-discrete-calculus)


In discrete calculus one defines the $h$-difference operator $$\Delta\_h[f(x)] = f(x+h) - f(x)$$
and we often define $\Delta = \Delta\_1.$ We can similarly define the [indefinite sums](https://en.wikipedia.org/wiki/Indefinite_sum) $\Delta\_h^{-1}$ and set $\Delta^{-1} = \Delta\_1^{-1}.$


Most books on discrete calculus only include results for $h=1.$
For example, the above link gives that
$$\Delta^{-1}\sin rx = \frac{-\cos(r[x-\frac{1}{2}])}{2\sin\frac{r}{2}}$$


Through an ugly computation I managed to work out that $$\Delta\_h^{-1}[\sin rx] = \frac{-\cos(r[x-\frac{h}{2}])}{2\sin\frac{rh}{2}}$$


What I'd like is a uniform procedure for recovering the "$h$-version" from such results. How could I have derived (or even just guessed) the second formula from the first? Is there something like the change of variable formula in calculus? Maybe some version of the 'dimensional analysis' heuristics in ordinary calculus that tells you where to insert the $h$'s? As a test, given a formula at random from the previous link like $\Delta[\log x] = \log(1 + \frac{1}{x})$ I'd like to be able to see immediately that $\Delta\_h[\log x] = \log(1 + \frac{h}{x})$ without carrying out the actual calculation.


A search trough the standard references didn't turn up much, so any help or references would be welcome.



Change of variable formulas in discrete calculus?
In terms of the differential operator $\partial\_x\equiv d/dx$ one has $f(x+h)=e^{h\partial\_x}f(x)$, hence
$$\Delta\_h=e^{h\partial\_x}-1.$$
Upon Fourier transformation, $\hat{f}(k)=\int\_{-\infty}^\infty e^{ikx}f(x)\,dx$, one has $\partial\_x\mapsto -ik$, hence
$$\hat{\Delta}\_h=e^{-ihk}-1\Leftrightarrow \hat{\Delta}\_h^{-1}=\frac{e^{ihk}}{1-e^{ihk}}.$$


In this way you can calculate, for any function $f(x)$ with Fourier transform $\hat{f}(k)$:
$$\Delta\_h^{-1}f(x)=\int\_{-\infty}^\infty\frac{e^{i(h-x)k}}{1-e^{ihk}}\hat{f}(k)\,\frac{dk}{2\pi}=\int\_{-\infty}^\infty F\_h(x'-x)f(x')\,dx'$$
with [Dirac-delta-function](https://www.google.com/search?client=safari&rls=en&q=dirac+delta+function&ie=UTF-8&oe=UTF-8) kernel
$$F\_h(x)=\int\_{-\infty}^\infty\frac{e^{i(h+x)k}}{1-e^{ihk}}\,\frac{dk}{2\pi}=\sum\_{n=1}^\infty\delta(x+nh).$$
The relation $F\_h(x)=(1/h)F\_1(x/h)$ gives the scaling with $h$ requested by the OP:



> 
> $\Delta\_h^{-1}f(x)$ follows from $\Delta\_1^{-1}f(rx)$ by rescaling
> $r\mapsto rh$ and $x\mapsto x/h$.
> 
> 
> 




---


The "test cases" mentioned by the OP include: $\Delta\_1^{-1}\log rx=x\log r+\log \Gamma(x)$, hence
$$\Delta\_h^{-1}\log rx=(x/h)\log(rh)+\log\Gamma(x/h).$$
Similarly, $\Delta\_1^{-1}a^{rx}=a^{rx}/(a^r-1)$, hence
$$\Delta\_h^{-1} a^{rx}=a^{rx}/(a^{rh}-1),$$
and $\Delta\_1^{-1}(rx)^a=r^aB\_{a+1}(x)/(a+1)$, hence
$$\Delta\_h^{-1}(rx)^a=(rh)^aB\_{a+1}(x/h)/(a+1).$$


