s_max = numpy.array(
                [utils.sigmas(S(s[p]), 'max') for p in range(points)])