def cdf(self, x): return 0.5 * (1. + erf.erf(( x - broadcast.broadcast_to(self.loc, x.shape)) / (2 ** 0.5 * broadcast.broadcast_to(self.scale, x.shape))))
def survival_function(self, x): return 0.5 * (1. - erf.erf( (x - broadcast.broadcast_to(self.loc, x.shape)) / (2 ** 0.5 * broadcast.broadcast_to(self.scale, x.shape))))