※ 引述《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]/*