Подбор кинетической модели для экспериментальных данных

[alex7851]

Имеется задача подбора оптимальных кинетических констант для имеющихся экспериментальных данных.
Например, для цепочки превращений трех веществ $A \leftrightarrow B \to C$ кинетические законы (изменение концентраций $A, B, C$ от времени $t$) имеют вид

$$\frac{dA}{dt} = -k_{1}A + k_{-1}B$$

$$\frac{dB}{dt} = k_{1}A - k_{-1}B - k_{2}C$$

$$\frac{dC}{dt} = k_{2}C$$

$$A(0)=A_{0}, \quad B(0)=0, \quad C(0)=0$$

и имеется набор экспериментальных точек вида

$$\left\lbrace t, A(t), B(t), C(t)\right\rbrace$$

где $t$ меняется с неодинаковым шагом. Нужно подобрать такие константы $k_{1}, k_{-1}, k_{2}$, чтобы согласие экспериментальных данных с получающейся кинетической моделью было наилучшим (например, в смысле наименьшей суммы квадратов между экспериментальными концентрациями и получающимися из модели). Для приведенной модели точные решения известны, но в общем случае это не так. В идеале требуется для одного набора данных попробовать разные модели (например, $A \leftrightarrow B \to C$ и $B \leftrightarrow A \to C$ или не первый порядок по веществу в каком-либо уравнении) с тем, чтобы выбрать наиболее подходящую.
Какими программами можно решать такую задачу?
Заранее спасибо.