function F0=expl_euler(h)
%% Parameter festlegen

%Anfangswert b

b = 3;

%Anfangswert x0
x0 = 2;

% Simulationszeit
t_start = 0;
t_stop = 5;

% Schrittweite h
%h = 0.05;
%% Berechnung

% Maximaler Index
i_max = round(t_stop/h);

% Speicherplatz reservieren (Vektoren anlegen)
t_sim = zeros(1, i_max+1);
x_sim = zeros(1, i_max+1);
x_exakt = zeros(1, i_max+1);

% Initialisierung
t_sim(1) = t_start;
x_sim(1) = x0;

% Schleife
%
% Hier entspricht
% i ~ t (Index i entspricht Zeitpunkt t "aktuelle Zeit")
% i+1 ~ t + h = t + 1*h
% t_sim(i) = t
% t_sim(i+1) = t + h
%
% x_sim(i) = x(t)
% x_sim(i+1) = x(t+h)
%


for i = 1:i_max
   
    % Neuer Simulationszeitpunkt
    % t+h = t + h
    % t_sim(i+1) = t        + h
    t_sim(i+1)   = t_sim(i) + h;
   
    % Funktionswert berechnen: f(x(t),t)
    % *************************************************************
    % *************************************************************
    % FUNKTIONSSPEZIFISCH !!!!
    %
    % f(x(t),t) =
    fxt = x_sim(i)+b*t_sim(i);
    %
    % *************************************************************
    % *************************************************************
   
    % Neuer x-Wert berechnen x(t+h)
    % x(t+h)   =   x(t)   + f(x(t),t)*h
    x_sim(i+1) = x_sim(i) +    fxt*h;
   
    
   
end

  x_exakt = (b+2)*exp(t_sim)-b*(t_sim+1);

% *************************************************************
% *************************************************************
%% Berechnung des Fehlers (Exakte Lösung - Simulationslösung)

sim_error = x_exakt - x_sim;
F_max=sim_error(end);

% *************************************************************
% *************************************************************
%% Berechnung des Fehlers (Exakte Lösung - Simulationslösung)

sim_error = x_exakt - x_sim;
F_max=sim_error(end);

F_soll=0.025;
F0=F_max-F_soll;



end