Sehr „harmonisch“ wird es so aber wohl nicht aussehen...
Da fehlt die Dämpfung und die Folgerichtigkeit der Bewegung.
Man könnte auch permanent Lissajous-Figuren berechnen (da ist das "harmonische" drin) und einfach den Radius mit einer natürlichen Dampfungsfunktion (exp(-ax)) multiplizieren. Ab Auslenkung < Epsilon abschneiden, dann isses wieder endlich.
Man könnte auch mit einer zeitlich umgekehrten Dämpfung einfaden, das wäre dann rundum harmonisch und glatt und hübsch.