各位高手大家好，遇到一个棘手的积分问题想请问大家!
如果有没表达好的部分请见谅。
是这样的，我想要将一包含 x v 两变量之函数 I 对 u 作图
其中 x的范围为 0~1 v的范围为0~Vd 以比较不同 Vd 时 I(u)的分布差异
而 J0为 0阶 Bessel Function 完整式子如下
I=∫(Vd~0) |∫(0~1) exp(1/2*i*u*x^2) J0(x*v) xdx |^2 * vdv
将绝对值内的积分式对x积分并取平方后，再接着乘上v后对v积分
想请问这个问题有办法用MATLAB解决吗?
因为绝对内的函数并无反导函数，仅能求出数值解
我有先试着不考虑顺序直接将上式做双重积分最后再取绝对值的平方
想当然尔结果会有问题，当Vd>4时 (J(3.81)=0) 会出错
以下为我目前的程式码提供给大家参考
还请各位高手抽空帮我看看 谢谢大家了。