fprintf('    j      Gauss-Legendre            error\n')
 fprintf('----------------------------------------------------\n')

 for j=1:15
    eval(sprintf('f=inline(''x.^%d'');',j));
    intf = gausslegendre(f,0,1);
    xact = 1/(j+1);
    fprintf('   %2d   %19.17f   %19.17f  \n', j, intf, xact-intf);
 end