% Beispiel0103
%
% Dieses m-file ist Teil der Arbeitsmaterialien des 
%   Schuelerseminars "Bewegungssimulation mit dem Computer"
% am Fachbereich Mathematik und Informatik der
% Martin-Luther-Universitaet Halle-Wittenberg, Deutschland
%
% Author :      Prof. Dr. M. Arnold, martin.arnold@mathematik.uni-halle.de
% Version of :  Sep 30, 2005
%
% Freier Fall
%
% Beispiel:
%   beispiel0103
%   print -dpng ../png/beispiel0103.png

% -> physikalische Parameter
t0 = 0.0;               % Anfangszeit
y0 = 3.0;               % Anfangshoehe
g  = 9.81;              % Fallbeschleunigung

% -> Ober- und Untergrenze Ordinatenachse
ymin = -10.0;   ymax = 10.0;

% -> Initialisierung der graphischen Darstellung
h = plot ( t0, y0, '.g', 'MarkerSize', 24 );
xlabel ( 'x' );
ylabel ( 'y' );
title ( 'Freier Fall' );
axis ( [ -1 1 ymin ymax ] );

% -> Berechnung der Ergebnisse und Animation
t = t0;   y = y0;       % Initialisierung
while y>ymin,
    t = t + 0.005;      % Zeitschritt von 5 ms
    y = y0 - g*t^2/2;   % neue Hoehe
    set ( h, 'YData', y );
    title ( sprintf('Freier Fall  t = %4.3f [s]',t) );
    drawnow;
end;

