 The likelihood equations
in these models are typically non-linear and estimates of the 
parameters are
obtained by numerical methods.Early in Fournier's work, it became clear that general numerical
solutions to these likelihood problems could only be reliably
achieved using function minimization algorithms that
incorporate accurate information about the gradients of the likelihood
surface.Computing the gradients (i.e.