为什么我在运行这个极坐标下的二重积分时会报错得到非数值呢：
xi = {1, 1, -1, -1}; yi = {-1, 1, 1, -1};
si[i_] :=
NIntegrate[
r Log ((0.2 + r Cos[\[Phi]] - xi[[i]])^2 + (0.2 + r Sin[\[Phi]] -
yi[[i]])^2)^(1/2), {r, 0, 0.5}, {\[Phi], 0, 2 \[Pi]}];
Table[si[i], {i, 1, Length[xi]}]

log后面用方括号