die Lösungsstrategie ist richitg gewählt, jedoch habe ich bedänken was die Lösung des homogenen Teils angeht.
Behauptung: yhom (x) = c1exp( ω0 * (sqrt(μ2 - ω02 ) -μ) ) + c2 exp( ω0 * (-sqrt(μ2 - ω02 ) -μ) )
Hier sollte jedoch der Definitionsbereich: 0 ≤ μ < ω0 bedacht werden.
Wähle als Ansatz y = exp(λx), y' = (y)', y" = (y')' (prüfe durch Rechnung)
Setzt man y, y', y" in den Homogenen Teil der DGL, so erkennt man:
exp(λt) * ( λ2 + 2μλ + ω02 ) = 0 dann wenn λ2 + 2μλ + ω02 = 0 (löse mit PQ-Formel)
λ1/2 = -μ +- sqrt(μ² - ω02 ) | ACHTUNG: μ² - ω02 < 0, da 0 ≤ μ < ω0
Daraus folgt, dass wir uns im ℂ befinden und hier geschickt Lösen müssen um wieder in ℝ zu landen.
Danach sollte der weitere Lösungsweg benutzt werden können.