如何在不使用嵌套函数的情况下求解ODE?

我有一些微分方程需要用MATLAB的ODE解算器来解。虽然微分方程本身相当简单,但它们依赖于许多“常数”。这些常量不是通用的,需要调用方提供。

此类ODE的一个示例是:

dx/dt = -j * (k + x) ./ (l + x)

其中j、k和l是常数,x是变量。

function [outputVector] = someFunction(x, y, j, k, l, m, n, o)
    function [output] = someFunctionStep(t, inputVector)
        x = inputVector(1);
        y = inputVector(2);
        dx = -j .* (k + x) ./ (l + x);
        dy = -m .* (n + y) ./ (o + y);
        output = [dx;dy]
    end
    outputVector = ode15s(@someFunctionStep, [0, endTime], [x,y]);
end

然而,随着变量数量和代码大小的增加,这变得越来越不优雅,并导致几乎无法读取的代码混乱。所以,我想做的是将每个系统的阶跃函数移动到它自己的文件中,而不必a)将常数传递给输入向量中的阶跃函数,或者b)使用全局变量。有没有什么合理的方法可以做到这一点,或者我应该接受它并编写丑陋的代码?