horizon radius by newton-raphson method