intermediate version: Homeier algorithm

http://dx.doi.org/10.1016/j.cam.2003.12.041
interesting modification of Newton,
but still does not solve issue #4