% Two first-order ODEs to solve the mass-spring-damper system
function yprime = fy1(t,y)

global m b k

% Define the forcing function, f(t), here
F = 10000;                % Units are Newtons

% Define the equations of state 
% y(1) = position (m)
% y(2) = speed (m/s)

yprime(1) = y(2);
yprime(2) = F/m - k/m * y(1) - b/m * y(2);

yprime = yprime';




