4
Método para la integración numérica de la integral oscilatoria difícil.
Necesito evaluar numéricamente la integral a continuación: ∫∞0sinc′(xr)rE(r)−−−−√dr∫0∞sinc′(xr)rE(r)dr\int_0^\infty \mathrm{sinc}'(xr) r \sqrt{E(r)} dr donde E(r)=r4(λκ2+r2−−−−−−√)−ν−5/2K−ν−5/2(λκ2+r2−−−−−−√)E(r)=r4(λκ2+r2)−ν−5/2K−ν−5/2(λκ2+r2)E(r) = r^4 (\lambda\sqrt{\kappa^2+r^2})^{-\nu-5/2} K_{-\nu-5/2}(\lambda\sqrt{\kappa^2+r^2}),x∈R+x∈R+x \in \mathbb{R}_+yλ,κ,ν>0λ,κ,ν>0\lambda, \kappa, \nu >0. AquíKKKes la función Bessel modificada del segundo tipo. En mi caso particular tengoλ=0.00313λ=0.00313\lambda = 0.00313,κ=0.00825κ=0.00825\kappa = 0.00825yν= 0,33ν=0,33\nu = 0.33. Estoy usando MATLAB, y he probado las …