code: def getInterval(means, var, p):
      alpha = abs(normal_CDF_inverse(min(1 - p, p) / 2))
      upper = [0] * len(means)
      lower = [0] * len(means)
      for i in range(0, len(means)):
          upper[i] = means[i] + alpha * math.sqrt(var[i])
          lower[i] = means[i] - alpha * math.sqrt(var[i])
  
      return (upper, lower)
  
file_path: pydlm/pydlm/base/tools.py
imports: import math
label_desc: Use image node_id 4 for calling a global function with example usage: getInterval(means, var, p) and returns: upper,  lower