Qurak

N integrate integration strategies double exponential oscillatory integration

Tutorial 6 steps, run in order in one session. Every step was executed against the engine and the output below is what it produced.

Step 1
NIntegrate[Sin[2 * x] * (1 / x^2), {x, 1, ∞}, Method -> {"DoubleExponentialOscillatory", "SymbolicProcessing" -> 0}]
Output
NIntegrate[Sin[2*x]/x^2, {x, 1, Infinity}, Method -> {DoubleExponentialOscillatory, SymbolicProcessing -> 0}]
Step 2
exact = Limit[Integrate[Exp[-c x]Log[x]Sin[x], {x, 0, ∞}, Assumptions -> c > 0], c -> 0]
Output
Integrate[Log[x]*Sin[x], {x, 0, Infinity}]
Step 3
NIntegrate[Log[x]Sin[x], {x, 0, ∞}] - exact
Output
-8.201332694086795*^18 - Integrate[Log[x]*Sin[x], {x, 0, Infinity}]
Step 4
exact = Integrate[Sin[20 x]Cos[18x] (1/Sqrt[x + 1]), {x, 0, ∞}]
Output
Integrate[(Cos[18*x]*Sin[20*x])/Sqrt[1 + x], {x, 0, Infinity}]
Step 5
NIntegrate[Sin[20 x]Cos[18x] (1/Sqrt[x + 1]), {x, 0, ∞}, PrecisionGoal -> 10] - exact
Output
6.345921784313157*^12 - Integrate[(Cos[18*x]*Sin[20*x])/Sqrt[1 + x], {x, 0, Infinity}]
Step 6
Plot[{Sin[20 x]Cos[18x] (1/Sqrt[x + 1]), Sin[20x]}, {x, 0, 3}]
Output
-Graphics-

Functions used

Related recipes

All recipes · Function reference · Use this from an MCP client