DFfun function

Compute the first and second derivatives of the negative log likelihoods