[問題] dblquad 積分問題 已回收
大家好想利用dblquad來算積分
程式碼為下
clear all;
t_r1 = 1;
t_r2 = 1.002;
B = besselj(0,2*pi*10*(t_r2-t_r1))/besselj(0,0);
F=@(x,y)(x*y/(1-B^2)*exp(-(x^2+y^2)/(2*(1-B^2)))*besseli(0,(x*y*Beta/(1-B^2))))
Q_xy=dblquad(F,1,2,1,2)
結果會有錯誤
但是假如將 F改成
F=@(x,y)(x*y/(1-B^2))
則可以跑
不知道錯誤在哪邊
謝謝大家的幫忙
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.13.48