code: def test_spherical_jn_yn_cross_product_2(self):
      # https://dlmf.nist.gov/10.50.E3
      n = np.array([1, 5, 8])
      x = np.array([0.1, 1, 10])
      left = spherical_jn(n + 2, x) * spherical_yn(n, x) - spherical_jn(
          n, x
      ) * spherical_yn(n + 2, x)
      right = (2 * n + 3) / x**3
      assert_allclose(left, right)
  
imports: import numpy
  from numpy.testing import assert_almost_equal, assert_allclose, assert_array_almost_equal, suppress_warnings
  import pytest
  from numpy import sin, cos, sinh, cosh, exp, inf, nan, r_, pi
  from scipy.special import spherical_jn, spherical_yn, spherical_in, spherical_kn
  from scipy.integrate import quad
label_desc: Use image node_id 2 for calling the TestSphericalJnYnCrossProduct obj's underlying member method code with example usage: obj.test_spherical_jn_yn_cross_product_2() without return types