global m rR A cw roh iA iG eta g FR x_profile h_profile
global tend regelung_aktiv vsoll

% +: Parameter initialisieren ************************
m   = 1600;     % kg Fahrzeugmasse
rR  = 0.291;    % m dynamischer Reifenradius
A   = 1.7;      % m² Stirnfläche
cw  = 0.32;     % Luftwiderstandsbeiwert
roh = 1.2;      % kg/m³ Luftdichte
iA  = 4.06;     % Achsantrieb
iG  = 1.3;      % Getriebeübersetzung 5. Gang
eta = 0.9;      % Wirkungsgrad Antriebsstrang
g   = 9.81;     % m/s² Erdbeschleunigung
FR  = 160;      % N Rollreibungskraft

% x_profile = input('Geben Sie den Vektor für den Weg x in Metern ein: ');
% h_profile = input('Geben Sie den Vektor für die Hoehe h in Metern ein: ');

x_profile = [0  1  10  20 30 10000]*1000; % Strecke [m]
h_profile = [0  0 300 300  0     0];      % Höhe [m]
% -: Parameter initialisieren ************************

% +: Anfangsbedingungen setzen ************************
x0 = 0; % Anfangsbedingung: Position
v0 = 0; % Anfangsbedingung: Geschwindigkeit
% -: Anfangsbedingungen setzen ************************

% +: Simulationsparameter setzen ************************
tend = 60*10; % 10 min.
regelung_aktiv = 0;
vsoll = 80/3.6;
% -: Simulationsparameter setzen ************************

[t,y] = ode23('vehicle_ohneReg',[0:tend]',[x0; v0]);
figure(1), plot(t,y(:,1)/1000), title('Strecke in km'),             xlabel('t [sec]');
figure(2), plot(t,y(:,2)*3.6),  title('Geschwindigkeit in km/h'),   xlabel('t [sec]');