教學共 6 個步驟，在同一個工作階段中依序執行。每一步都在引擎上實際執行過，下方的輸出就是它產生的結果。

步驟 1

NIntegrate[Sin[2 * x] * (1 / x^2), {x, 1, ∞}, Method -> {"DoubleExponentialOscillatory", "SymbolicProcessing" -> 0}]

輸出

NIntegrate[Sin[2*x]/x^2, {x, 1, Infinity}, Method -> {DoubleExponentialOscillatory, SymbolicProcessing -> 0}]

步驟 2

exact = Limit[Integrate[Exp[-c x]Log[x]Sin[x], {x, 0, ∞}, Assumptions -> c > 0], c -> 0]

輸出

Integrate[Log[x]*Sin[x], {x, 0, Infinity}]

步驟 3

NIntegrate[Log[x]Sin[x], {x, 0, ∞}] - exact

輸出

-8.201332694086795*^18 - Integrate[Log[x]*Sin[x], {x, 0, Infinity}]

步驟 4

exact = Integrate[Sin[20 x]Cos[18x] (1/Sqrt[x + 1]), {x, 0, ∞}]

輸出

Integrate[(Cos[18*x]*Sin[20*x])/Sqrt[1 + x], {x, 0, Infinity}]

步驟 5

NIntegrate[Sin[20 x]Cos[18x] (1/Sqrt[x + 1]), {x, 0, ∞}, PrecisionGoal -> 10] - exact

輸出

6.345921784313157*^12 - Integrate[(Cos[18*x]*Sin[20*x])/Sqrt[1 + x], {x, 0, Infinity}]

步驟 6

Plot[{Sin[20 x]Cos[18x] (1/Sqrt[x + 1]), Sin[20x]}, {x, 0, 3}]

輸出

-Graphics-

N integrate integration strategies double exponential oscillatory integration

使用到的函式