Implemented crude optimization of likelihood function from Ellison
(2007), constraints still missing. Gradients needed for performance.