matlab ode45函数 常微分方程中参数为矩阵不是固定值怎么调整代码

[chen10217]

新手入门Matlab 实践练习完成一个存在竞争的两个种群增长的模型：

m_1为耐缺氧的细胞种群A的数量
m_2为不耐缺氧细胞种群B的的数量
r_1为A种群的固定增长率，
r_2为B种群的固定增长率
N_x为距离球心x处种群A,B的最大容量（目前设为定值）
σ_1 (O)为B种群相对于A的相对竞争参数，
σ_2 (O)为A种群相对于B的相对竞争参数，它们是含氧量O(s)的函数

我使用ode45解常微分方程，关键问题是O是求解函数后得到的1X3501的矩阵= = 实在改了半天不知道怎么调整。。
代码：

其中的fun函数：

现在出来的图实在是#￥%*&。。。：


直接设为常数画出来的图：


求大神指点！！1