U(\beta) = \sum_{i=1}^N \frac{\partial \mu_{ij}}{\partial \beta_k} V_i^{-1} \{ Y_i - \mu_i(\beta)\} \,\! .
The parameters \beta_k are estimated by solving U(\beta)=0 and are typically obtained via the Newton-Raphson algorithm.