plot[2/(pi-1)^2*(u^3-2*pi*u^2+pi^2*u), {u, 0, pi}]