麦克斯韦·斯特凡·奥德

Question

我正在尝试求解膜上的麦克斯韦·斯特凡方程，以获取膜厚度'z'上的瞬时摩尔分数分布。但是以某种方式我无法使用ODE45对其进行编码，更有可能是我无法编写使用ODE45解决的系统。如果有人可以在主要语法和功能上帮助我，那将是非常不错的。我想解决的方程是

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

其中c是浓度，D_（i，j）是二元扩散系数。

谢谢。