% Esimerkkej? Gaussin-Legendren kvadratuurista
n=10000;
k=1:n;

beta=sqrt(k.^2 ./ (4 * k.^2 - 1));
alpha=zeros(size(beta)+[0,1]);

[x,w]=gausslegendre(alpha,beta);

integral(@(x) x.^2,-1,1)
integral(@(x) x.^5,-1,1)

w*x.^2
w*x.^5

integral(@(x) exp(sin(pi * x)),-1,1)
w*exp(sin(pi*x))