code: def _expected_jacobian_sphere_immersion(point):
      theta = point[..., 0]
      phi = point[..., 1]
      jacobian = gs.array(
          [
              [
                  gs.cos(phi) * gs.cos(theta),
                  -gs.sin(phi) * gs.sin(theta),
              ],
              [
                  gs.sin(phi) * gs.cos(theta),
                  gs.cos(phi) * gs.sin(theta),
              ],
              [-gs.sin(theta), 0.0],
          ]
      )
      return jacobian
  
imports: import geomstats.backend
  from geomstats.test.data import TestData
  from .riemannian_metric import RiemannianMetricComparisonTestData, RiemannianMetricTestData
label_desc: Use image node_id 2 for calling a global function with example usage: _expected_jacobian_sphere_immersion(point) and returns: jacobian