Thanks for your great article, but one flaw I should mention is that in your implementation, the `sqrt` has been placed on wrong position compared to your formula.