ode - 麦克斯韦斯特凡常微分方程

Question

我试图在膜上求解 Maxwell Stefan 方程，以获得膜厚度“z”上的瞬态摩尔分数分布。但不知何故，我无法使用 ODE45 对其进行编码，更有可能我无法编写使用 ODE45 解决的系统。如果有人可以帮助我了解主要的语法和功能，那就太好了。我试图解决的方程是

(dy_i)/dz=1/(cD_(i,j) ) [y_i (N_i+N_j )-N_i ]

其中 c 是浓度，D_(i,j) 是二元扩散系数。

先感谢您。