plot (x(pi - acos(x)) + sqrt(1-x^2)*(x^2+2)/3)/(pi (x+1)), from -1 to 1