code: def spherical_bessel_formulas(n):
      """
      Computes the sympy formulas for the spherical bessel functions up to order n (excluded)
      """
      x = sym.symbols("x")
      # j_i = (-x)^i * (1/x * d/dx)^î * sin(x)/x
      j = [sym.sin(x) / x]  # j_0
      a = sym.sin(x) / x
      for i in range(1, n):
          b = sym.diff(a, x) / x
          j += [sym.simplify(b * (-x) ** i)]
          a = sym.simplify(b)
      return j
  
file_path: DIG/dig/threedgraph/method/pronet/features.py
imports: import torch
  import sympy
  import numpy
  from scipy.optimize import brentq
  from scipy import special
label_desc: Use image node_id 3 for calling a global function with example usage: spherical_bessel_formulas(n) and returns: j