code: def bode_plot(
      system,
      initial_exp=-5,
      final_exp=5,
      grid=True,
      show_axes=False,
      show=True,
      freq_unit="rad/sec",
      phase_unit="rad",
      phase_unwrap=True,
      **kwargs,
  ):
      r"""
      Returns the Bode phase and magnitude plots of a continuous-time system.
  
      Parameters
      ==========
  
      system : SISOLinearTimeInvariant type
          The LTI SISO system for which the Bode Plot is to be computed.
      initial_exp : Number, optional
          The initial exponent of 10 of the semilog plot. Defaults to -5.
      final_exp : Number, optional
          The final exponent of 10 of the semilog plot. Defaults to 5.
      show : boolean, optional
          If ``True``, the plot will be displayed otherwise
          the equivalent matplotlib ``plot`` object will be returned.
          Defaults to True.
      prec : int, optional
          The decimal point precision for the point coordinate values.
          Defaults to 8.
      grid : boolean, optional
          If ``True``, the plot will have a grid. Defaults to True.
      show_axes : boolean, optional
          If ``True``, the coordinate axes will be shown. Defaults to False.
      freq_unit : string, optional
          User can choose between ``'rad/sec'`` (radians/second) and ``'Hz'`` (Hertz) as frequency units.
      phase_unit : string, optional
          User can choose between ``'rad'`` (radians) and ``'deg'`` (degree) as phase units.
  
      Examples
      ========
  
      .. plot::
          :context: close-figs
          :format: doctest
          :include-source: True
  
          >>> from sympy.abc import s
          >>> from sympy.physics.control.lti import TransferFunction
          >>> from sympy.physics.control.control_plots import bode_plot
          >>> tf1 = TransferFunction(1*s**2 + 0.1*s + 7.5, 1*s**4 + 0.12*s**3 + 9*s**2, s)
          >>> bode_plot(tf1, initial_exp=0.2, final_exp=0.7)   # doctest: +SKIP
  
      See Also
      ========
  
      bode_magnitude_plot, bode_phase_plot
  
      """
      plt.subplot(211)
      mag = bode_magnitude_plot(
          system,
          initial_exp=initial_exp,
          final_exp=final_exp,
          show=False,
          grid=grid,
          show_axes=show_axes,
          freq_unit=freq_unit,
          **kwargs,
      )
      mag.title(f"Bode Plot of ${latex(system)}$", pad=20)
      mag.xlabel(None)
      plt.subplot(212)
      bode_phase_plot(
          system,
          initial_exp=initial_exp,
          final_exp=final_exp,
          show=False,
          grid=grid,
          show_axes=show_axes,
          freq_unit=freq_unit,
          phase_unit=phase_unit,
          phase_unwrap=phase_unwrap,
          **kwargs,
      ).title(None)
  
      if show:
          plt.show()
          return
  
      return plt
  
imports: from sympy.core.numbers import I, pi
  from sympy.functions.elementary.exponential import exp, log
  from sympy.polys.partfrac import apart
  from sympy.core.symbol import Dummy
  from sympy.external import import_module
  from sympy.functions import arg, Abs
  from sympy.integrals.laplace import _fast_inverse_laplace
  from sympy.physics.control.lti import SISOLinearTimeInvariant
  from sympy.plotting.series import LineOver1DRangeSeries
  from sympy.polys.polytools import Poly
  from sympy.printing.latex import latex
label_desc: Use image node_id 14 for calling a global function with example usage: bode_plot(system, initial_exp, final_exp, grid, show_axes, show, freq_unit, phase_unit, phase_unwrap) and returns: plt