function xout = mass(x)
%
% pendulum example
%
k=10;
n=10;
xout    = x;   % make sure shapes agree
xout(1) = x(2);
xout(2) = -k*k*x(1)-2*n*x(2);