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

global m b k

% define the forcing function, f(t), here

f = 10;

% define the rate eauations 
%
% 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';
