# Differenzenquotient

Der Differenzenquotient ist ein Begriff aus der Mathematik. Er beschreibt das Verhältnis der Veränderung einer Größe zu der Veränderung einer anderen, wobei die erste Größe von der zweiten abhängt. In der Analysis verwendet man Differenzenquotienten, um die Ableitung einer Funktion zu definieren. In der numerischen Mathematik werden sie zum Lösen von Differentialgleichungen und für die näherungsweise Bestimmung der Ableitung einer Funktion (numerische Differentiation) benutzt.
Das gilt auch für Übertragungsfunktionen $G(s)$ der Systemtheorie, der Steuerungs- und Regelungstechnik für dynamische Systeme mit dem Ausgangs-Eingangsverhältnis der Laplace-transformierten gewöhnlichen Differenzialgleichungen (mit Störfunktion). Sie werden mit der inversen Laplace-Transformation auf gewöhnliche Differenzialgleichungen zurückgeführt und können mit Hilfe des Differenzenquotienten näherungsweise numerisch gelöst werden.

## Definition

Ist $f\colon D_{f}\to \mathbb {R}$ eine reellwertige Funktion, die im Bereich $D_{f}\subset \mathbb {R}$ definiert ist, und ist $[x_{0};x_{1}]\subset D_{f}$, so nennt man den Quotienten
    $\varphi (x_{1},x_{0})={\frac {f(x_{1})-f(x_{0})}{x_{1}-x_{0}}}$
Differenzenquotient von $f$ im Intervall $[x_{0};x_{1}]$.
Schreibt man $\Delta x:=x_{1}-x_{0}$ und $\Delta y:=f\left(x_{1}\right)-f\left(x_{0}\right)$, dann ergibt sich die alternative Schreibweise
    ${\frac {\Delta y}{\Delta x}}={\frac {f(x_{1})-f(x_{0})}{x_{1}-x_{0}}}$.
Setzt man $h=x_{1}-x_{0}$, also $x_{1}=x_{0}+h$, so erhält man die Schreibweise
    ${\frac {f(x_{0}+h)-f(x_{0})}{h}}$.
Geometrisch entspricht der Differenzenquotient der Steigung der Sekante des Graphen von $f$ durch die Punkte $(x_{0},f(x_{0}))$ und $(x_{1},f(x_{1}))$. Für $x_{1}\rightarrow x_{0}$ bzw. $h\rightarrow 0$ wird aus der Sekante eine Tangente an der Stelle $x_{0}$.

## Differentialquotient

Differenzenquotienten bilden zusammen mit dem Grenzwertbegriff eine Grundlage der Differentialrechnung. Den Grenzwert des Differenzenquotienten für $\displaystyle x_{1}\rightarrow x_{0}$ bezeichnet man als Differentialquotienten oder Ableitung der Funktion an der Stelle $x_{0}$, sofern dieser Grenzwert existiert.
Die Tabelle zeigt die Ableitungen einiger Funktionen. Dabei stimmt der Differenzenquotient jeweils nur für $x_{1}\neq x_{0}$.
| Funktion | $\displaystyle f(x)$ | Differenzenquotient ${\tfrac {f(x_{1})-f(x_{0})}{x_{1}-x_{0}}}$ | Differentialquotient $\lim _{x_{1}\rightarrow x_{0}}{\tfrac {f(x_{1})-f(x_{0})}{x_{1}-x_{0}}}$ |
| Konstante Funktion | $\displaystyle c$ | $\displaystyle 0$ | $\displaystyle 0$ |
| Lineare Funktion | $\displaystyle a\cdot x$ | $\displaystyle a$ | $\displaystyle a$ |
| Quadratfunktion | $\displaystyle x^{2}$ | $\displaystyle x_{1}+x_{0}$ | $\displaystyle 2\cdot x_{0}$ |
| Kubikfunktion | $\displaystyle x^{3}$ | $\displaystyle x_{1}^{2}+x_{1}\cdot x_{0}+x_{0}^{2}$ | $\displaystyle 3\cdot x_{0}^{2}$ |
| Allgemeine Potenz | $\displaystyle x^{n}$ | $\displaystyle \sum _{i=0}^{n-1}{x_{1}^{i}\cdot x_{0}^{n-1-i}}$ | $\displaystyle n\cdot x_{0}^{n-1}$ |
| Exponentialfunktion | $\displaystyle \exp(x)$ | $\displaystyle \exp(x_{0})\cdot {\frac {\exp(x_{1}-x_{0})-1}{x_{1}-x_{0}}}$ | $\displaystyle \exp(x_{0})$ |

## Definition der Varianten von Differenzenquotienten erster Ableitungsordnung

In der numerischen Mathematik werden zur Behandlung und Lösung von meist gewöhnlichen Differenzialgleichungen mit konstanten Koeffizienten die kontinuierlichen Funktionswerte mit Hilfe einer Differenzengleichung in Abhängigkeit von konstanten Intervallen ($\Delta x$, $\Delta t$, h) hintereinander berechnet. Die numerische Lösung einer gewöhnlichen Differenzengleichung erfolgt rekursiv über viele Berechnungsfolgen und stellt sich meist als eine tabellarisch geordnete Aufstellung von System-Ausgangsfolgen $y_{k}$ (Stützstellen, Knoten) in Abhängigkeit von der unabhängigen Variablen $x_{k}$, oder bei zeitabhängigen Systemen $t_{k}$ dar.

### Begriffsklärungen

* Gewöhnliche Differenzialgleichung mit konstanten Koeffizienten
    Differenzialgleichungen beschreiben dynamische Vorgänge unserer Umwelt wie Technik, Natur, Wirtschaft. Eine Differentialgleichung enthält außer der gesuchten Funktion z. B. $y(t)$ auch eine Ableitung der gesuchten Funktion $y'(t)$. Eine Differentialgleichung nennt man gewöhnlich, wenn die gesuchte Funktion nur von einer Veränderlichen (Variable) abhängt.
$\to$ siehe Gewöhnliche Differenzialgleichung#Begriffsdefinitionen der Differenzialgleichungen (DGL)
* Differenzengleichung
    Numerische Berechnungen von gewöhnlichen Differenzialgleichungen erfolgen über Differenzengleichungen. Dabei entsteht anstelle der kontinuierlichen Funktion eine endliche Anzahl von nummerierten Folgegliedern (Wertefolgen). Jede Folge bezieht sich bei einer Differenzengleichung 1. Ordnung rekursiv auf eine zurückliegende Folge.
    Für die numerische Lösung von Differenzialgleichungen existieren zahlreiche Varianten von Differenzengleichungen. Mit steigender Komplexität der Differenzengleichungen wird erreicht, dass für eine gleiche Approximation an den analytischen Verlauf der Ausgangsgröße $y(x)$ die Zahl der Folgeglieder $k_{\text{max}}$ erheblich reduziert wird.
* Mathematische Folge
    Bei der Lösung einer Differenzialgleichung mit Differenzengleichungen handelt es sich um eine rekursive Folge von nummerierten Elementen, also um eine Aufzählung von meist Funktionen oder auch von Zahlen. Je nach Art der gewöhnlichen Differenzialgleichung und der zugehörigen Differenzengleichung erhalten die Eingangs- und Ausgangsfolgeglieder der Differenzengleichung für die Nummerierung die Indizierung $k$.
       $k=[0,1,2,3,\dotsc ,k_{\text{max}}]$
    Arithmetische Folge:
    Der gegebene Folgewert wächst oder sinkt mit jedem Folgeglied um einen festen Betrag. Beispiel: „Sparschwein“.
    Exponentielle Folge:
    Der gegebene Folgewert wächst oder sinkt mit jedem Folgeglied um einen gleichen Prozentsatz bzw. um den gleichen relativen Anteil. Beispiel: Zinseszins.
* Indizierung der Folgen
    Die abhängige Variable $y_{k+1}$ entspricht dem nächsten beliebig nummerierten Folgeglied $y_{k}$ nach einem Rechenschritt $\Delta x=h$. Die abhängige Variable $y_{k-1}$ entspricht einem zurückliegenden beliebig nummeriertem Folgeglied $y_{k}$ vor einem Rechenschritt $\Delta x=h$.
* Differenzenverfahren
    Das Ergebnis der Lösung einer gewöhnlichen Differenzialgleichung ergibt eine kontinuierliche Funktion. Durch Überführung der Differenzialgleichung in eine Differenzengleichung ergibt die Lösung eine diskrete Funktion.
    Mit Hilfe des Differenzenverfahrens lassen sich häufig mit geringem Aufwand Differenzengleichungen aufstellen, die je nach der gewählten Schrittweite $\Delta x$ eine mehr oder weniger gute diskrete Annäherung an die analytische Lösung bilden.
    Gewöhnliche lineare Differenzialgleichungen, die z. B. ein dynamisches System 1. Ordnung beschreiben:
       $a_{1}\cdot y'(t)+a_{0}\cdot y(t)=b_{0}\cdot u(t)$
    können nach dem Differenzenverfahren relativ einfach in eine Differenzengleichung überführt werden. Dies geschieht dadurch, dass die Differenzialquotienten der Differenzialgleichung direkt durch die verschiedenen Formen der Differenzenquotienten ausgetauscht werden. Damit entsteht automatisch die rekursive Differenzengleichung.
    Das klassische Verfahren der Lösung von Differentialgleichungen mit Differenzengleichungen ist das eulersche Streckenzugverfahren zur numerischen Lösung eines Anfangswertproblems.

### Differenzenquotient

In der Numerik versteht man unter einem Differenzenquotienten die zeitdiskrete Form des Differentialquotienten einer gewöhnlichen Differentialgleichung. Die Differenzenquotienten mit der abhängigen Variablen $y(x)$ und der unabhängigen Variablen $x$ bzw. bei zeitabhängigen Funktionen der Variable $t$ unterscheiden folgende Methoden:
Der Vorwärts-Differenzenquotient für eine Funktion $y'=f(x,y)$ bezieht sich auf die linke Intervallgrenze laut Diagramm $x_{(k)}$ nach $x_{(k+1)}$ mit dem Intervall $h$.
    $y'\approx {\frac {\Delta y}{\Delta x}}={\frac {y(x+h)-y(x)}{(x+h)-x}}={\frac {y_{(k+1)}-y_{(k)}}{h}}$
Der Rückwärts-Differenzenquotient bezieht sich auf die rechte Intervallgrenze rückwärts nach dem Intervall $h$ von $x_{(k)}$ nach $x_{(k-1)}$.
    $y'\approx {\frac {\Delta y}{\Delta x}}={\frac {y(x)-y(x-h)}{x-(x-h)}}={\frac {y_{(k)}-y_{(k-1)}}{h}}$
Der zentrale Differenzenquotient bezieht sich auf die rechte und linke Intervallgrenze: ${y_{(k+1)}{\text{ nach }}y_{(k-1)}}{\text{ mit }}{2\cdot h}$
    $y'\approx {\frac {\Delta y}{\Delta x}}={\frac {y(x+h)-y(x-h)}{(x+h)-(x-h)}}={\frac {y_{(k+1)}-y_{(k-1)}}{2\cdot h}}$
Wird der zentrale Differenzenquotient in eine Differenzialgleichung eingesetzt, handelt es sich nicht um einen arithmetischen Mittelwert zweier Verfahren. Die hohe Genauigkeit der Annäherung an eine analytische Funktion steigt nicht mit fallendem Wert von $h$, sondern mit dem Quadrat des fallenden Wertes von $h$.

### Numerische Behandlung einer Differenzialgleichung mit Hilfe des Differenzenquotienten

* Folgende einfache Differenzialgleichung ist gegeben:
       $y'(x)=y(x)+e^{x}$.
    Die Lösung einer gewöhnlichen Differenzialgleichung 1. Ordnung ergibt in der Regel eine allgemeine Lösung in Form einer Funktionenschar mit unendlich vielen Lösungen mit ähnlichem Verhalten.
* Anfangswertproblem der Differenzialgleichung:
    Die Lösung eines Anfangswertproblems ist die Lösung einer Differenzialgleichung unter Berücksichtigung eines vorgegebenen Anfangswertes.
    Der Anfangswert $y_{0}$ für $x=0$ wird stets vorgegeben.
* Analytische Funktion (falls zu Vergleichszwecken benötigt):
    Für eine geschlossene Lösung einer gegebenen Differenzialgleichung wird die Stammfunktion (Integration)
       $\int f(x)dx=F(x)+C$
    gebildet.
    Für die Bestimmung der analytischen Funktion wird die Integrationskonstante $C$ berechnet, indem in die Gleichung der Stammfunktion anstelle der Größe $y(x)$ der Anfangswert $y_{0}$ gesetzt wird.
    Bei komplizierteren Differenzialgleichungen kann nicht immer die analytische Funktion durch Integration bestimmt werden. Die Konstante $C$ tritt häufig in der gesuchten Lösung nicht immer additiv, sondern auch faktoriell auf.
* Differenzengleichung mit dem Vorwärts-Differenzenquotient:
    Mit Hilfe von Differenzengleichungen kann die Differentialgleichung gelöst werden.
    Wird die Ableitung einer gewöhnlichen Differenzialgleichung $y'(x)=f(x,y)$ durch den Vorwärts-Differenzenquotienten ersetzt,
       $y'(x)=\lim _{x\to 0}{\frac {\Delta y}{\Delta x}}\approx {\frac {y(x+h)-y(x)}{(x+h)-x}}:={\frac {y_{k+1}-y_{k}}{h}}$.
    entsteht die explizite Differenzengleichung
       ${\frac {y_{k+1}-y_{k}}{h}}=f(x_{k},y_{k})$.
    Allgemeine Form der Differenzengleichung 1. O. nach dem Vorwärts-Differenzenquotienten (entspricht: „Euler-Vorwärts“):
* Entwicklung der Differenzengleichung für die oben gegebene Differenzialgleichung:
       ${\frac {y_{k+1}-y_{k}}{h}}=y_{k}+e^{x_{k}}$
    Differenzengleichung:
    $y_{k+1}=y_{k}+h\cdot (y_{k}+e^{x_{k}})$
    Der gewünschte Anfangswert wird in der nummerierten Tabelle für $y_{0}$ bei $k=0$ eingegeben.
$\to$ siehe auch Anwendung Differenzengleichung (Differenzenverfahren)
$\to$ siehe auch Artikel Explizites Euler-Verfahren

## Gewöhnliche Differenzenquotienten höherer Ableitungs- und Fehlerordnung

Neben der Approximation der Ableitung erster Ordnung, existieren auch Differenzenquotienten zur numerischen Berechnung höherer Ableitungen. Dazu werden in diesem Abschnitt ausschließlich zentrale Differenzenquotienten betrachtet. Analoge Überlegungen existieren auch für den Vorwärts- und der Rückwärtsdifferenzenquotienten. Die Grundlage zur Herleitung solcher Differenzenquotienten ist die Taylor-Reihe. Weiterhin existieren auch Differenzenquotienten mit einer höheren Fehlerordnung.
Für die zweite Ableitung kann zum Beispiel der Zusammenhang
    ${\frac {\Delta ^{2}y}{\Delta x^{2}}}:={\frac {y_{i+1}-2y_{i}+y_{i-1}}{\Delta x^{2}}}={\frac {f(x+\Delta x)-2f(x)+f(x-\Delta x)}{\Delta x^{2}}}=f''(x)+{\mathcal {O}}(\Delta x^{2})$
verwendet werden. Die hinter der ${\mathcal {O}}$-Notation stehende Wert kann dabei von $x$ abhängig sein. In der nachfolgenden Tabelle sind einige gewöhnliche, zentrale Differenzenquotienten höherer Ableitungsordnung angegeben. Die Tatsache, dass bei ungerader Ableitungsordnung der Funktionswert $y_{i}$ nicht vorhanden ist, geht auf das Prinzip der zentralen Differenzenquotienten zurück, bei welchem durch Mittelwertbildung die Fehlerordnung erhöht ist. Die Differenzenquotienten mit gerader Ableitungsordnung sind hier mit der minimalen Fehlerordnung angegeben. Diese lässt sich durch hinzunahme weiterer Funktionswerte erhöhen.
Gewöhnliche Zentrale Differenzenquotienten höherer Ableitungsordnung :
| Ableitungsordnung | Formel des Differenzenquotienten |
| $1$ | ${\frac {y_{i+1}-y_{i-1}}{2\Delta x}}$ |
| $2$ | ${\frac {y_{i+1}-2y_{i}+y_{i-1}}{\Delta x^{2}}}$ |
| $3$ | ${\frac {y_{i+2}-2y_{i+1}+2y_{i-1}-y_{i-2}}{2\Delta x^{3}}}$ |
| $4$ | ${\frac {y_{i+2}-4y_{i+1}+6y_{i}-4y_{i-1}+y_{i-2}}{\Delta x^{4}}}$ |
| $5$ | ${\frac {y_{i+3}-4y_{i+2}+5y_{i+1}-5y_{i-1}+4y_{i-2}-y_{i-3}}{2\Delta x^{5}}}$ |
| $6$ | ${\frac {y_{i+3}-6y_{i+2}+15y_{i+1}-20y_{i}+15y_{i-1}-6y_{i-2}+y_{i-3}}{\Delta x^{6}}}$ |

### Rekursionsgleichung

Die Berechnung der höheren gewöhnlichen, zentralen Differenzenquotienten kann mit Hilfe der nachfolgenden Rekursionsgleichung durchgeführt werden. Dabei repräsentiert $i$ den Index der Ortskoordinate $x_{i}$ und $n$ den Index der aktuellen Ableitungsordnung. Gestartet wird mit $n=1$ und folglich mit der Rekursionsgleichung für ungerade $n$.
    $\left.{\frac {\mathrm {d} ^{n}y}{\mathrm {d} x^{n}}}\right|_{x=x_{i}}\approx y_{i}^{(n)}={\begin{cases}\displaystyle {\frac {y_{i+1}^{(n-1)}-y_{i-1}^{(n-1)}}{2\Delta x}}&n:{\text{ungerade}}\\\\\displaystyle {\frac {y_{i+1}^{(n-2)}-2y_{i}^{(n-2)}+y_{i-1}^{(n-2)}}{\Delta x^{2}}}&n:{\text{gerade}}\end{cases}}\quad {\text{mit}}\quad n\in \mathbb {N} ,\;n>0$

### Summendarstellung

Die gewöhnlichen, zentralen Differenzenquotienten können weiterhin mit einer endlichen Summe dargestellt werden. Die Struktur dieser Formel besitzt eine direkte Verbindung zum pascalschen Dreieck beziehungsweise den Binomialkoeffizienten. Die Summendarstellung lässt sich mittels der obigen Rekursionsgleichung herleiten. Der Index $i$ repräsentiert die Ortskoordinate, zu welcher der Differenzenquotient ausgewertet wird. Die Summendarstellung von Ableitungen ungerader Ordnung $n$ beinhaltet die Methode des zentralen Differenzenquozentien, daher der Vorfaktor $1/2$.
    $\left.{\frac {\mathrm {d} ^{n}y}{\mathrm {d} x^{n}}}\right|_{x=x_{i}}\approx y_{i}^{(n)}={\frac {1}{\Delta x^{n}}}\cdot {\begin{cases}\displaystyle \sum _{k=0}^{n}\left[(-1)^{k}{\begin{pmatrix}n\\k\end{pmatrix}}\cdot y_{i+k-n/2}\right]&\quad n{\text{ ist gerade}}\\\\\displaystyle {\frac {1}{2}}\sum _{k=0}^{n-1}\left[(-1)^{k}{\begin{pmatrix}n-1\\k\end{pmatrix}}\cdot \left(y_{i+k+1-(n-1)/2}-y_{i+k-1-(n-1)/2}\right)\right]&\quad n{\text{ ist ungerade}}\end{cases}}$
mit $n\in \mathbb {N}$ und ${\begin{pmatrix}n\\k\end{pmatrix}}:={\frac {n!}{k!\cdot (n-k)!}}$.

### Produktdarstellung

Ausgehend von der obigen Rekursionsgleichung zur Berechnung von gewöhnlichen, zentralen Differenzenquotienten lässt sich eine Matrix-Produkt-Darstellung herleiten. Im ersten Schritt ist dazu eine Produktgleichung für die geradzahligen Ableitungen zu bestimmen, da in diesem Fall die zugehörige Rekursionsgleichung im Gegensatz zu den ungeraden Ableitungen eine geschlossene Kette bildet. Die Elemente der Matrizen $A^{l}$ sind wie folgt definiert und von der Dimension $[(2l+1)\times (2l+3)]$. Die Matrizen $A^{l}$ entsprechen genau der Signatur der obigen Rekursionsgleichung für gerade $n$.
    $A_{ij}^{l}:={\begin{cases}1&(i,j)\in \left\{(\nu ,\nu )\lor (\nu ,\nu +2)\;|\;1\leq \nu \leq 2l+1\right\}\\-2&(i,j)\in \left\{(\nu ,\nu +1)\;|\;1\leq \nu \leq 2l+1\right\}\\0&{\text{sonst}}\end{cases}}$
Der nachfolgende Vektor ${\underline {y}}_{k}$ beinhaltet die Funktionswerte $y_{m}=y(x_{i}+m\cdot \Delta x)$.
    ${\underline {y}}_{k}:={\begin{bmatrix}y_{i+k}&\cdots &y_{i}&\cdots &y_{i-k}\end{bmatrix}}^{T}$
Damit lässt sich die Näherung der $2k$-ten Ableitung im Punkt $x=x_{i}$ wie folgt darstellen.
    $\left.{\frac {{\text{d}}^{2k}y}{{\text{d}}x^{2k}}}\right|_{x=x_{i}}\approx {\frac {1}{\Delta x^{2k}}}\prod _{l=0}^{k-1}A^{l}\cdot {\underline {y}}_{k}\quad ,k\in \mathbb {N}$
Mit Hilfe der Matrizen $B^{k}$, mit der Dimension $[(2k-1)\times (2k+1)]$, findet sich ebenfalls eine Produktdarstellung für ungerade Ableitungsordnungen. Die Matrizen $B^{k}$ entsprechen genau der Signatur der obigen Rekursionsgleichung für ungerade $n$.
    $B_{ij}^{k}:={\begin{cases}1&(i,j)\in \left\{(\nu ,\nu )\;|\;1\leq \nu \leq 2k-1\right\}\\-1&(i,j)\in \left\{(\nu ,\nu +2)\;|\;1\leq \nu \leq 2k-1\right\}\\0&{\text{sonst}}\end{cases}}$
    $\left.{\frac {{\text{d}}^{2k-1}y}{{\text{d}}x^{2k-1}}}\right|_{x=x_{i}}\approx {\frac {1}{2\Delta x^{2k-1}}}\prod _{l=0}^{k-2}A^{l}\cdot B^{k}\cdot {\underline {y}}_{k}$

### Gewöhnliche Differenzenquotienten höherer Ableitungs- und Fehlerordnung
#### Zentrale Differenzenquotienten

Durch geschickte Anwendung der Taylor-Reihe (bzw. Taylor-Polynome) findet sich eine Matrizen-Gleichung zur Berechnung von Differenzenquotienten. Als Ansatz dient dazu die folgende Taylor-Approximation einer $2N$-fach differenzierbaren Funktion $y(x)$. Die Verwendung der oberen Grenze $2N$ der Summe bietet sich aufgrund der größeren Symmetrie an.
    $y(x)\approx \sum _{n=0}^{2N}{\frac {(x-x')^{n}}{n!}}\cdot {\frac {{\text{d}}^{n}y}{{\text{d}}x^{n}}}(x')$
Ausgehend von dieser Näherung von $y(x)$ sind die Substitutionen $x'=x_{i},\quad x=x_{i}+\nu \cdot \Delta x$ durchzuführen. Dies hat, wie in der nachfolgenden Gleichung zu sehen, zur Folge, dass die gesuchten Ableitungen der Funktion $y(x)$ am Ort $x_{i}$ vorhanden sind. Weiterhin ist hier zur Verkürzung die Index-Notation $y_{i,\nu }:=y(x_{i}+\nu \cdot \Delta x)\;,\quad y_{i}^{(n)}:={\tfrac {{\text{d}}^{n}y}{{\text{d}}x^{n}}}(x_{i})$ verwendet.
    $y_{i,\nu }\approx \sum _{n=0}^{2N}{\frac {(\nu \cdot \Delta x)^{n}}{n!}}\cdot y_{i}^{(n)}$
Durch Verschiebung des Index $\nu$ findet sich schlussendlich das nachfolgende lineare Gleichungssystem zur Berechnung der Differenzenquotienten bis zur Ableitungsordnung $2N$. Interessant ist dabei die enge Verwandtschaft der System-Matrix zur Vandermonde-Matrix, welche z. B. von der Polynominterpolation bekannt ist.
    ${\begin{bmatrix}y_{i,-N}\\\vdots \\y_{i,N}\end{bmatrix}}={\begin{bmatrix}{\frac {(-N\cdot \Delta x)^{0}}{0!}}&\cdots &{\frac {(-N\cdot \Delta x)^{2N}}{(2N)!}}\\\vdots &&\vdots \\{\frac {(N\cdot \Delta x)^{0}}{0!}}&\cdots &{\frac {(N\cdot \Delta x)^{2N}}{(2N)!}}\\\end{bmatrix}}\cdot {\begin{bmatrix}y_{i}^{0}\\\vdots \\y_{i}^{(2N)}\end{bmatrix}}$
In der nachfolgenden Tabelle sind einige Lösungen dieses Gleichungssystems angegeben. Zu beachten ist, dass für große $N$ die Matrix singulär wird und folglich die Matrix-Inversion am Rechner nicht mehr durchführbar ist. Neben den hier angegebenen Differenzenquotienten, welche in die Klasse der zentralen DZQ's einzuordnen sind, existieren auch andere Varianten.
Gewöhnliche zentrale Differenzenquotienten höherer Ableitungs- und Fehlerordnung :
| | $N=1$ | $N=2$ | $N=3$ |
| $n=1$ | ${\frac {y_{i+1}-y_{i-1}}{2\Delta x}}$ | ${\frac {-y_{i+2}+8y_{i+1}-8y_{i-1}+y_{i-2}}{12\Delta x}}$ | ${\frac {y_{i+3}-9y_{i+2}+45y_{i+1}-45y_{i-1}+9y_{i-2}-y_{i-3}}{60\Delta x}}$ |
| $n=2$ | ${\frac {y_{i+1}-2y_{i}+y_{i-1}}{\Delta x^{2}}}$ | ${\frac {-y_{i+2}+16y_{i+1}-30y_{i}+16y_{i-1}-y_{i-2}}{12\Delta x^{2}}}$ | ${\frac {2y_{i+3}-27y_{i+2}+270y_{i+1}-490y_{i}+270y_{i-1}-27y_{i-2}+2y_{i-3}}{180\Delta x^{2}}}$ |
| $n=3$ | | ${\frac {y_{i+2}-2y_{i+1}+2y_{i-1}-y_{i-2}}{2\Delta x^{3}}}$ | ${\frac {-y_{i+3}+8y_{i+2}-13y_{i+1}+13y_{i-1}-8y_{i-2}+y_{i-3}}{8\Delta x^{3}}}$ |
| $n=4$ | | ${\frac {y_{i+2}-4y_{i+1}+6y_{i}-4y_{i-1}+y_{i-2}}{\Delta x^{4}}}$ | ${\frac {-y_{i+3}+12y_{i+2}-39y_{i+1}+56y_{i}-39y_{i-1}+12y_{i-2}-y_{i-3}}{6\Delta x^{4}}}$ |
| $n=5$ | | | ${\frac {y_{i+3}-4y_{i+2}+5y_{i+1}-5y_{i-1}+4y_{i-2}-y_{i-3}}{2\Delta x^{5}}}$ |
| $n=6$ | | | ${\frac {y_{i+3}-6y_{i+2}+15y_{i+1}-20y_{i}+15y_{i-1}-6y_{i-2}+y_{i-3}}{\Delta x^{6}}}$ |

#### Differenzenquotienten zu beliebigen Stützstellen

Weiterhin besteht die Möglichkeit, Differenzenquotienten mit beliebige Stützstellen zu berechnen. Generell lässt sich ein Differenzenquotient mit der nachfolgenden Summe darstellen. Die Konstanten $s_{n}\in \mathbb {R} ,\;s_{i}\neq s_{j}$ entsprechen dabei den Stützstellen mit Verschiebung um $x_{i}\in \mathbb {R}$. Der Index $m$ entspricht der Ableitungsordnung. Die kleinste Genauigkeit ergibt sich bei $N\equiv m+1$. Durch Hinzunahme weiterer Stützstellen kann die Genauigkeit erhöht werden. Die weiter oben angegebenen zentralen Differenzenquotienten sind ein Spezialfall der hiesigen Betrachtung.
    $\left.{\frac {{\text{d}}^{m}y}{{\text{d}}x^{m}}}\right|_{x=x_{i}}\approx \sum _{n=0}^{N-1}C_{n}^{(m,N)}\cdot y(x_{i}+s_{n})\;,\quad N>m$
Die Koeffizienten $C_{n}^{(m,N)}$ berechnen sich durch Lösung des folgenden linearen Gleichungssystems, wobei $\delta _{i,j}$ das Kronecker-Delta repräsentiert.
    ${\begin{bmatrix}s_{1}^{0}&\cdots &s_{N}^{0}\\\vdots &&\vdots \\s_{1}^{N-1}&\cdots &s_{N}^{N-1}\end{bmatrix}}\cdot {\begin{bmatrix}C_{0}^{(m,N)}\\\vdots \\C_{N-1}^{(m,N)}\end{bmatrix}}=m!\cdot {\begin{bmatrix}\delta _{0,m}\\\vdots \\\delta _{N-1,m}\end{bmatrix}}$
Werden äquidistante Stützstellen $s_{n}=k_{n}\cdot \Delta x\;,\;k_{n}\in \mathbb {Z}$ gewählt, so stellt sich das lineare Gleichungssystem wie folgt dar.
    ${\begin{bmatrix}k_{1}^{0}&\cdots &k_{N}^{0}\\\vdots &&\vdots \\k_{1}^{N-1}&\cdots &k_{N}^{N-1}\end{bmatrix}}\cdot {\begin{bmatrix}C_{0}^{(m,N)}\\\vdots \\C_{N-1}^{(m,N)}\end{bmatrix}}={\frac {m!}{\Delta x^{m}}}\cdot {\begin{bmatrix}\delta _{0,m}\\\vdots \\\delta _{N-1,m}\end{bmatrix}}$