In der folgenden Beispielimplementierung wird zum Lösen der quartischen Gleichung die Funktion <tt>magnet::math::quarticSolve</tt> verwendet, deren Implementierung unter der unten angegebenen Referenz<ref>[https://github.com/toastedcrumpets/DynamO/blob/master/src/magnet/magnet/math/quartic.hpp ''dynamo - Event driven molecular dynamics simulator'']. Marcus Bannerman. 2011.</ref> zu finden ist. Es wird außerdem davon ausgegangen, dass <tt>Vector</tt> eine Vektorklasse ist, auf der die üblichen Operationen definiert sind. Die statische Methode <tt>dot</tt> muss das Skalarprodukt zweier Vektoren berechnen, die Methode <tt>lengthSq()</tt> das Quadrat der Vektorlänge (das Skalarprodukt mit sich selbst).