Интегрирование в Mathematica

[Андрей123]

Здравствуйте! Интересует такой вопрос. Мне нужно взять такой вот интеграл

$$\int\limits_0^{y_0} \int\limits_{\frac{x_0}{y_0}\eta}^{x_0-\frac{(b-x_0)(\eta -x_0)}{y_0}} 1/\sqrt{(x_1-\xi)^2+(x_2-\eta)^2+x_3^2}d\xi d\eta$$

. Mathematica такой интеграл берет Судя по всему, при $x_3$ отличном от нуля этот интеграл должен быть дифференцируемой функцией. Но по графику видно, что есть разрывы. Почему это происходит и как можно от этого избавиться? Мои предположения, что это из-за вычисления различных вето комплексного логарифма, который присутствует в окончательном выражении. Вот код.

R[x1_, x2_, x3_] = Sqrt[x1^2 + x2^2 + x3^2];
res = ReplaceAll[
   1/R[x[1], x[2], x[3]], {x[1] -> x1 - \[Xi], x[2] -> x2 - \[Eta], 
    x[3] -> x3}];
res = Integrate[res, \[Xi]];
res = ReplaceAll[res, \[Xi] -> x0/y0 \[Eta]] - 
   ReplaceAll[res, \[Xi] -> x0 - (b - x0)/y0 (\[Eta] - y0)];
res = Integrate[res, \[Eta]];
res = ReplaceAll[res, \[Eta] -> y0] - ReplaceAll[res, \[Eta] -> 0];
\[CapitalPsi][x1_, x2_, x3_, x0_, y0_, b_] = res;
ContourPlot[
 Re[\[CapitalPsi][x, y, 0.01, 1 , 1, 2]], {x, -2, 2}, {y, -2, 2}, 
 PlotRange -> All, ContourLabels -> All, ColorFunction -> "Rainbow"]

[Андрей123]

Вопрос снят.