3: compute σ using approximation (23)