用这段程序进行数值积分一直出现下面这样的报错。该怎么解决啊
这是程序
f1=@(x,y)((y-DY1)/((x-DX1)^2+(y-DY1)^2)^0.5)*0.9*(3*x+2*y)
f1 =
包含以下值的 function_handle:
@(x,y)((y-DY1)/((x-DX1)^2+(y-DY1)^2)^0.5)*0.9*(3*x+2*y)
>> DY1==2,DX1==3,BX==1,L==3
ans =
DY1 == 2
ans =
DX1 == 3
ans =
BX == 1
ans =
L == 3
>>FX1=dblquad(f1,-BX/2,BX/2,-L/2,L/2,1.0e-3)
错误使用 dblquad>innerintegral (第 76 行)
输入必须为单精度或双精度浮点值。
出错 quad (第 67 行)
y = f(x, varargin{:});
出错 dblquad (第 55 行)
Q = quadf(@innerintegral, ymin, ymax, tol,trace, intfcn, ...