在运行matlab程序时，出现如下所示报错，求助各位大佬该如何解决？感谢！

警告：已达到函数计算的最大数目(10000)。这个结果将使全局误差测试失败。

In funfun\private\integral2Calc>integral2t at 129

  In funfun\private\integral2Calc at 9

  In integral2 at 106

  In fourierbmn at 13

  In parallel_function>make_general_channel/channel_general at 914

  In remoteParallelFunction at 38

%主程序文件
clear all;
T=277;
syms x y;
nG=20;
nG2=2*nG+1;
amn=zeros(nG2,nG2);
bmn=zeros(nG2,nG2);
parfor m=1:nG2
for n=1:nG2
amn(m,n)=fourieramn(m-nG-1,n-nG-1,T);
end
end
parfor m=1:nG2
for n=1:nG2
bmn(m,n)=fourierbmn(m-nG-1,n-nG-1,T);
end
end


%自定义函数fourieramn.m 
function amn=fourieramn(p,q,T)
T=277;
a=40;
b=40;
c=sqrt(a^2-b^2);
c1=[-c/sqrt(2),-c/sqrt(2)];c2=[c/sqrt(2),c/sqrt(2)];
R=40;
d=T/2;
g=@(x,y)cos(p*2*pi/T*x+q*2*pi/T*y).*(8.13^(-1)*heaviside(2*a-sqrt((x-c1(1)).^2+(y-c1(2)).^2)-sqrt((x-c2(1)).^2+(y-c2(2)).^2))+8.30^(-1)*heaviside(R^2-(x-d).^2-(y+d).^2)+12.46^(-1)*heaviside(heaviside(sqrt((x-c1(1)).^2+(y-c1(2)).^2)+sqrt((x-c2(1)).^2+(y-c2(2)).^2)-2*a)+heaviside((x-d).^2+(y+d).^2-R^2)-1.1));
amn=integral2(g,-T/2,T/2,-T/2,T/2,'AbsTol',1e-2)*(1/T^2);


%自定义函数fourierbmn.m 
function bmn=fourierbmn(p,q,T)
T=277;
a=40;
b=40;
c=sqrt(a^2-b^2);
c1=[-c/sqrt(2),-c/sqrt(2)];c2=[c/sqrt(2),c/sqrt(2)];
R=40;
d=T/2;
g=@(x,y)(8.13^(-1)*heaviside(2*a-sqrt((x-c1(1)).^2+(y-c1(2)).^2)-sqrt((x-c2(1)).^2+(y-c2(2)).^2))+8.30^(-1)*heaviside(R^2-(x-d).^2-(y+d).^2)+12.46^(-1)*heaviside(heaviside(sqrt((x-c1(1)).^2+(y-c1(2)).^2)+sqrt((x-c2(1)).^2+(y-c2(2)).^2)-2*a)+heaviside((x-d).^2+(y+d).^2-R^2)-1.1)).*sin(p*2*pi/T*x+ q*2*pi/T*y);
bmn=integral2(g,-T/2,T/2,-T/2,T/2,'AbsTol',1e-2)*(1/T^2);

您好，我是有问必答小助手，您的问题已经有小伙伴解答了，您看下是否解决，可以追评进行沟通哦~

如果有您比较满意的答案 / 帮您提供解决思路的答案，可以点击【采纳】按钮，给回答的小伙伴一些鼓励哦～～

ps:问答VIP仅需29元，即可享受5次/月 有问必答服务，了解详情>>>https://vip.csdn.net/askvip?utm_source=1146287632

非常感谢您使用有问必答服务，为了后续更快速的帮您解决问题，现诚邀您参与有问必答体验反馈。您的建议将会运用到我们的产品优化中，希望能得到您的支持与协助！

速戳参与调研>>>https://t.csdnimg.cn/Kf0y

你根据报错信息行数查找相应代码