function xp = linearf(t,x) %evaluates a particular linear ODE/vector field global A; xp = A*x;