PROS - Ergasia 1 - Thema 1 - msd.m

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