PROS_Exercises1_msd

function dx = msd(t, x, m, b, k, u)

% Παραγωγος θεσης
dx(1) = x(2);
% Παραγωγος ταχυτητας
dx(2) = 1/m * (-k*x(1) - b * (x(2))^2 + u(t));
% Για να επιστρεφει στηλη
dx = dx';

end