code: def compute_rdp(q, noise_multiplier, steps, orders):
      """Computes RDP of the Sampled Gaussian Mechanism.
  
      Args:
        q: The sampling rate.
        noise_multiplier: The ratio of the standard deviation of the Gaussian noise
          to the l2-sensitivity of the function to which it is added.
        steps: The number of steps.
        orders: An array (or a scalar) of RDP orders.
  
      Returns:
        The RDPs at all orders. Can be `np.inf`.
      """
      if np.isscalar(orders):
          rdp = _compute_rdp(q, noise_multiplier, orders)
      else:
          rdp = np.array(
              [
                  _compute_rdp(q, noise_multiplier, order)
                  for order in orders
              ]
          )
  
      return rdp * steps
  
imports: from __future__ import absolute_import
  from __future__ import division
  from __future__ import print_function
  import math
  import sys
  import numpy
  from scipy import special
  import six
label_desc: Use image node_id 15 for calling a global function with example usage: compute_rdp(q, noise_multiplier, steps, orders) and returns: unknown