How to use the quadpy.t2._helpers.s3 function in quadpy

def lyness_jespersen_17():
    weights, points = concat(
        s3(-2.834183851113958e-01),
        s2(
            [2.097208857979572e-01 / 3, 4.766654393821525e-01],
            [5.127273801480265e-02 / 3, 3.377184405448033e-02],
            [6.564896469913508e-01 / 3, 2.703478891654040e-01],
        ),
        s1([3.659351143072855e-01 / 6, 5.146433548666149e-02, 7.458294907672514e-01]),
    )
    return T2Scheme("Lyness-Jespersen 17", weights, points, 8, source)

def liu_vinokur_05():
    weights, points = concat(
        s3(-frac(9, 16)),
        # ERR Incorrectly specified in the article as 25 (instead of 2/5).
        alpha([frac(25, 48), frac(2, 5)]),
    )
    return T2Scheme("Liu-Vinokur 5", weights, points, 3, source)

def laursen_gellert_13():
    weights, points = concat(
        s3(0.097135796282799),
        s2(
            [0.031334700227139, 0.489682519198738],
            [0.077827541004774, 0.437089591492937],
            [0.079647738927210, 0.188203535619033],
            [0.025577675658698, 0.044729513394453],
        ),
        s1([0.043283539377289, 0.036838412054736, 0.221962989160766]),
    )
    return T2Scheme("Laursen-Gellert 13", weights, points, 9, source)

def williams_shunn_jameson_7():
    weights, points = concat(
        s3(0.083608212215637),
        s2(
            [0.005272170280495, 0.019977187122193],
            [0.044552936679504, 0.131721767529998],
            [0.033815712804198, 0.485135346793461],
        ),
        s1(
            [0.015710461340183, 0.107951981846011, 0.867911210117951],
            [0.028205136280616, 0.270840772921567, 0.700872570380723],
            [0.066995957127830, 0.316549598844617, 0.536654684206138],
        ),
    )
    return T2Scheme("Williams-Shunn-Jameson 7", weights, points, 10, source, 4.041e-13)

def liu_vinokur_01():
    weights, points = s3(1)
    return T2Scheme("Liu-Vinokur 1", weights, points, 1, source)

def laursen_gellert_15a():
    weights, points = concat(
        s3(0.079894504741240),
        s2(
            [0.071123802232377, 0.425086210602091],
            [0.008223818690464, 0.023308867510000],
        ),
        s1(
            [0.045430592296170, 0.147925626209534, 0.223766973576973],
            [0.037359856234305, 0.029946031954171, 0.358740141864431],
            [0.030886656884564, 0.035632559587504, 0.143295370426867],
        ),
    )
    return T2Scheme("Laursen-Gellert 15a", weights, points, 10, source)

def laursen_gellert_15b():
    weights, points = concat(
        s3(0.081743329146286),
        s2(
            [0.045957963604745, 0.142161101056564],
            [0.013352968813150, 0.032055373216944],
        ),
        s1(
            [0.063904906396424, 0.148132885783821, 0.321812995288835],
            [0.034184648162959, 0.029619889488730, 0.369146781827811],
            [0.025297757707288, 0.028367665339938, 0.163701733737182],
        ),
    )
    return T2Scheme("Laursen-Gellert 15b", weights, points, 10, source)

def laursen_gellert_12():
    weights, points = concat(
        s3(0.144315607677787),
        s2(
            [0.103217370534718, 0.170569307751761],
            [0.032458497623198, 0.050547228317031],
            [0.095091634267284, 0.459292588292723],
        ),
        s1([0.027230314174435, 0.008394777409958, 0.263112829634638]),
    )
    return T2Scheme("Laursen-Gellert 12", weights, points, 8, source)

def liu_vinokur_12():
    sqrt15 = sqrt(15)
    weights, points = concat(
        s3(frac(9, 40)),
        alpha(
            [(155 - sqrt15) / 1200, (1 + sqrt15) / 7],
            [(155 + sqrt15) / 1200, (1 - sqrt15) / 7],
        ),
    )
    return T2Scheme("Liu-Vinokur 12", weights, points, 5, source)

def dunavant_01():
    weights, points = s3(1)
    return T2Scheme("Dunavant 1", weights, points, 1, source, 7.850e-17)