function Aufg28_5
global   m c my a g
m = 20 ; c = 300 ; my = 0.8 ; a = 0.5 ; g = 9.81 ;

x0 = [2*a ; 0] ; tspan = [0 10] ;       % Anfangswerte, Zeitintervall
options = odeset ('MaxStep' , 0.1) ;    % Integration des AWPs: 
[tres xvres]  = ode23s (@RechteSeite , tspan , x0 , options) ;

x = xvres(:,1) ;  v = xvres(:,2) ;      % Sortieren der Ergebnisse 
subplot(2,1,1) ; plot (tres , x) , grid on , title  ('x-t-Diagramm:')
subplot(2,1,2) ; plot (tres , v) , grid on , title  ('v-t-Diagramm:')
End_x = x(end)
End_v = v(end)

% ============================================================
% Function, die die "rechte Seite" des Dgl.-Systems definiert: 
function f = RechteSeite (t , xv)
global m c my a g

x = xv(1) ;
v = xv(2) ;
if (x < a)  myq = 0 ; else myq = my ; end
if (v > 0)  r   = 1 ; else r   = -1 ; end

xPunkt = v ;
vPunkt = - c*x/m - r*myq*g ;

f = [xPunkt ; vPunkt] ;