diff --git a/GPy/kern/kern.py b/GPy/kern/kern.py
index bf8ba612..37a18f04 100644
--- a/GPy/kern/kern.py
+++ b/GPy/kern/kern.py
@@ -861,12 +861,14 @@ def kern_test(kern, X=None, X2=None, output_ind=None, verbose=False, X_positive=
         if X_positive:
             X = abs(X)
         if output_ind is not None:
+            assert(output_ind<kern.input_dim)
             X[:, output_ind] = np.random.randint(low=0,high=kern.parts[0].output_dim, size=X.shape[0])
     if X2==None:
         X2 = np.random.randn(20, kern.input_dim)
         if X_positive:
             X2 = abs(X2)
         if output_ind is not None:
+            assert(output_ind<kern.input_dim)
             X2[:, output_ind] = np.random.randint(low=0, high=kern.parts[0].output_dim, size=X2.shape[0])
 
     if verbose:
diff --git a/GPy/testing/kernel_tests.py b/GPy/testing/kernel_tests.py
index 5d2fbeec..0fceac60 100644
--- a/GPy/testing/kernel_tests.py
+++ b/GPy/testing/kernel_tests.py
@@ -36,7 +36,7 @@ class KernelTests(unittest.TestCase):
     def test_eq_sympykernel(self):
         if SYMPY_AVAILABLE:
             kern = GPy.kern.eq_sympy(5, 3)
-            self.assertTrue(GPy.kern.kern_test(kern, output_ind=3, verbose=verbose))
+            self.assertTrue(GPy.kern.kern_test(kern, output_ind=4, verbose=verbose))
 
     def test_ode1_eqkernel(self):
         if SYMPY_AVAILABLE: