function dzdt = ODEcar(z) 
%% simulation ODE   dv/dt = (F - k * v )/ m;
% Constant
m = 1000;  % mass  % mass

% Parameters
k = 50; % Gain (m/s % pedal)

F = z(1);
x = z(2);
v = z(3);

% ODE equation 
dxdt = v;
dvdt =  (F - k * v )/ m;

dzdt = [dxdt, dvdt];

end