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

def laursen_gellert_11():
    weights, points = concat(
        s2([0.053077801790233, 0.064930513159165]),
        s1(
            [0.070853083692136, 0.284575584249173, 0.517039939069325],
            [0.069274682079415, 0.313559184384932, 0.043863471792371],
        ),
    )
    return T2Scheme("Laursen-Gellert 11", weights, points, 7, source, 1.306e-14)

def dunavant_19():
    weights, points = concat(
        s3(0.032906331388919),
        s2(
            [0.010330731891272, 0.489609987073006],
            [0.022387247263016, 0.454536892697893],
            [0.030266125869468, 0.401416680649431],
            [0.030490967802198, 0.255551654403098],
            [0.024159212741641, 0.177077942152130],
            [0.016050803586801, 0.110061053227952],
            [0.008084580261784, 0.055528624251840],
            [0.002079362027485, 0.012621863777229],
        ),
        s1(
            [0.003884876904981, 0.395754787356943, 0.600633794794645],
            [0.025574160612022, 0.307929983880436, 0.557603261588784],
            [0.008880903573338, 0.264566948406520, 0.720987025817365],
            [0.016124546761731, 0.358539352205951, 0.594527068955871],
            [0.002491941817491, 0.157807405968595, 0.839331473680839],
            [0.018242840118951, 0.075050596975911, 0.701087978926173],
            [0.010258563736199, 0.142421601113383, 0.822931324069857],
            [0.003799928855302, 0.065494628082938, 0.924344252620784],
        ),
    )
    return T2Scheme("Dunavant 19", weights, points, 19, source, 1.670e-14)

def williams_shunn_jameson_5():
    weights, points = concat(
        s2(
            [0.017915455012303, 0.035870877695734],
            [0.127712195881265, 0.241729395767967],
            [0.076206062385535, 0.474308787777079],
        ),
        s1([0.055749810027115, 0.201503881881800, 0.751183631106484]),
    )
    return T2Scheme("Williams-Shunn-Jameson 5", weights, points, 7, source, 5.652e-14)

def lyness_jespersen_06():
    a0, a1 = [(3 + i * sqrt(3)) / 6 for i in [+1, -1]]

    weights, points = concat(
        s3(+frac(9, 20)), s2([-frac(1, 60), 0]), s1([+frac(1, 10), a0, a1])
    )
    return T2Scheme("Lyness-Jespersen 6", weights, points, 4, source)

def lyness_jespersen_19():
    weights, points = concat(
        s3(1.133624844599192e-01),
        s2(
            [1.062573789846330e-03 / 3, 0],
            [4.803411513859279e-02 / 3, 0.5],
            [2.524243006337300e-01 / 3, 4.497793381870162e-01],
            [7.819254371487040e-02 / 3, 4.694744319909033e-02],
            [2.472227459993048e-01 / 3, 1.918719127374489e-01],
        ),
        s1([2.597012362637364e-01 / 6, 3.683841205473626e-02, 7.411985987844980e-01]),
    )
    return T2Scheme("Lyness-Jespersen 19", weights, points, 9, source)

def lyness_jespersen_10():
    weights, points = concat(
        s2(
            [3.503588271790222e-01 / 3, 2.492867451709329e-01],
            [1.525347191106164e-01 / 3, 6.308901449150177e-02],
        ),
        s1([4.971064537103375e-01 / 6, 6.365024991213939e-01, 5.314504984483216e-02]),
    )
    return T2Scheme("Lyness-Jespersen 10", weights, points, 6, source, 1.240e-13)

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 laursen_gellert_04():
    roots = numpy.polynomial.polynomial.polyroots([-1, 15, -60, 60])
    weights, points = s1([1.0 / 6.0, roots[2], roots[1]])
    return T2Scheme("Laursen-Gellert 4", weights, points, 3, source)

def lyness_jespersen_15():
    weights, points = concat(
        s3(1.443156076777862e-01),
        s2(
            [2.852749028018549e-01 / 3, 4.592925882927229e-01],
            [9.737549286959440e-02 / 3, 5.054722831703103e-02],
            [3.096521116041552e-01 / 3, 1.705693077517601e-01],
        ),
        s1([1.633818850466092e-01 / 6, 8.394777409957211e-03, 7.284923929554041e-01]),
    )
    return T2Scheme("Lyness-Jespersen 15", weights, points, 8, source)

def dunavant_18():
    weights, points = concat(
        s3(+0.030809939937647),
        s2(
            [+0.009072436679404, 0.493344808630921],
            [+0.018761316939594, 0.469210594241957],
            [+0.019441097985477, 0.436281395887006],
            [+0.027753948610810, 0.394846170673416],
            [+0.032256225351457, 0.249794568803157],
            [+0.025074032616922, 0.161432193743843],
            [+0.015271927971832, 0.076598227485371],
            [+0.006793922022963, 0.024252439353450],
            [-0.002223098729920, 0.043146367216965],
        ),
        s1(
            [+0.006331914076406, 0.358911494940944, 0.632657968856636],
            [+0.027257538049138, 0.294402476751957, 0.574410971510855],
            [+0.017676785649465, 0.325017801641814, 0.624779046792512],
            [+0.018379484638070, 0.184737559666046, 0.748933176523037],
            [+0.008104732808192, 0.218796800013321, 0.769207005420443],
            [+0.007634129070725, 0.101179597136408, 0.883962302273467],
            [+0.000046187660794, 0.020874755282586, 1.014347260005363],
        ),
    )
    return T2Scheme("Dunavant 18", weights, points, 18, source, 5.089e-13)