Although this maximization cannot be done in closed form, it can be optimized numerically.
effective are gradientbased methods that use approximate second-order information, such as conjugate gradient and limited-memory BFGS (Byrd et al., 1994).