diff --git a/misc/finitesqlatzsym-scatter.py b/misc/finitesqlatzsym-scatter.py
index f10082b..4da7ae7 100755
--- a/misc/finitesqlatzsym-scatter.py
+++ b/misc/finitesqlatzsym-scatter.py
@@ -236,7 +236,7 @@ N = positions.shape[0]
 kx = np.linspace(pargs.kxmin, pargs.kxmax, num=pargs.kxdensity, endpoint=True) * 2*np.pi / dx
 ky = np.linspace(pargs.kymin, pargs.kymax, num=pargs.kydensity, endpoint=True) * 2*np.pi / dy
 kx, ky = np.meshgrid(kx, ky, indexing='ij', copy=False)
-kz = np.sqrt(k_0 - (kx ** 2 + ky ** 2))
+kz = np.sqrt(k_0**2 - (kx ** 2 + ky ** 2))
 
 klist_full = np.stack((kx,ky,kz), axis=-1).reshape((-1,3))
 TMatrices_om = TMatrices_interp(freq)
diff --git a/setup.py b/setup.py
index feac6ff..5ab8365 100644
--- a/setup.py
+++ b/setup.py
@@ -34,7 +34,7 @@ qpms_c = Extension('qpms_c',
         )
 
 setup(name='qpms',
-        version = "0.2.14",
+        version = "0.2.15",
         packages=['qpms'],
 #        setup_requires=['setuptools_cython'],
         install_requires=['cython>=0.21','quaternion','spherical_functions','py_gmm'],