Make Newton-Raphson iterations more robust

[johnlees]

Should check that parameters aren't diverging, and that likelihood is constantly increasing (using step-halving if it isn't)