code: def main():
      """Plot a geodesic on SE3."""
      n_steps = 40
  
      se3_group = SpecialEuclidean(n=3, point_type="vector")
  
      initial_point = se3_group.identity
      initial_tangent_vec = gs.array([1.8, 0.2, 0.3, 3.0, 3.0, 1.0])
      geodesic = se3_group.metric.geodesic(
          initial_point=initial_point,
          initial_tangent_vec=initial_tangent_vec,
      )
  
      t = gs.linspace(-3.0, 3.0, n_steps)
  
      points = geodesic(t)
  
      visualization.plot(points, space="SE3_GROUP")
      plt.show()
  
imports: import matplotlib.pyplot
  import geomstats.backend
  import geomstats.visualization
  from geomstats.geometry.special_euclidean import SpecialEuclidean
label_desc: Use image node_id 1 for calling a global function with example usage: main() without return types