вычислить символьный градиент в R

Question

Мне нужно вычислить градиент следующей функции в R .Затем я хочу использовать его в функции ODE45.

f(a,b)=log(|Q|)+log(Z'*Inv(Q)*Z)

, где:

1. Z - известный вектор (n by 1) с действительными элементами
2. Q является n by n симметричной матрицей.его элементы q[i,j] являются функциями a, b, который получается из другой функции, такой как q[i,j]=a*exp(-h[i,j]/b)), где известно h[i,j].
3. Inv(Q) является инверсией Q
4. |Q| - это определитель Q
5. n может быть больше 100.

Я могу сделать это в Matlab (с символической панелью инструментов), но я незнать, как это сделать в R .

Спасибо за все предложения!