Help:=proc(): print(`GaussQ(d) `):end: #GaussQ(d):conjectures Gaussian Quadrature #that integrates exactly polynomials of #degree d GaussQ:=proc(d) local x,A,var,eq,sol: var:={seq(A[i],i=1..d), seq(x[i],i=1..d)}: eq:= { seq( 1/(i+1)=add(A[j]*x[j]^i,j=1..d), i=0..2*d-1)}: sol:=solve(eq,var): subs(sol,[[seq(x[i],i=1..d)], [seq(A[i],i=1..d)]]): end: