code: def Jn_zeros(n, k):
      zerosj = np.zeros((n, k), dtype="float32")
      zerosj[0] = np.arange(1, k + 1) * np.pi
      points = np.arange(1, k + n) * np.pi
      racines = np.zeros(k + n - 1, dtype="float32")
      for i in range(1, n):
          for j in range(k + n - 1 - i):
              foo = brentq(Jn, points[j], points[j + 1], (i,))
              racines[j] = foo
          points = racines
          zerosj[i][:k] = racines[:k]
  
      return zerosj
  
imports: import numpy
  from scipy.optimize import brentq
  from scipy import special
  import torch
  from math import sqrt, pi
label_desc: Use image node_id 2 for calling a global function with example usage: Jn_zeros(n, k) and returns: zerosj