code: def compute_rdp_of_smooth_sensitivity_gaussian(beta, sigma, order):
      """Computes the RDP curve for the GNSS mechanism.
  
      Implements Theorem 23 (https://arxiv.org/pdf/1802.08908.pdf).
      """
      if beta > 0 and not 1 < order < 1 / (2 * beta):
          raise ValueError("Order outside the (1, 1/(2*beta)) range.")
  
      return order * math.exp(2 * beta) / sigma**2 + (
          -0.5 * math.log(1 - 2 * order * beta) + beta * order
      ) / (order - 1)
  
imports: from __future__ import absolute_import
  from __future__ import division
  from __future__ import print_function
  import math
  from absl import app
  import numpy
  import scipy
  import sympy
  import core
label_desc: Use image node_id 15 for calling a global function with example usage: compute_rdp_of_smooth_sensitivity_gaussian(beta, sigma, order) and returns: unknown