pi*(integrate ((y-b)/(-b^3))^2 dy from y=0..b)