我有一个Octave程序，可以将实验数据融入不同ODE的模型中。 Octave程序使用leasqr（）函数来拟合数据。
我的任务是在C ++中重新实现它，但我没有找到任何提供类似功能的库。 ALG和GSL都提供了在没有实验数据的情况下模拟ODE系统的功能，或者将数据拟合到单个ODE中。
ODE系统看起来像这样：
dA/dt = -k1*A
dB/dt = k1*A-k2*B-k3*B
dC/dt = k3*B-k4*C
有谁知道如何解决这个问题或使用什么库？我无法想象在Octave（或R）中这是可能的，但是在C ++之前没有人能够实现这一点。
感谢您的帮助！