Create "gradient.R"

A private function.
It computes the gradient of ODE solution with respect to ODE parameters.