Bài học 9 bước, chạy theo thứ tự trong một phiên. Từng bước đều được thực thi trên bộ máy và đầu ra bên dưới là những gì nó tạo ra.

Bước 1

NIntegrate[Sin[200 x ^ 2 + 5](1/(x + 1) ^ 2), {x, 0, ∞}, Method -> "ExtrapolatingOscillatory"]

Đầu ra

NIntegrate[Sin[5 + 200*x^2]/(1 + x)^2, {x, 0, Infinity}, Method -> ExtrapolatingOscillatory]

Bước 2

Clear[ω, k, f];
ω = 20;
k[x_] := Sin[ω x];
f[x_] := (1/(x + 1) ^ 2);

Không có đầu ra - bước này chuẩn bị cho bước kế tiếp.

Bước 3

Plot[k[x]f[x], {x, 0, 10}, PlotPoints -> 1000, PlotRange -> All]

Đầu ra

-Graphics-

Bước 4

psum[i_ ? NumberQ] := NIntegrate[k[x]f[x], {x, i(1/ω) π, (i + 1)(1/ω) π}]

Không có đầu ra - bước này chuẩn bị cho bước kế tiếp.

Bước 5

res = NSum[psum[i], {i, 0, ∞}, Method -> "AlternatingSigns", "VerifyConvergence" -> False]

Đầu ra

0.04928412772506762

Bước 6

exact = Integrate[k[x]f[x], {x, 0, ∞}]//N

Đầu ra

0.05692157040109628

Bước 7

(Abs[exact - res]/Abs[exact])

Đầu ra

0.13417484131607102

Bước 8

resEO = NIntegrate[Sin[ω x]f[x], {x, 0, ∞}, Method -> "ExtrapolatingOscillatory"]

Đầu ra

NIntegrate[Sin[20*x]/(1 + x)^2, {x, 0, Infinity}, Method -> ExtrapolatingOscillatory]

Bước 9

(Abs[exact - resEO]/Abs[exact])

Đầu ra

17.568032521828325*Abs[0.05692157040109628 - NIntegrate[Sin[20*x]/(1 + x)^2, {x, 0, Infinity}, Method -> ExtrapolatingOscillatory]]

N integrate integration strategies extrapolating oscillatory strategy

Các hàm được dùng