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