Creating a parameter with the arcsine distribution provides incorrect weights (nodes are correct).
cheb = eq.Parameter(distribution='arcsine', order=9, lower=-1, upper=1, endpoints='both')
points, weights = cheb._get_local_quadrature()
1 Like
Working on this…will have a fix soon 