Next phase to implement BHHH algorithm