Finite difference
A finite difference is a mathematical expression of the form f (x + b) − f (x + a). If a finite difference is divided by b − a, one gets a difference quotient. The approximation of derivatives by finite differences plays a central role in finite difference methods for the numerical solution of differential equations, especially boundary value problems.
Certain recurrence relations can be written as difference equations by replacing iteration notation with finite differences.
Today, the term "finite difference" is often taken as synonymous with finite difference approximations of derivatives, especially in the context of numerical methods.[1][2][3] Finite difference approximations are finite difference quotients in the terminology employed above.
Finite differences have also been the topic of study as abstract self-standing mathematical objects, such as in works by George Boole (1860), L. M. Milne-Thomson (1933), and Károly Jordan (1939), tracing its origins back to one of Jost Bürgi's algorithms (c. 1592) and others including Isaac Newton. In this viewpoint, the formal calculus of finite differences is an alternative to the calculus of infinitesimals.[4]
Contents
1 Forward, backward, and central differences
2 Relation with derivatives
3 Higher-order differences
3.1 Arbitrarily sized kernels
3.2 Properties
4 Finite difference methods
5 Newton's series
6 Calculus of finite differences
7 Rules for calculus of finite difference operators
8 Generalizations
9 Finite difference in several variables
10 See also
11 References
12 External links
Forward, backward, and central differences
Three forms are commonly considered: forward, backward, and central differences.[1][2][3]
A forward difference is an expression of the form
- Δh[f](x)=f(x+h)−f(x).{displaystyle Delta _{h}[f](x)=f(x+h)-f(x).}
Depending on the application, the spacing h may be variable or constant. When omitted, h is taken to be 1: Δ[ f ](x) = Δ1[ f ](x).
A backward difference uses the function values at x and x − h, instead of the values at x + h and x:
- ∇h[f](x)=f(x)−f(x−h).{displaystyle nabla _{h}[f](x)=f(x)-f(x-h).}
Finally, the central difference is given by
- δh[f](x)=f(x+12h)−f(x−12h).{displaystyle delta _{h}[f](x)=fleft(x+{tfrac {1}{2}}hright)-fleft(x-{tfrac {1}{2}}hright).}
Relation with derivatives
The derivative of a function f at a point x is defined by the limit.
- f′(x)=limh→0f(x+h)−f(x)h.{displaystyle f'(x)=lim _{hto 0}{frac {f(x+h)-f(x)}{h}}.}
If h has a fixed (non-zero) value instead of approaching zero, then the right-hand side of the above equation would be written
- f(x+h)−f(x)h=Δh[f](x)h.{displaystyle {frac {f(x+h)-f(x)}{h}}={frac {Delta _{h}[f](x)}{h}}.}
Hence, the forward difference divided by h approximates the derivative when h is small. The error in this approximation can be derived from Taylor's theorem. Assuming that f is differentiable, we have
- Δh[f](x)h−f′(x)=O(h)→0as h→0.{displaystyle {frac {Delta _{h}[f](x)}{h}}-f'(x)=O(h)to 0quad {text{as }}hto 0.}
The same formula holds for the backward difference:
- ∇h[f](x)h−f′(x)=O(h)→0as h→0.{displaystyle {frac {nabla _{h}[f](x)}{h}}-f'(x)=O(h)to 0quad {text{as }}hto 0.}
However, the central (also called centered) difference yields a more accurate approximation. If f is twice differentiable,
- δh[f](x)h−f′(x)=O(h2).{displaystyle {frac {delta _{h}[f](x)}{h}}-f'(x)=Oleft(h^{2}right).}
The main problem[citation needed] with the central difference method, however, is that oscillating functions can yield zero derivative. If f (nh) = 1 for n odd, and f (nh) = 2 for n even, then f ′(nh) = 0 if it is calculated with the central difference scheme. This is particularly troublesome if the domain of f is discrete.
Authors for whom finite differences mean finite difference approximations define the forward/backward/central differences as the quotients given in this section (instead of employing the definitions given in the previous section).[1][2][3]
Higher-order differences
In an analogous way, one can obtain finite difference approximations to higher order derivatives and differential operators. For example, by using the above central difference formula for f ′(x + h/2) and f ′(x − h/2) and applying a central difference formula for the derivative of f ′ at x, we obtain the central difference approximation of the second derivative of f:
- Second-order central
- f″(x)≈δh2[f](x)h2=f(x+h)−f(x)h−f(x)−f(x−h)hh=f(x+h)−2f(x)+f(x−h)h2.{displaystyle f''(x)approx {frac {delta _{h}^{2}[f](x)}{h^{2}}}={frac {{frac {f(x+h)-f(x)}{h}}-{frac {f(x)-f(x-h)}{h}}}{h}}={frac {f(x+h)-2f(x)+f(x-h)}{h^{2}}}.}
Similarly we can apply other differencing formulas in a recursive manner.
- Second order forward
- f″(x)≈Δh2[f](x)h2=f(x+2h)−f(x+h)h−f(x+h)−f(x)hh=f(x+2h)−2f(x+h)+f(x)h2.{displaystyle f''(x)approx {frac {Delta _{h}^{2}[f](x)}{h^{2}}}={frac {{frac {f(x+2h)-f(x+h)}{h}}-{frac {f(x+h)-f(x)}{h}}}{h}}={frac {f(x+2h)-2f(x+h)+f(x)}{h^{2}}}.}
- Second order backward
- f″(x)≈∇h2[f](x)h2=f(x)−f(x−h)h−f(x−h)−f(x−2h)hh=f(x)−2f(x−h)+f(x−2h)h2.{displaystyle f''(x)approx {frac {nabla _{h}^{2}[f](x)}{h^{2}}}={frac {{frac {f(x)-f(x-h)}{h}}-{frac {f(x-h)-f(x-2h)}{h}}}{h}}={frac {f(x)-2f(x-h)+f(x-2h)}{h^{2}}}.}
More generally, the nth order forward, backward, and central differences are given by, respectively,
- Forward
- Δhn[f](x)=∑i=0n(−1)i(ni)f(x+(n−i)h),{displaystyle Delta _{h}^{n}[f](x)=sum _{i=0}^{n}(-1)^{i}{binom {n}{i}}f{bigl (}x+(n-i)h{bigr )},}
or for h = 1,
- Δn[f](x)=∑k=0n(nk)(−1)n−kf(x+k){displaystyle Delta ^{n}[f](x)=sum _{k=0}^{n}{binom {n}{k}}(-1)^{n-k}f(x+k)}
- Backward
- ∇hn[f](x)=∑i=0n(−1)i(ni)f(x−ih),{displaystyle nabla _{h}^{n}[f](x)=sum _{i=0}^{n}(-1)^{i}{binom {n}{i}}f(x-ih),}
- Central
- δhn[f](x)=∑i=0n(−1)i(ni)f(x+(n2−i)h).{displaystyle delta _{h}^{n}[f](x)=sum _{i=0}^{n}(-1)^{i}{binom {n}{i}}fleft(x+left({frac {n}{2}}-iright)hright).}
These equations use binomial coefficients after the summation sign shown as (n
i). Each row of Pascal's triangle provides the coefficient for each value of i.
Note that the central difference will, for odd n, have h multiplied by non-integers. This is often a problem because it amounts to changing the interval of discretization. The problem may be remedied taking the average of δn[ f ](x − h/2) and δn[ f ](x + h/2).
Forward differences applied to a sequence are sometimes called the binomial transform of the sequence, and have a number of interesting combinatorial properties. Forward differences may be evaluated using the Nörlund–Rice integral. The integral representation for these types of series is interesting, because the integral can often be evaluated using asymptotic expansion or saddle-point techniques; by contrast, the forward difference series can be extremely hard to evaluate numerically, because the binomial coefficients grow rapidly for large n.
The relationship of these higher-order differences with the respective derivatives is straightforward,
- dnfdxn(x)=Δhn[f](x)hn+O(h)=∇hn[f](x)hn+O(h)=δhn[f](x)hn+O(h2).{displaystyle {frac {d^{n}f}{dx^{n}}}(x)={frac {Delta _{h}^{n}[f](x)}{h^{n}}}+O(h)={frac {nabla _{h}^{n}[f](x)}{h^{n}}}+O(h)={frac {delta _{h}^{n}[f](x)}{h^{n}}}+Oleft(h^{2}right).}
Higher-order differences can also be used to construct better approximations. As mentioned above, the first-order difference approximates the first-order derivative up to a term of order h. However, the combination
- Δh[f](x)−12Δh2[f](x)h=−f(x+2h)−4f(x+h)+3f(x)2h{displaystyle {frac {Delta _{h}[f](x)-{frac {1}{2}}Delta _{h}^{2}[f](x)}{h}}=-{frac {f(x+2h)-4f(x+h)+3f(x)}{2h}}}
approximates f ′(x) up to a term of order h2. This can be proven by expanding the above expression in Taylor series, or by using the calculus of finite differences, explained below.
If necessary, the finite difference can be centered about any point by mixing forward, backward, and central differences.
Arbitrarily sized kernels
Using linear algebra one can construct finite difference approximations which utilize an arbitrary number of points to the left and a (possibly different) number of points to the right of the evaluation point, for any order derivative. This involves solving a linear system such that the Taylor expansion of the sum of those points around the evaluation point best approximates the Taylor expansion of the desired derivative. Such formulas can be represented graphically on a hexagonal or diamond-shaped grid.[5]
This is useful for differentiating a function on a grid, where, as one approaches the edge of the grid, one must sample fewer and fewer points on one side.
The details are outlined in these notes.
The Finite Difference Coefficients Calculator constructs finite difference approximations for non-standard (and even non-integer) stencils given an arbitrary stencil and a desired derivative order.
Properties
- For all positive k and n
- Δkhn(f,x)=∑i1=0k−1∑i2=0k−1⋯∑in=0k−1Δhn(f,x+i1h+i2h+⋯+inh).{displaystyle Delta _{kh}^{n}(f,x)=sum limits _{i_{1}=0}^{k-1}sum limits _{i_{2}=0}^{k-1}cdots sum limits _{i_{n}=0}^{k-1}Delta _{h}^{n}left(f,x+i_{1}h+i_{2}h+cdots +i_{n}hright).}
Leibniz rule:
- Δhn(fg,x)=∑k=0n(nk)Δhk(f,x)Δhn−k(g,x+kh).{displaystyle Delta _{h}^{n}(fg,x)=sum limits _{k=0}^{n}{binom {n}{k}}Delta _{h}^{k}(f,x)Delta _{h}^{n-k}(g,x+kh).}
Finite difference methods
An important application of finite differences is in numerical analysis, especially in numerical differential equations, which aim at the numerical solution of ordinary and partial differential equations respectively. The idea is to replace the derivatives appearing in the differential equation by finite differences that approximate them. The resulting methods are called finite difference methods.
Common applications of the finite difference method are in computational science and engineering disciplines, such as thermal engineering, fluid mechanics, etc.
An open Python package of the finite difference method for arbitrary accuracy and order in any dimension on uniform and non-uniform grids is the Findiff project.
Newton's series
The Newton series consists of the terms of the Newton forward difference equation, named after Isaac Newton; in essence, it is the Newton interpolation formula, first published in his Principia Mathematica in 1687,[6] namely the discrete analog of the continuum Taylor expansion,
f(x)=∑k=0∞Δk[f](a)k!(x−a)k=∑k=0∞(x−ak)Δk[f](a),{displaystyle f(x)=sum _{k=0}^{infty }{frac {Delta ^{k}[f](a)}{k!}},(x-a)_{k}=sum _{k=0}^{infty }{binom {x-a}{k}},Delta ^{k}[f](a),}
which holds for any polynomial function f and for most (but not all) analytic functions. Here, the expression
- (xk)=(x)kk!{displaystyle {binom {x}{k}}={frac {(x)_{k}}{k!}}}
is the binomial coefficient, and
- (x)k=x(x−1)(x−2)⋯(x−k+1){displaystyle (x)_{k}=x(x-1)(x-2)cdots (x-k+1)}
is the "falling factorial" or "lower factorial", while the empty product (x)0 is defined to be 1. In this particular case, there is an assumption of unit steps for the changes in the values of x, h = 1 of the generalization below.
Note the formal correspondence of this result to Taylor's theorem. Historically, this, as well as the Chu–Vandermonde identity,
- (x+y)n=∑k=0n(nk)(x)n−k(y)k,{displaystyle (x+y)_{n}=sum _{k=0}^{n}{binom {n}{k}}(x)_{n-k},(y)_{k},}
(following from it, and corresponding to the binomial theorem), are included in the observations that matured to the system of umbral calculus.
To illustrate how one may use Newton's formula in actual practice, consider the first few terms of doubling the Fibonacci sequence f = 2, 2, 4, ... One can find a polynomial that reproduces these values, by first computing a difference table, and then substituting the differences that correspond to x0 (underlined) into the formula as follows,
- xf=Δ0Δ1Δ212_0_222_234f(x)=Δ0⋅1+Δ1⋅(x−x0)11!+Δ2⋅(x−x0)22!(x0=1)=2⋅1+0⋅x−11+2⋅(x−1)(x−2)2=2+(x−1)(x−2){displaystyle {begin{matrix}{begin{array}{|c||c|c|c|}hline x&f=Delta ^{0}&Delta ^{1}&Delta ^{2}\hline 1&{underline {2}}&&\&&{underline {0}}&\2&2&&{underline {2}}\&&2&\3&4&&\hline end{array}}&quad {begin{aligned}f(x)&=Delta ^{0}cdot 1+Delta ^{1}cdot {dfrac {(x-x_{0})_{1}}{1!}}+Delta ^{2}cdot {dfrac {(x-x_{0})_{2}}{2!}}quad (x_{0}=1)\\&=2cdot 1+0cdot {dfrac {x-1}{1}}+2cdot {dfrac {(x-1)(x-2)}{2}}\\&=2+(x-1)(x-2)\end{aligned}}end{matrix}}}
For the case of nonuniform steps in the values of x, Newton computes the divided differences,
- Δj,0=yj,Δj,k=Δj+1,k−1−Δj,k−1xj+k−xj∋{k>0,j≤max(j)−k},Δ0k=Δ0,k{displaystyle Delta _{j,0}=y_{j},qquad Delta _{j,k}={frac {Delta _{j+1,k-1}-Delta _{j,k-1}}{x_{j+k}-x_{j}}}quad ni quad left{k>0,;jleq max left(jright)-kright},qquad Delta 0_{k}=Delta _{0,k}}
the series of products,
- P0=1,Pk+1=Pk⋅(ξ−xk),{displaystyle {P_{0}}=1,quad quad P_{k+1}=P_{k}cdot left(xi -x_{k}right),}
and the resulting polynomial is the scalar product,[7]
f(ξ)=Δ0⋅P(ξ){displaystyle f(xi )=Delta 0cdot Pleft(xi right)} .
In analysis with p-adic numbers, Mahler's theorem states that the assumption that f is a polynomial function can be weakened all the way to the assumption that f is merely continuous.
Carlson's theorem provides necessary and sufficient conditions for a Newton series to be unique, if it exists. However, a Newton series does not, in general, exist.
The Newton series, together with the Stirling series and the Selberg series, is a special case of the general difference series, all of which are defined in terms of suitably scaled forward differences.
In a compressed and slightly more general form and equidistant nodes the formula reads
- f(x)=∑k=0(x−ahk)∑j=0k(−1)k−j(kj)f(a+jh).{displaystyle f(x)=sum _{k=0}{binom {frac {x-a}{h}}{k}}sum _{j=0}^{k}(-1)^{k-j}{binom {k}{j}}f(a+jh).}
Calculus of finite differences
The forward difference can be considered as an operator, called the difference operator, which maps the function f to Δh[ f ].[8][9] This operator amounts to
- Δh=Th−I,{displaystyle Delta _{h}=T_{h}-I,}
where Th is the shift operator with step h, defined by Th[ f ](x) = f (x + h), and I is the identity operator.
The finite difference of higher orders can be defined in recursive manner as Δn
h ≡ Δh(Δn − 1
h). Another equivalent definition is Δn
h = [Th − I]n.
The difference operator Δh is a linear operator, as such it satisfies Δh[αf + βg](x) = α Δh[ f ](x) + β Δh[g](x).
It also satisfies a special Leibniz rule indicated above,
Δh(f (x)g(x)) = (Δhf (x)) g(x+h) + f (x) (Δhg(x)). Similar statements hold for the backward and central differences.
Formally applying the Taylor series with respect to h, yields the formula
- Δh=hD+12!h2D2+13!h3D3+⋯=ehD−I,{displaystyle Delta _{h}=hD+{frac {1}{2!}}h^{2}D^{2}+{frac {1}{3!}}h^{3}D^{3}+cdots =mathrm {e} ^{hD}-I,}
where D denotes the continuum derivative operator, mapping f to its derivative f ′. The expansion is valid when both sides act on analytic functions, for sufficiently small h. Thus, Th = ehD, and formally inverting the exponential yields
- hD=log(1+Δh)=Δh−12Δh2+13Δh3+⋯.{displaystyle hD=log(1+Delta _{h})=Delta _{h}-{tfrac {1}{2}}Delta _{h}^{2}+{tfrac {1}{3}}Delta _{h}^{3}+cdots .}
This formula holds in the sense that both operators give the same result when applied to a polynomial.
Even for analytic functions, the series on the right is not guaranteed to converge; it may be an asymptotic series. However, it can be used to obtain more accurate approximations for the derivative. For instance, retaining the first two terms of the series yields the second-order approximation to f ′(x) mentioned at the end of the section Higher-order differences.
The analogous formulas for the backward and central difference operators are
- hD=−log(1−∇h)andhD=2arsinh(12δh).{displaystyle hD=-log(1-nabla _{h})quad {text{and}}quad hD=2operatorname {arsinh} left({tfrac {1}{2}}delta _{h}right).}
The calculus of finite differences is related to the umbral calculus of combinatorics. This remarkably systematic correspondence is due to the identity of the commutators of the umbral quantities to their continuum analogs (h → 0 limits),
[Δhh,xTh−1]=[D,x]=I.{displaystyle left[{frac {Delta _{h}}{h}},x,T_{h}^{-1}right]=[D,x]=I.}
A large number of formal differential relations of standard calculus involving
functions f (x) thus map systematically to umbral finite-difference analogs involving f (xT−1
h).
For instance, the umbral analog of a monomial xn is a generalization of the above falling factorial (Pochhammer k-symbol),
- (x)n≡(xTh−1)n=x(x−h)(x−2h)⋯(x−(n−1)h),{displaystyle ~(x)_{n}equiv left(xT_{h}^{-1}right)^{n}=x(x-h)(x-2h)cdots {bigl (}x-(n-1)h{bigr )},}
so that
- Δhh(x)n=n(x)n−1,{displaystyle {frac {Delta _{h}}{h}}(x)_{n}=n(x)_{n-1},}
hence the above Newton interpolation formula (by matching coefficients in the expansion of an arbitrary function f (x) in such symbols), and so on.
For example, the umbral sine is
- sin(xTh−1)=x−(x)33!+(x)55!−(x)77!+⋯{displaystyle sin left(x,T_{h}^{-1}right)=x-{frac {(x)_{3}}{3!}}+{frac {(x)_{5}}{5!}}-{frac {(x)_{7}}{7!}}+cdots }
As in the continuum limit, the eigenfunction of Δh/h also happens to be an exponential,
- Δhh(1+λh)xh=Δhheln(1+λh)xh=λeln(1+λh)xh,{displaystyle {frac {Delta _{h}}{h}}(1+lambda h)^{frac {x}{h}}={frac {Delta _{h}}{h}}e^{ln(1+lambda h){frac {x}{h}}}=lambda e^{ln(1+lambda h){frac {x}{h}}},}
and hence Fourier sums of continuum functions are readily mapped to umbral Fourier sums faithfully, i.e., involving the same Fourier coefficients multiplying these umbral basis exponentials.[10] This umbral exponential thus amounts to the exponential generating function of the Pochhammer symbols.
Thus, for instance, the Dirac delta function maps to its umbral correspondent, the cardinal sine function,
- δ(x)↦sin[π2(1+xh)]π(x+h),{displaystyle delta (x)mapsto {frac {sin left[{frac {pi }{2}}left(1+{frac {x}{h}}right)right]}{pi (x+h)}},}
and so forth.[11]Difference equations can often be solved with techniques very similar to those for solving differential equations.
The inverse operator of the forward difference operator, so then the umbral integral, is the indefinite sum or antidifference operator.
Rules for calculus of finite difference operators
Analogous to rules for finding the derivative, we have:
Constant rule: If c is a constant, then
- Δc=0{displaystyle Delta c=0}
Linearity: if a and b are constants,
- Δ(af+bg)=aΔf+bΔg{displaystyle Delta (af+bg)=a,Delta f+b,Delta g}
All of the above rules apply equally well to any difference operator, including ∇ as to Δ.
Product rule:
- Δ(fg)=fΔg+gΔf+ΔfΔg∇(fg)=f∇g+g∇f−∇f∇g{displaystyle {begin{aligned}Delta (fg)&=f,Delta g+g,Delta f+Delta f,Delta g\nabla (fg)&=f,nabla g+g,nabla f-nabla f,nabla gend{aligned}}}
Quotient rule:
- ∇(fg)=1gdet[∇f∇gfg](det[g∇g11])−1{displaystyle nabla left({frac {f}{g}}right)={frac {1}{g}}det {begin{bmatrix}nabla f&nabla g\f&gend{bmatrix}}left(det {begin{bmatrix}g&nabla g\1&1end{bmatrix}}right)^{-1}}
- or
- ∇(fg)=g∇f−f∇gg⋅(g−∇g){displaystyle nabla left({frac {f}{g}}right)={frac {g,nabla f-f,nabla g}{gcdot (g-nabla g)}}}
Summation rules:
- ∑n=abΔf(n)=f(b+1)−f(a)∑n=ab∇f(n)=f(b)−f(a−1){displaystyle {begin{aligned}sum _{n=a}^{b}Delta f(n)&=f(b+1)-f(a)\sum _{n=a}^{b}nabla f(n)&=f(b)-f(a-1)end{aligned}}}
See references.[12][13][14][15]
Generalizations
- A generalized finite difference is usually defined as
- Δhμ[f](x)=∑k=0Nμkf(x+kh),{displaystyle Delta _{h}^{mu }[f](x)=sum _{k=0}^{N}mu _{k}f(x+kh),}
where μ = (μ0,… μN) is its coefficient vector. An infinite difference is a further generalization, where the finite sum above is replaced by an infinite series. Another way of generalization is making coefficients μk depend on point x: μk = μk(x), thus considering weighted finite difference. Also one may make the step h depend on point x: h = h(x). Such generalizations are useful for constructing different modulus of continuity.
- The generalized difference can be seen as the polynomial rings R[Th]. It leads to difference algebras.
- Difference operator generalizes to Möbius inversion over a partially ordered set.
- As a convolution operator: Via the formalism of incidence algebras, difference operators and other Möbius inversion can be represented by convolution with a function on the poset, called the Möbius function μ; for the difference operator, μ is the sequence (1, −1, 0, 0, 0, ...).
Finite difference in several variables
Finite differences can be considered in more than one variable. They are analogous to partial derivatives in several variables.
Some partial derivative approximations are:
- fx(x,y)≈f(x+h,y)−f(x−h,y)2hfy(x,y)≈f(x,y+k)−f(x,y−k)2kfxx(x,y)≈f(x+h,y)−2f(x,y)+f(x−h,y)h2fyy(x,y)≈f(x,y+k)−2f(x,y)+f(x,y−k)k2fxy(x,y)≈f(x+h,y+k)−f(x+h,y−k)−f(x−h,y+k)+f(x−h,y−k)4hk.{displaystyle {begin{aligned}f_{x}(x,y)&approx {frac {f(x+h,y)-f(x-h,y)}{2h}}\f_{y}(x,y)&approx {frac {f(x,y+k)-f(x,y-k)}{2k}}\f_{xx}(x,y)&approx {frac {f(x+h,y)-2f(x,y)+f(x-h,y)}{h^{2}}}\f_{yy}(x,y)&approx {frac {f(x,y+k)-2f(x,y)+f(x,y-k)}{k^{2}}}\f_{xy}(x,y)&approx {frac {f(x+h,y+k)-f(x+h,y-k)-f(x-h,y+k)+f(x-h,y-k)}{4hk}}.end{aligned}}}
Alternatively, for applications in which the computation of f is the most costly step, and both first and second derivatives must be computed, a more efficient formula for the last case is
- fxy(x,y)≈f(x+h,y+k)−f(x+h,y)−f(x,y+k)+2f(x,y)−f(x−h,y)−f(x,y−k)+f(x−h,y−k)2hk,{displaystyle f_{xy}(x,y)approx {frac {f(x+h,y+k)-f(x+h,y)-f(x,y+k)+2f(x,y)-f(x-h,y)-f(x,y-k)+f(x-h,y-k)}{2hk}},}
since the only values to compute that are not already needed for the previous four equations are f (x + h, y + k) and f (x − h, y − k).
See also
- Arc elasticity
- Carlson's theorem
- Central differencing scheme
- Divided differences
- Finite difference coefficients
- Finite difference method
- Finite volume method
- Five-point stencil
- Gilbreath's conjecture
- Lagrange polynomial
- Modulus of continuity
- Multiplicative calculus
- Newton polynomial
- Nörlund–Rice integral
- Numerical differentiation
- Sheffer sequence
- Summation by parts
- Table of Newtonian series
- Taylor series
- Time scale calculus
- Umbral calculus
- Upwind differencing scheme for convection
References
^ abc Paul Wilmott; Sam Howison; Jeff Dewynne (1995). The Mathematics of Financial Derivatives: A Student Introduction. Cambridge University Press. p. 137. ISBN 978-0-521-49789-3..mw-parser-output cite.citation{font-style:inherit}.mw-parser-output .citation q{quotes:"""""""'""'"}.mw-parser-output .citation .cs1-lock-free a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/6/65/Lock-green.svg/9px-Lock-green.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .citation .cs1-lock-limited a,.mw-parser-output .citation .cs1-lock-registration a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/d/d6/Lock-gray-alt-2.svg/9px-Lock-gray-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .citation .cs1-lock-subscription a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/a/aa/Lock-red-alt-2.svg/9px-Lock-red-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration{color:#555}.mw-parser-output .cs1-subscription span,.mw-parser-output .cs1-registration span{border-bottom:1px dotted;cursor:help}.mw-parser-output .cs1-ws-icon a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/4/4c/Wikisource-logo.svg/12px-Wikisource-logo.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output code.cs1-code{color:inherit;background:inherit;border:inherit;padding:inherit}.mw-parser-output .cs1-hidden-error{display:none;font-size:100%}.mw-parser-output .cs1-visible-error{font-size:100%}.mw-parser-output .cs1-maint{display:none;color:#33aa33;margin-left:0.3em}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration,.mw-parser-output .cs1-format{font-size:95%}.mw-parser-output .cs1-kern-left,.mw-parser-output .cs1-kern-wl-left{padding-left:0.2em}.mw-parser-output .cs1-kern-right,.mw-parser-output .cs1-kern-wl-right{padding-right:0.2em}
^ abc Peter Olver (2013). Introduction to Partial Differential Equations. Springer Science & Business Media. p. 182. ISBN 978-3-319-02099-0.
^ abc M Hanif Chaudhry (2007). Open-Channel Flow. Springer. p. 369. ISBN 978-0-387-68648-6.
^ Jordán, op. cit., p. 1 and Milne-Thomson, p. xxi.
Milne-Thomson, Louis Melville (2000): The Calculus of Finite Differences (Chelsea Pub Co, 2000)
ISBN 978-0821821077
^ Fraser, Duncan C. (1 January 1909). "On the Graphic Delineation of Interpolation Formulæ". Journal of the Institute of Actuaries. 43 (2): 235–241. Retrieved 17 April 2017.
^ Newton, Isaac, (1687). Principia, Book III, Lemma V, Case 1
^ Richtmeyer, D. and Morton, K.W., (1967). Difference Methods for Initial Value Problems, 2nd ed., Wiley, New York.
^ Boole, George, (1872). A Treatise On The Calculus of Finite Differences, 2nd ed., Macmillan and Company. On line. Also, [Dover edition 1960]
^ Jordan, Charles, (1939/1965). "Calculus of Finite Differences", Chelsea Publishing. On-line: [1]
^ Zachos, C. (2008). "Umbral Deformations on Discrete Space-Time". International Journal of Modern Physics A. 23 (13): 2005–2014. arXiv:0710.2306. Bibcode:2008IJMPA..23.2005Z. doi:10.1142/S0217751X08040548.
^ Curtright, T. L.; Zachos, C. K. (2013). "Umbral Vade Mecum". Frontiers in Physics. 1. arXiv:1304.0429. Bibcode:2013FrP.....1...15C. doi:10.3389/fphy.2013.00015.
^ Levy, H.; Lessman, F. (1992). Finite Difference Equations. Dover. ISBN 0-486-67260-3.
^ Ames, W. F., (1977). Numerical Methods for Partial Differential Equations, Section 1.6. Academic Press, New York.
ISBN 0-12-056760-1.
^ Hildebrand, F. B., (1968). Finite-Difference Equations and Simulations, Section 2.2, Prentice-Hall, Englewood Cliffs, New Jersey.
^ Flajolet, Philippe; Sedgewick, Robert (1995). "Mellin transforms and asymptotics: Finite differences and Rice's integrals" (PDF). Theoretical Computer Science. 144 (1–2): 101–124. doi:10.1016/0304-3975(94)00281-M.
- Richardson, C. H. (1954): An Introduction to the Calculus of Finite Differences (Van Nostrand (1954) online copy
- Mickens, R. E. (1991): Difference Equations: Theory and Applications (Chapman and Hall/CRC)
ISBN 978-0442001360
External links
Hazewinkel, Michiel, ed. (2001) [1994], "Finite-difference calculus", Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 978-1-55608-010-4
Table of useful finite difference formula generated using Mathematica- D. Gleich (2005), Finite Calculus: A Tutorial for Solving Nasty Sums
- Discrete Second Derivative from Unevenly Spaced Points