% 28: Kinetik des Massenpunktes
% Beispiel 2: Berechnung eines Gleitsteines S.501
% Author:  Eugen Lass
% Revision: Prof. Dr. Thomas Frischgesell
% Datum:   03.04.2003

% Lösung der Differentialgleichung

clear all;

options = odeset('MaxStep',0.1); % Ohne Schrittweitensteuerung erhalten wir ein falsches Ergebnis

% setze Anfangswerte, Zeitspanne
x0 = [-4.493; 0]; % Anfangswerte [xi ; xipunkt]
tspan = [0 20];  % Zeitspanne

% Integration der DGL
[t x] = ode45('S501_xpunkt', tspan, x0,options);

% grafische Ausgabe von xi und xipunkt
xi = [1 0]*x';
xipunkt = [0 1]*x';

plot(t,xi,t,xipunkt)

% legt die Achsen fest
grid on
xlabel('Zeit [s]');title('Vertikale Führung')
legend('xi [rad]','xipunkt [rad/s]',4)