function xdot = first( t, x)
% First: simple first-order dynamics

a = -1;
xdot = a*x;