2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Интегрирование в Mathematica
Сообщение31.03.2014, 14:28 
Здравствуйте! Интересует такой вопрос. Мне нужно взять такой вот интеграл
$$\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"]

 
 
 
 Re: Интегрирование в Mathematica
Сообщение02.04.2014, 23:15 
Вопрос снят.

 
 
 [ Сообщений: 2 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group