function  Solve_Feder_3( w, d, F, w_ext, s, tm )
   % Anfangswertproblem lösen
   [t,y] = ode45( @Fun_Feder_3, [0,tm], [s,0], [], w*w, d, F, w_ext );
   z = y( :, 1 );                  % Funktionswert = 1. Komponente

   frequenz = w_ext / (2*pi);      % Frequenz der äußeren Schwingung
   x = F*sin(2*pi*frequenz*t)/50;  % äußere Schwingung als Vergleich
   plot( t, z, t, x, 'r:' );       % alles plotten
end