function Anlaufvorgang1

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

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

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

x = xx (:,1) ; % ... numerisch ermitteltes Ergebnis x(t)  

% Analytische Lösung:
omega  = sqrt(c/(mM + mF)) ;
eta    = Omega/omega ;
D      = k/(2*sqrt((mM+mF)*c)) ;

phi    = atan((2*D*eta)/(1-eta^2)) ;
if (phi < 0) phi = phi+pi ; end

omegaD = omega*sqrt(1-D^2) ;
A      = mue*Omega^2/(c*sqrt((1-eta^2)^2+4*D^2*eta^2)) ;
A1     = - A*cos(phi) ;
A2     = - A*(D*omega*cos(phi)+Omega*sin(phi))/omegaD ;

xhom   = exp(-D*omega*t).*(A1*cos(omegaD*t)+A2*sin(omegaD*t)) ;
xpart  = A*cos(Omega*t-phi) ;
xan    = xhom + xpart ;        % ... analytisch ermitteltes Ergebnis x(t)

plot (t , xan , t , x) , grid on , title  ('x(t), analytische und numerische Lösung')

% ============================================================
% Funktion, die die "rechte Seite" des Dgl.-Systems definiert: 
function f = RechteSeite (t , xx)
global mM mF mue c k Omega
x   = xx(1) ;  
v   = xx(2) ;  
xp   = v ;  
vp   = (mue*Omega^2*cos(Omega*t) - k*v - c*x)/(mF+mM) ;  
f = [xp ; vp] ;