численное решение интеграла от экспоненты в Mathematica

[Leierkastenmann]

1. Куда пропало определение t?
2. Так и не избавились от P(t).

c = 3*10^8;
H = 7.2*10^5;
\[Tau] = 47*10^-6;
Fd = 7200;
Tp = 50*10^-6;
f = 3.5*10^8;
v = 7500;
Q = 18000;
L = 100;
s2 = 0.003;
s1 = 2;
\[CapitalTheta]0 = 0.018;
m = 256;
\[Lambda] = 0.022;
n = 3;
lm = 25;
d = 1.2;
\[Lambda] = 0.022;
n = 3;
t = List[ -10^-8, -10^-9, 0, 10^-16, 10^-13, 10^-9, 10^-8, 20*10^-8];
p = Sum[NIntegrate[
   Exp[x^2*(\[Pi]*Fd^2*\[Tau]^2*(m - Abs[k])^2)/
       H^2 - (0.00005*n^2*\[Pi]^2*d^2*
        y^2)/(\[CapitalTheta]0^2*\[Lambda]^2*H^2) - (100*lm^2*x^2*
        s1^2)/(H^2*L^2*\[Pi]^2*s2^2) - (100*lm^2*y^2*s1^2)/(H^2*
        L^2*\[Pi]^2*
        s2^2) - (5.55*(y^2 + 
          x^2))/(H^2*\[CapitalTheta]0^2) - \[Pi]*(f^2*(t - k*Tp - 
            x^2/(c*H) - y^2/(c*H))^2 + 
        2*Q*(v*k*Tp)/(H*
            d) + \[Tau]^2*((v*k*Tp)/(H*d))^2)], {x, -Infinity, 
    Infinity}, {y, -Infinity, Infinity}, MaxRecursion -> 40, 
   AccuracyGoal -> 60, Method -> "AdaptiveMonteCarlo"], {k, -255, 
   255}]
ListPlot[Transpose[{t, p}]]

У меня этот код работает без ошибок. Правда с таким масштабом по t и p лучше все же, наверное, строить логарифмические графики. И еще одно замечание - я-то привел метод Монте-Карло как тот, который проинтегрирует без ошибок, но с ним надо быть аккуратным, он может для одной и той же функции выдавать различные результаты.

[salang]

Спасибо, получилось, но гарфик корявый.
Как сделать плавную, т.е. непрерывную кривую на графике?
Вообще возможно задание построение графиков вида y=x(t)?

[Leierkastenmann]

Ну какая дискретность, такой в общем-то и график. Сделайте разбиение по t чаще, получите более плавную кривую (наверное). Графики функции конечно же строить можно, например, если переделать в функцию

pp[tt_] := 
 Sum[NIntegrate[
   Exp[x^2*(\[Pi]*Fd^2*\[Tau]^2*(m - Abs[k])^2)/
       H^2 - (0.00005*n^2*\[Pi]^2*d^2*
        y^2)/(\[CapitalTheta]0^2*\[Lambda]^2*H^2) - (100*lm^2*x^2*
        s1^2)/(H^2*L^2*\[Pi]^2*s2^2) - (100*lm^2*y^2*s1^2)/(H^2*
        L^2*\[Pi]^2*
        s2^2) - (5.55*(y^2 + 
          x^2))/(H^2*\[CapitalTheta]0^2) - \[Pi]*(f^2*(tt - k*Tp - 
            x^2/(c*H) - y^2/(c*H))^2 + 
        2*Q*(v*k*Tp)/(H*
            d) + \[Tau]^2*((v*k*Tp)/(H*d))^2)], {x, -Infinity, 
    Infinity}, {y, -Infinity, Infinity}, MaxRecursion -> 40, 
   AccuracyGoal -> 60, Method -> "AdaptiveMonteCarlo"], {k, -255, 
   255}]
Plot[pp[t], {t, -10^(-8), 20*10^(-8)}]

Но для такой функции график строиться будет так долго, что не факт, что и дождетесь результата.

[salang]

а как сделать логарифмический масштаб по вертикали?

[Aritaborian]

LogPlot.

[salang]

Leierkastenmann в сообщении #807641 писал(а):

график строиться будет так долго

как-то можно задать шаг аргумента для графика?

[Aritaborian]

Опцией PlotPoints.

[salang]

это правильно, что в начале объявлено pp[tt_], а построение идет Plot[pp[t] ?

[Aritaborian]

А что вас смущает?

[salang]

отличие tt_ от t

[Aritaborian]

salang, вы, я смотрю, не очень сильны в Mathematica.
Объясняю. Вот мы определяем функцию: f[x_] := x^2.
А вот мы её вычисляем: f[3 a]. Как вы думаете, что получим? Правильно, 9 а^2 ;-) Если в определении стоит, скажем, икс, это не значит, что и в вычислении должен фигурировать икс. Можете, впрочем, в своём коде заменить tt на t или наоборот, всё по-прежнему будет работать (или не работать ;-)
А знак подчёркивания _ в определении функции как раз указывает, что на месте икса может быть что угодно.
Загляните в эти разделы справки: Defining Functions, Functions as Procedures, How to Work with Variables and Functions.

[salang]

PlotPoints задает число выборок, а есть возможность задать абсолютный шаг?
В ListLinePlot не получается вставить PlotPoints. Это ошибка?
Cколько примерно для моих данных занимает расчет непрерывной кривой на средней машине?

c = 3*10^8;
H = 7.2*10^5;
\[Tau] = 47*10^-6;
Fd = 7200;
Tp = 50*10^-6;
f = 3.5*10^8;
v = 7500;
Q = 18000;
L = 100;
s2 = 0.003;
s1 = 2;
\[CapitalTheta]0 = 0.018;
m = 256;
\[Lambda] = 0.022;
n = 3;
lm = 25;
d = 1.2;
\[Lambda] = 0.022;
n = 3;
t = List[ -10^-8, -10^-9, 0, 10^-11, 10^-10, 10^-9, 10^-8, 20*10^-7];

pp[tt_] = Sum[
  NIntegrate [
   Exp[x^2*(\[Pi]*Fd^2*\[Tau]^2*(m - Abs[k])^2 )/H^2 - (
     0.00005*n^2*\[Pi]^2*d^2*y^2 )/(\[CapitalTheta]0^2*\[Lambda]^2 *H^2) - (
     100*lm^2*x^2*s1^2)/(H^2*L^2*\[Pi]^2*s2^2) - (100*lm^2*y^2*s1^2)/(
     H^2*L^2*\[Pi]^2*s2^2) - (5.55*(y^2 + x^2))/(
     H^2*\[CapitalTheta]0^2) - \[Pi] * (f^2 *(t - k*Tp - x^2/(c*H) - y^2/(
           c*H))^2 + 
        2*Q*(v*k*Tp)/(H*d) + \[Tau]^2*((v*k*Tp)/(H*d))^2 )], {x, -Infinity, 
    Infinity}, {y, -Infinity, Infinity}, MaxRecursion -> 40, 
   AccuracyGoal -> 60, Method -> "AdaptiveMonteCarlo"], {k, -255, 255}]


In[1]:= Plot[pp[t], {t, -10^(-8), 20*10^(-8)}, PlotPoints -> 20]
ListLogPlotPlot[Transpose[{t, p}]]

Out[1]= \!\(\*
GraphicsBox[{},
AspectRatio->0.6180339887498948,
Axes->True,
AxesLabel->{None, None},
AxesOrigin->{0, 0},
Method->{},
PlotRange->{{-1.*^-8, 2.*^-7}, {0., 0.}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02], 
Scaled[0.02]}]\)

During evaluation of In[1]:= Transpose::nmtx: The first two levels of the one-dimensional list {t,p} cannot be transposed. >>

Out[2]= ListLogPlotPlot[Transpose[{t, p}]]

Sum[NIntegrate[
  Exp[x^2*(\[Pi]*Fd^2*\[Tau]^2*(m - Abs[k])^2)/
      H^2 - (0.00005*n^2*\[Pi]^2*d^2*y^2)/(\[CapitalTheta]0^2*\[Lambda]^2*
       H^2) - (100*lm^2*x^2*s1^2)/(H^2*L^2*\[Pi]^2*s2^2) - (100*lm^2*y^2*
       s1^2)/(H^2*L^2*\[Pi]^2*
       s2^2) - (5.55*(y^2 + 
         x^2))/(H^2*\[CapitalTheta]0^2) - \[Pi]*(f^2*(t - k*Tp - x^2/(c*H) - 
           y^2/(c*H))^2 + 
       2*Q*(v*k*Tp)/(H*d) + \[Tau]^2*((v*k*Tp)/(H*d))^2)], {x, -Infinity, 
   Infinity}, {y, -Infinity, Infinity}, MaxRecursion -> 40, 
  AccuracyGoal -> 60, Method -> "AdaptiveMonteCarlo"], {k, -255, 255}]

{1.5983*10^-12, 335063., 836068., 936968., 90404., 0., 0., 702430.}

ListLinePlot[%, AxesLabel -> {"t", "P(t)"}, PlotStyle -> PointSize[0.01]]

\!\(\*
GraphicsBox[{{}, {{}, {}, 
{RGBColor[0.24720000000000017`, 0.24, 0.6], PointSize[0.01], 
      LineBox[{{1., 1.598299392712959*^-12}, {2., 335062.66072612564`}, {3., 
       836067.9196629806}, {4., 936968.342828458}, {5., 90403.98299123366}, {
       6., 0.}, {7., 0.}, {8., 702429.7247402478}}]}}, {}},
AspectRatio->0.6180339887498948,
Axes->True,
AxesLabel->{
FormBox["\"t\"", TraditionalForm], 
FormBox["\"P(t)\"", TraditionalForm]},
AxesOrigin->{0, 0},
Method->{},
PlotRange->{{0, 8.}, {0, 936968.342828458}},
PlotRangeClipping->True,
PlotRangePadding->{{0.16, 0.16}, {18739.36685656916, 18739.36685656916}}]\)

t = List[ -10^-8, -10^-9, 0, 10^-9, 10^-8, 10^-7, 10^-6, 20*10^-5];

Файл: http://zalil.ru/34881989/1310a13d.52c89978/6.nb

[Aritaborian]

salang в сообщении #809443 писал(а):

есть возможность задать абсолютный шаг?

Используйте дополнительно опцию MaxRecursion -> 0.

salang в сообщении #809443 писал(а):

В ListLinePlot не получается вставить PlotPoints. Это ошибка?

Разумеется, ошибка. ListLinePlot строит график по уже имеющемуся набору точек.

salang в сообщении #809443 писал(а):

Сколько примерно для моих данных занимает расчет непрерывной кривой на средней машине?

А что такое средняя машина?

[salang]

MaxRecursion вставляет без стрелки-указателя
Под средней машиной я имел в виду обычный Prescott 2,8+4Гб RAM

c = 3*10^8;
H = 7.2*10^5;
\[Tau] = 47*10^-6;
Fd = 7200;
Tp = 50*10^-6;
f = 3.5*10^8;
v = 7500;
Q = 18000;
L = 100;
s2 = 0.003;
s1 = 2;
\[CapitalTheta]0 = 0.018;
m = 256;
\[Lambda] = 0.022;
n = 3;
lm = 25;
d = 1.2;
\[Lambda] = 0.022;
n = 3;
t = List[ -10^-8, -10^-9, 0, 10^-11, 10^-10, 10^-9, 10^-8, 20*10^-7];

In[1]:= pp[tt_] = Sum[
  NIntegrate [
   Exp[x^2*(\[Pi]*Fd^2*\[Tau]^2*(m - Abs[k])^2 )/H^2 - (
     0.00005*n^2*\[Pi]^2*d^2*y^2 )/(\[CapitalTheta]0^2*\[Lambda]^2 *H^2) - (
     100*lm^2*x^2*s1^2)/(H^2*L^2*\[Pi]^2*s2^2) - (100*lm^2*y^2*s1^2)/(
     H^2*L^2*\[Pi]^2*s2^2) - (5.55*(y^2 + x^2))/(
     H^2*\[CapitalTheta]0^2) - \[Pi] * (f^2 *(t - k*Tp - x^2/(c*H) - y^2/(
           c*H))^2 + 
        2*Q*(v*k*Tp)/(H*d) + \[Tau]^2*((v*k*Tp)/(H*d))^2 )], {x, -Infinity, 
    Infinity}, {y, -Infinity, Infinity}, MaxRecursion -> 40, 
   AccuracyGoal -> 60, Method -> "AdaptiveMonteCarlo"], {k, -255, 255}]


During evaluation of In[1]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-<<1>>/<<1>>-<<1>>/<<1>>+(Fd^2 <<3>> <<1>>)/H^2-\[Pi] (-((510 Q Tp v)/(d H))+f^2 (t+Times[<<2>>]+Times[<<4>>]+Times[<<4>>])^2+(65025 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[1]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-<<1>>/<<1>>-<<1>>/<<1>>+(Fd^2 <<3>> <<1>>)/H^2-\[Pi] (-((508 Q Tp v)/(d H))+f^2 (t+Times[<<2>>]+Times[<<4>>]+Times[<<4>>])^2+(64516 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[1]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-<<1>>/<<1>>-<<1>>/<<1>>+(Fd^2 <<3>> <<1>>)/H^2-\[Pi] (-((506 Q Tp v)/(d H))+f^2 (t+Times[<<2>>]+Times[<<4>>]+Times[<<4>>])^2+(64009 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[1]:= General::stop: Further output of NIntegrate::inumr will be suppressed during this calculation. >>

Out[1]= 511 NIntegrate[
  Exp[(x^2 (\[Pi] Fd^2 \[Tau]^2 (m - Abs[k])^2))/H^2 - (
    0.00005 n^2 \[Pi]^2 d^2 y^2)/(\[CapitalTheta]0^2 \[Lambda]^2 H^2) - (
    100 lm^2 x^2 s1^2)/(H^2 L^2 \[Pi]^2 s2^2) - (100 lm^2 y^2 s1^2)/(
    H^2 L^2 \[Pi]^2 s2^2) - (5.55 (y^2 + x^2))/(
    H^2 \[CapitalTheta]0^2) - \[Pi] (f^2 (t - k Tp - x^2/(c H) - y^2/(
          c H))^2 + (2 Q (v k Tp))/(
       H d) + \[Tau]^2 ((v k Tp)/(
         H d))^2)], {x, -\[Infinity], \[Infinity]}, {y, -\[Infinity], \
\[Infinity]}, MaxRecursion -> 40, AccuracyGoal -> 60, 
  Method -> "AdaptiveMonteCarlo"]

In[2]:= Plot[pp[t], {t, -10^(-8), 20*10^(-8)}, PlotPoints -> 20]
ListLogPlotPlot[Transpose[{t, p}]]

During evaluation of In[2]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-(5.55 (<<1>>))/(H^2 \[CapitalTheta]0^2)-<<1>>/<<1>>-\[Pi] ((2 k Q Tp v)/(d H)+f^2 (<<1>>)^2+(k^2 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))+(Fd^2 \[Pi] x^2 \[Tau]^2 (m-Abs[k])^2)/H^2) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[2]:= NIntegrate::inumr: The integrand E^(<<9>>+(3.14159 Fd^2 x^2 \[Tau]^2 (m-1. <<1>>)^2)/H^2) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[2]:= NIntegrate::inumr: The integrand E^(<<9>>+(3.14159 Fd^2 x^2 \[Tau]^2 (m-1. <<1>>)^2)/H^2) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[2]:= General::stop: Further output of NIntegrate::inumr will be suppressed during this calculation. >>

Out[2]= \!\(\*
GraphicsBox[{},
AspectRatio->0.6180339887498948,
Axes->True,
AxesLabel->{None, None},
AxesOrigin->{0, 0},
Method->{},
PlotRange->{{-1.*^-8, 2.*^-7}, {0., 0.}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02], 
Scaled[0.02]}]\)

During evaluation of In[2]:= Transpose::nmtx: The first two levels of the one-dimensional list {t,p} cannot be transposed. >>

ListLogPlotPlot[Transpose[{t, p}], MaxRecursion]

In[4]:= p[tt_] = Sum[
  NIntegrate[
   Exp[x^2*(\[Pi]*Fd^2*\[Tau]^2*(m - Abs[k])^2)/
       H^2 - (0.00005*n^2*\[Pi]^2*d^2*y^2)/(\[CapitalTheta]0^2*\[Lambda]^2*
        H^2) - (100*lm^2*x^2*s1^2)/(H^2*L^2*\[Pi]^2*s2^2) - (100*lm^2*y^2*
        s1^2)/(H^2*L^2*\[Pi]^2*
        s2^2) - (5.55*(y^2 + 
          x^2))/(H^2*\[CapitalTheta]0^2) - \[Pi]*(f^2*(t - k*Tp - x^2/(c*H) - 
            y^2/(c*H))^2 + 
        2*Q*(v*k*Tp)/(H*d) + \[Tau]^2*((v*k*Tp)/(H*d))^2)], {x, -Infinity, 
    Infinity}, {y, -Infinity, Infinity}, MaxRecursion -> 40, 
   AccuracyGoal -> 60, Method -> "AdaptiveMonteCarlo"], {k, -255, 255}]

During evaluation of In[4]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-<<1>>/<<1>>-<<1>>/<<1>>+(Fd^2 <<3>> <<1>>)/H^2-\[Pi] (-((510 Q Tp v)/(d H))+f^2 (t+Times[<<2>>]+Times[<<4>>]+Times[<<4>>])^2+(65025 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[4]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-<<1>>/<<1>>-<<1>>/<<1>>+(Fd^2 <<3>> <<1>>)/H^2-\[Pi] (-((508 Q Tp v)/(d H))+f^2 (t+Times[<<2>>]+Times[<<4>>]+Times[<<4>>])^2+(64516 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[4]:= NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-<<1>>/<<1>>-<<1>>/<<1>>+(Fd^2 <<3>> <<1>>)/H^2-\[Pi] (-((506 Q Tp v)/(d H))+f^2 (t+Times[<<2>>]+Times[<<4>>]+Times[<<4>>])^2+(64009 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

During evaluation of In[4]:= General::stop: Further output of NIntegrate::inumr will be suppressed during this calculation. >>

Out[4]= 511 NIntegrate[
  Exp[(x^2 (\[Pi] Fd^2 \[Tau]^2 (m - Abs[k])^2))/H^2 - (
    0.00005 n^2 \[Pi]^2 d^2 y^2)/(\[CapitalTheta]0^2 \[Lambda]^2 H^2) - (
    100 lm^2 x^2 s1^2)/(H^2 L^2 \[Pi]^2 s2^2) - (100 lm^2 y^2 s1^2)/(
    H^2 L^2 \[Pi]^2 s2^2) - (5.55 (y^2 + x^2))/(
    H^2 \[CapitalTheta]0^2) - \[Pi] (f^2 (t - k Tp - x^2/(c H) - y^2/(
          c H))^2 + (2 Q (v k Tp))/(
       H d) + \[Tau]^2 ((v k Tp)/(
         H d))^2)], {x, -\[Infinity], \[Infinity]}, {y, -\[Infinity], \
\[Infinity]}, MaxRecursion -> 40, AccuracyGoal -> 60, 
  Method -> "AdaptiveMonteCarlo"]

ListLinePlot[%, AxesLabel -> {"t", "P(t)"}, PlotStyle -> PointSize[0.01]]

NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-(5.55 (<<1>>))/(H^2 \[CapitalTheta]0^2)-<<1>>/<<1>>-\[Pi] ((2 k Q Tp v)/(d H)+f^2 (t+<<3>>)^2+(k^2 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))+(Fd^2 \[Pi] x^2 \[Tau]^2 (m-Abs[k])^2)/H^2) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-(5.55 (<<1>>))/(H^2 \[CapitalTheta]0^2)-<<1>>/<<1>>-\[Pi] ((2 k Q Tp v)/(d H)+f^2 (t+<<3>>)^2+(k^2 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))+(Fd^2 \[Pi] x^2 \[Tau]^2 (m-Abs[k])^2)/H^2) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

NIntegrate::inumr: The integrand E^(-((100 lm^2 s1^2 x^2)/(H^2 L^2 \[Pi]^2 s2^2))-(100 lm^2 s1^2 y^2)/(H^2 L^2 \[Pi]^2 s2^2)-(5.55 (<<1>>))/(H^2 \[CapitalTheta]0^2)-<<1>>/<<1>>-\[Pi] ((2 k Q Tp v)/(d H)+f^2 (t+<<3>>)^2+(k^2 Tp^2 v^2 \[Tau]^2)/(d^2 H^2))+(Fd^2 \[Pi] x^2 \[Tau]^2 (m-Abs[k])^2)/H^2) has evaluated to non-numerical values for all sampling points in the region with boundaries {{0,1},{0,1}}. >>

General::stop: Further output of NIntegrate::inumr will be suppressed during this calculation. >>

ListLinePlot::lpn: 511 NIntegrate[Exp[(x^2 (\[Pi] Power[<<2>>] Power[<<2>>] Power[<<2>>]))/Power[<<2>>]-(0.00005 Power[<<2>>] Power[<<2>>] Power[<<2>>] Power[<<2>>])/Times[<<3>>]-(100 Power[<<2>>] Power[<<2>>] Power[<<2>>])/Times[<<4>>]-(100 Power[<<2>>] Power[<<2>>] Power[<<2>>])/Times[<<4>>]-(5.55 Plus[<<2>>])/Times[<<2>>]-\[Pi] (Times[<<2>>]+Times[<<3>>]+Times[<<2>>])],{x,-\[Infinity],\[Infinity]},<<3>>,Method->AdaptiveMonteCarlo] is not a list of numbers or pairs of numbers. >>

[Aritaborian]

salang в сообщении #809852 писал(а):

MaxRecursion вставляет без стрелки-указателя

Простите, не понял вас.
И прячьте код под спойлеры, что ли, а то как-то раздражает листать эти простыни вверх-вниз.