code: def plot_distribution(X, Y, loc, scale):
      plt.figure(figsize=(15, 5))
      x = X.squeeze()
      for k in (1, 2):
          lb = (loc - k * scale).squeeze()
          ub = (loc + k * scale).squeeze()
          plt.fill_between(
              x, lb, ub, color="silver", alpha=1 - 0.05 * k**3
          )
      plt.plot(x, lb, color="silver")
      plt.plot(x, ub, color="silver")
      plt.plot(X, loc, color="black")
      plt.scatter(X, Y, color="gray", alpha=0.8)
      plt.show()
      plt.close()
  
imports: import matplotlib.pyplot
  import numpy
  import tensorflow
  import tensorflow_probability
  import gpflow
label_desc: Use image node_id 1 for calling a global function with example usage: plot_distribution(X, Y, loc, scale) without return types