It uses U0 but never explains what it is. It also calls a the "shift" but never explains what it is. Perhaps some notes explaining this stuff would make the solver more transparent. And a reference to further material that uses this approach would be great if it exists. Barry