Added Sb Boston-Sullivan calculation
Added indetermination handling for _Rlj(0)
Broyden second review