Note that if I don’t register the derivatives, the gradient will not evaluate the partial derivatives of the lower incomplete gamma function with respect to x, y, and z, and the ODESystem
constructor will fail due to the appearance of multiple independent variables.