code: def MikotaPair(n):
      # Mikota pair acts as a nice test since the eigenvalues
      # are the squares of the integers n, n=1,2,...
      x = np.arange(1, n + 1)
      B = diag(1.0 / x)
      y = np.arange(n - 1, 0, -1)
      z = np.arange(2 * n - 1, 0, -2)
      A = diag(z) - diag(y, -1) - diag(y, 1)
      return A, B
  
imports: from __future__ import division, print_function, absolute_import
  import itertools
  import numpy
  from numpy.testing import assert_almost_equal, assert_equal, assert_allclose, assert_array_less
  from scipy import ones, rand, r_, diag, linalg, eye
  from scipy.linalg import eig, eigh, toeplitz
  import scipy.sparse
  from scipy.sparse.linalg.eigen.lobpcg import lobpcg
  from scipy.sparse.linalg import eigs
  from scipy.sparse import spdiags
  import pytest
label_desc: Use image node_id 2 for calling a global function with example usage: MikotaPair(n) and returns: A,  B