function y=legendreQ3(x)
y=2/3-5/2*x.^2-3/4*x.*(1-5/3*x.^2).*log((1+x)./(1-x));
