matlab 求解微分方程时，系数矩阵随状态变化而变化怎么解？
 

例如：

主程序调用 ode45：
[t y] = ode45('linearODE',tspan,x0,options);

function dx = linearODE(t,x,flag)
% 如果 | x(i) - x(j) | < 1 则 A(i,j) 元素是1 否则为0
dx = A*x;

这里A(t)与方程演化的状态有关，请问有办法解吗？
另外 如果函数写成
function dx = linearODE(t,x,flag)
A = 某个矩阵;
dx = A*x;
x

matlab会输出一些向量x，与解出的 y 只有最后一维相同，剩下的不同，
请问这些输出的x是什么量呢？

谢谢~