It is a very good approximation based on general relativity solution :) More specifically, it is approximation of Schwarzschild black hole potential by Paczynski and Wiita (1980):
F = -G * M / (R - Rg) 2
It is not a Kerr metric approximation, so spin is ignored. We still didn't find easy to ...