function  Solve_Feder_2( w, d, s, tm )
   y0 = [ s, 0 ];   % Anfangswerte
   tspan = [0,tm];  % Zeitintervall
   % Lösung des Anfangswertproblems
   [t,y] = ode45( @Fun_Feder_2, tspan, y0, [], w*w, d );
   z = y( :, 1 ); % Funktionswert = 1. Komponente
   plot( t, z ); grid on; % alles plotten
end