Frage zur Diskretisierung einer partiellen Differentialgleichung

[soerenP]

Hallo, Leute!
Ich versuche die Diskretisierung der partiellen Differentialgleichung einer biegesteifen Saite nachzuvollziehen, wie sie in einer wissenschaftlichen Arbeit durchgeführt wurde.
Sie sieht kontinuierlich so aus:
∂²y/∂t²=
c²∂²y/∂x²+
-εc²L²∂⁴y/∂x⁴
-2b1∂y/∂t
+2b3∂³/∂t³
+f(x,t)
Bis auf die dritte Ableitung nach der Zeit versteh ich auch alles. Die wiederum wird aber in der Arbeit diskret scheinbar unter anderem mit samples von anderen Orten gebildet. Also so:
d∂³y/∂t³=1/2*(-2y[x][t+1]+2y[x][t-1]+y[x+2][t]+y[x-2][t]-y[x+1][t-1]-y[x-1][t-1]-y[x][t-2])/deltaT^3
Zumindest sind das die Werte, die zur Berechnung herangezogen werden.
Weiß jemand, was da passiert ist?
Vielen Dank und viele Grüße
Sören