%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                                                                      %
%   Gewoehnliche Differentialgleichungen                               %
%                                                                      %
%    -- Eine Einfuehrung aus der Perspektive der Dynamischen Systeme   %
%                                                                      %
%                                             Vieweg + Teubner, 2009   %
%                                                                      %
%                                                                      %
%                                                                      %
%       MATLAB-Beispiele aus Anhang B.2 - Numerische Loesungen         %
%                                                                      %
%         (benoetigt zusaetzlich 'bsp1_f.m' und 'bsp2_f.m')            %
%                                                                      %
%                                   (C) Lars Gruene, Oliver Junge 2008 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


% Loesung von Beispiel 1

ode45(@bsp1_f, [0,1],0.5)

figure;

% Loesung von Beispiel 1 mit Ausgabe auf groeberem Gitter

ode45(@bsp1_f, [0:0.1:1],0.5)

figure;

% Loesung von Beispiel 2

ode45(@bsp2_f, [0,7],[1, 0])

figure;

% Beispiel fuer die Verwendung von "options"

options = odeset('RelTol',1e-13);
options = odeset(options,'AbsTol',1e-13);

ode45(@bsp1_f, [0,1],0.5, options);

figure;

% Vergleich ode113 - ode45 (mit Ausgabegitter = Rechengitter)

options = odeset(options,'refine',1);

ode113(@bsp1_f, [0,1],0.5, options);

figure;

ode45(@bsp1_f, [0,1],0.5, options);


% Ausgabe von Gitter und Gitterfunktionen in Variablen

[t1,x1] = ode45(@bsp1_f, [0:0.1:1],0.5);
[t2,x2] = ode45(@bsp2_f, [0,7],[1, 0]);


% Die Beispiele zur grafischen Darstellung finden sich in
% "plotbeispiele.m"
