※ 引述《obelisk0114 (追風箏的孩子)》之銘言：
: 我想要用數值積分求下面這積分的 FWHM(半高全寬)
: http://i.imgur.com/LgocICB.png
: 我想求出最大值發生處以及其數值,然後用 FindRoot 解
: 但是第一步用 MaxValue 就跑出什麼 max precision 的問題
: 想問有什麼簡單的數值方法大致求出 FWHM
: 附上積分程式碼
: Integrate[(
: Abs[Integrate[
: E^(I*2*Pi/0.7*z*0.5^2*p^2)*BesselJ[0, 2*Pi/0.7*r*0.5*p]*p, {p,
: 0, 1}]])^2*r, {r, 0, 3552/15}]
後來用 Simpson's 3/8 rule (for n intervals)
http://en.wikipedia.org/wiki/Simpson%27s_rule
在 NIntegrate 出現精確度的問題
他說甚麼 bisection 在 9 次遞迴無法精確收斂之類的
f[z_?NumericQ, r_?NumericQ] :=
[email protected]/*