(e^(productlog((ln[n!/√(2π)]/e))+1))-1/2