Dieses JavaScript berechnet die partielle 1.Ableitung einer Funktion von drei Variablen
nach der ersten Variable x mit dem Romberg-Schema. Der Funktionswert wird mit ausgegeben.
Das Schema verwendet zentrale Differenzenquotienten
Dj(0)(f) := [f(x0+hj - f(x0-hj)]/[2hj]
und bildet mit ihnen Linearkombinationen
Dj(k)(f) := [4kDj+1(k-1)(f) - Dj(k-1)(f)]/[4k-1]
wobei j,k = 1,..,n. Der Ausdruck D0(n)(f) soll die Ableitung
approximieren. Das grösstmögliche n ergibt sich aus der Voraussetzung
f in C2n[a,b]. Ist der Fehler δ=|D0(n)-D1(n-1)|≤ε (Vorgabe),
wird das Ergebnis akzeptiert. Der gesuchte Näherungswert für f'(x0) ist
dann von der Güte
f'(x0) = D0(n) + O(h2n+2),
wobei hier n die Schrittzahl und h die Schrittweite bezeichnen.