function dxdt = linsystem(t, x)
% LINSYSTEM Implements the linear ODE:
%
%   x'=Ax
%
% where
%      [  998    1998 ]
%  A = [              ]
%      [ -999   -1999 ]

A = [998 1998; -999 -1999];
dxdt = A*x;