code: def ceil(x):
      
      if isinstance(x, UncertainFunction):
          mcpts = np.ceil(x._mcpts)
          return UncertainFunction(mcpts)
      else:
          return np.ceil(x)
docs: Ceiling function (round towards positive infinity)