diff --git a/README.md b/README.md
index 241be25..d7125ca 100644
--- a/README.md
+++ b/README.md
@@ -121,7 +121,7 @@ result = miex.shexqnn2(x=1.0, m=complex(1.5, 0.0), nang=91, doSA=True)
 | ------------- | -------------------------------------------------------- | --------------- |
 | `x`           | Size parameter                                           | float           |
 | `m`           | Complex refractive index                                 | complex         |
-| `nang=1`      | Half number of scattering angles in the interval [0, 90] | int, optional   |
+| `nang=2`      | Half number of scattering angles in the interval [0, 90] | int, optional   |
 | `doSA=False`  | Calculate scattering amplitude functions                 | bool, optional  |
 | `nterm=2e7`   | Maximum number of terms to be considered                 | int, optional   |
 | `eps=1.0e-20` | Accuracy to be achieved                                  | float, optional |
diff --git a/src/miex/miex.py b/src/miex/miex.py
index 8492a9d..b4d5a4a 100644
--- a/src/miex/miex.py
+++ b/src/miex/miex.py
@@ -73,7 +73,7 @@ def aa2(ax, ri, num, ru):
 
 
 @njit(cache=True)
-def shexqnn2(x, ri, nang=1, doSA=False, nterm=2e7, eps=1.0e-20, xmin=1.0e-06):
+def shexqnn2(x, ri, nang=2, doSA=False, nterm=2e7, eps=1.0e-20, xmin=1.0e-06):
     """Derive quantities for a single size parameter and chemical composition.
     This routine is based on the routine 'shexqnn' published by
     N.V.Voshchinnikov: 'Optics of Cosmic Dust', Astrophysics and Space Physics Review 12, 1 (2002)
@@ -86,7 +86,7 @@ def shexqnn2(x, ri, nang=1, doSA=False, nterm=2e7, eps=1.0e-20, xmin=1.0e-06):
     ri : complex float
         complex refractive index
 
-    nang : int, optional, default = 1
+    nang : int, optional, default = 2
         half number of scattering angles theta in the intervall 0...pi/2 (equidistantly distributed)
 
     doSA : bool, optional, default = False