restart;
approx:=f->sum(w[k]*f(x[k]),k=1..n);
exact:=f->int(f(x),x=-1..1);
approx(t->t^2);
exact(t->t^2);
n:=5;
for j from 0 to n-1
do
    eq[j]:=approx(t->t^j)=exact(t->t^j);
end do;
eqs:={seq(eq[k],k=0..n-1)};
h:=2/(n-1);
for j from 1 to n
do
    x[j]:=-1+h*(j-1);
end do;
eqs;
vbls:={seq(w[k],k=1..n)};
solve(eqs,vbls);
[seq(x[k],k=1..n)];