function [dxdt] = linode(t,x,flag,A)
% an odefile for a standard linear system dxdt = Ax

dxdt = A*x;