function Anlaufvorgang2

global mM mF mue c k n Omega mT cT
mM = 50 ; mF = 100 ; mue = 0.1*0.01 ; c = 2.6e5 ; k = 3.75e3 ; n = 400/60 ; mT = 75 ; 
Omega = 2*pi*n ;
cT    = mT * Omega^2 ;

x0    = [0 ; 0 ; 0 ; 0] ;          % Anfangswerte
tspan = [0 1.5] ;                  % Zeitintervall

% Numerische Integration des Anfangswertproblems:
options = odeset('MaxStep' , 0.01) ; 
[t xx]  = ode45 (@RechteSeite , tspan , x0 , options) ;

% Grafische Ausgabe, zunaechst "Sortieren der Ergebnisse":
x  = xx (:,1) ;  
xT = xx (:,3) ;  

subplot(2,1,1) ; plot (t , x ) , grid on , title  ('x(t): Bewegung des Fundaments')
subplot(2,1,2) ; plot (t , xT) , grid on , title  ('x_T(t): Bewegung des Tilgers')

% ============================================================
% Funktion, die die "rechte Seite" des Dgl.-Systems definiert: 
function f = RechteSeite (t , xx)

global mM mF mue c k n Omega mT cT

x  = xx(1) ;  
v  = xx(2) ;  
xT = xx(3) ;  
vT = xx(4) ;  

xp  = v ;  
vp  = (mue*Omega^2*cos(Omega*t) - k*v - c*x + cT*(xT-x))/(mF+mM) ;
xTp = vT ;
vTp = -cT*(xT-x)/mT ;

f = [xp ; vp ; xTp ; vTp] ;