Then, however, I see that the der(f, …) allows one or more input arguments to be listed, and then I would have expected that the result of Gibbs_T would be an array of the same size as the number of ...