function d=cclin(t,x)
% CCLIN(T,X) provides the right side for the constant coefficient linear 
% system x'=A*x, for purposes of solving with ODE23, ODE45, or BDODE45.  

   global A 
   d=A*x(:);