怎样matlab运算下列二重积分

上下限分别是3191.27，3165.769，3191.27，3151.27

你好同学，计算代码如下：

x2 = 3191.27;
x1 = 3165.769;
y2 = 3191.27;
y1 = 3151.27;
f = @(x,y) (155.44 + ((x-3165.039)./3165.039-0.002179)*320.56+...
    ((y-3152.0087)./3152.0087-0.002179)*320.56).*sqrt((x-3165.039).^2+...
    (y-3152.0087).^2);
S = integral2(f,x1,x2,y1,y2)

结果为

S =

     4.074974512741888e+06