求助---求解状态方程中遇到的问题

samson888

小弟最近在求解状态方程时遇到问题，还请大虾们帮忙指点；
问题是：状态方程 dx/dt=Ax+Bu，其中dx/dt表示x对时间t的一阶导数，A是6X6阶矩阵，B是6X3阶矩阵，U是输入变量（常数）但是矩阵A中包含有状态变量x，也就是说方程中的系数A是状态变量的函数，我考虑用4阶龙格库塔法（自编程序）解状态方程，采用定步长，因此，在计算中，每循环一次，方程系数必须重新更新计算一次，把新的状态变量赋到A中，所以必须把与状态变量有关的系数放在循环内，但是如何编写这个循环子程序，使得每次更新的系数释放出来重新回到循环中去呢？？？请高手指导，万分感谢