Nello schema centrato si effettua un passo a sinistra ed uno a destra. In pratica la formula generale della sommatoria assume la seguente forma: $$ {dy \over dx}\lvert_{x_j} = \sum_{i=-1}^{1} a_iy(x_{j+i}) + \epsilon(h^p) \hspace{2cm} (\heartsuit) $$ $$ \downarrow $$ $$ {dy \over dx}\lvert_{x_j} = a_{-1}y(x_{j-1}) + a_0y(x_{j}) + a_1y(x_{j+1}) + \epsilon $$

Ora, noi sappiamo che gli sviluppi di Taylor della \(y\) a destra ed a sinistra sono rispettivamente i seguenti:

• \( y(x_{j-1}) = y_{j-1} = \left[ y_j -y_j'h +{1\over 2}y''_j h^2 -{1\over 6}y'''_j h^3 + O(h^4) \right] \)
• \( y(x_{j+1}) = y_{j+1} = \left[ y_j +y_j'h +{1\over 2}y''_j h^2 +{1\over 6}y'''_j h^3 + O(h^4) \right] \)

Se ora andiamo a sostituire questi sviluppi nella sommatoria generale \( (\heartsuit) \) otteniamo la seguente espressione: $$ \small {dy \over dx}\lvert_{x_j} = a_{-1}\left[ y_j -y_j'h +{1\over 2}y''_j h^2 -{1\over 6}y'''_j h^3 + O(h^4) \right] + a_0y_{j} + a_1\left[ y_j +y_j'h +{1\over 2}y''_j h^2 +{1\over 6}y'''_j h^3 + O(h^4) \right] + \epsilon $$ $$ \small {dy \over dx}\lvert_{x_j} = a_{-1}y_j - a_{-1}y_j'h + a_{-1}{1\over 2}y''_j h^2 -a_{-1}{1\over 6}y'''_j h^3 + a_{-1}O(h^4) + a_0y_{j} + a_1y_j + a_1y_j'h + a_1{1\over 2}y''_j h^2 +a_1{1\over 6}y'''_j h^3+ a_1O(h^4) + \epsilon $$ Rielaborando (ossia, raggruppando \( y_j \), \(y'_j\) ecc), otteniamo: $$ \small {dy \over dx}\lvert_{x_j} = (a_{-1} + a_0 + a_1)y_j + (a_1 - a_{-1})y_j'h + (a_{-1} + a_1){1\over 2}y''_j h^2 + (a_1-a_{-1}){1\over 6}y'''_jh^3 + (a_{-1}+a_1)O(h^4) + \epsilon $$

In questa espressione per essere verificata i coefficienti \(a_{-1}\), \(a_0\) ed \( a_1 \) devono soddisfare al seguente sistema: $$ \begin{cases} a_{-1} + a_0 + a_1 = 0 \\ a_1 - a_{-1} = {1 \over h} \\ a_1 + a_{-1} = 0 \\ \end{cases} \rightarrow \begin{cases} a_{-1} + a_0 -a_{-1} = 0 \\ a_1 + a_1 = {1 \over h} \\ a_1 = -a_{-1} \\ \end{cases} \rightarrow \begin{cases} a_0 = 0 \\ a_1 = {1 \over 2h} \\ a_{-1} = -{1 \over 2h} \\ \end{cases} $$

Sostituendo i coefficienti nell'espressione, troviamo l'errore \( \epsilon\): $$ \small {dy \over dx}\lvert_{x_j} = (0)y_j + {1 \over h}y_j'h + (0){1\over 2}y''_j h^2 + {1 \over h}{1\over 6}y'''_jh^3 + (0)O(h^4) + \epsilon $$ $$ {dy \over dx}\lvert_{x_j} = y_j' + {1\over 6}y'''_jh^2 + \epsilon $$ E quindi l'errore ha una precisione \(p=2\) $$ \epsilon = -{1\over 6}y'''_jh^2 + O(h^3) $$

Ed in fine la formula dello schema centrato si ottiene, sostituendo il valore dei coefficienti nella formula generale \( (\heartsuit) \): $$ {dy \over dx}\lvert_{x_j} = a_{-1}y(x_{j-1}) + a_0y(x_{j}) + a_1y(x_{j+1}) + \epsilon $$ $$ \downarrow $$ $$ {dy \over dx}\lvert_{x_j} = -{1 \over 2h}y_{j-1} + 0y_j+ {1 \over 2h}y_{j+1} + \epsilon(h^2) = {y_{j+1}-y_{j-1} \over 2h} + \epsilon(h^2) $$