Error with ’quad‘！！！Doubleintegral?

Question

wang 2013-1-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/60735-error-with-quad-doubleintegral

L = 3.128438e+003；
B= 16685；
X = 2.085625e+003；
Y = 0；
Z = 0；
Quad_1 = @(Zeta,Xi) exp(-sqrt(Xi.^2 + Zeta.^2 + Z^2))./(-sqrt(Xi.^2 + Zeta.^2 + Z^2));
T_1 = quad(@(Zeta) Quad_1,Y-B,Y+B);
T_2 = @(Xi) exp(-Xi)*T_1;
T = quad(@(Xi) T_2(Xi), X-L, X+L);
disp(T);
 Error using ==> quad at 80
   The integrand function must return an output vector of the same length as the input vector
Error in ==> stationaer at 75
 T_1 = quad(@(Zeta) Quad_1,Y-B,Y+B);

Can anyone help me to correct the errors? Thanks.

Jack

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

bym 2013-2-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/60735-error-with-quad-doubleintegral#answer_73115

在 MATLAB Online 中打开

how about this

L = 3.128438e+003;
B= 16685;
X = 2.085625e+003;
Y = 0;
Z = 0;
Quad_1 = @(Zeta,Xi) exp(-sqrt(Xi.^2 + Zeta.^2 + Z^2))./(-sqrt(Xi.^2 + Zeta.^2 + Z^2));
%T_1 = quad(@(Zeta) Quad_1,Y-B,Y+B);
%T_2 = @(Xi) exp(-Xi)*T_1;
%T = quad(@(Xi) T_2(Xi), X-L, X+L);
T = dblquad(Quad_1,Y-B,Y+B,X-L,X+L);
disp(T);