急求助：新手求助关于二重积分的MATLAB问题，谢谢

jingsuxuyi

clear;
clc;
lam=0.6551/2;
f0=inline('1/sqrt(120*pi).*exp(-(e-30).^2/120)','e');  %发射机OFF
f1=inline('1/sqrt(4000*pi/27).*exp(-27*(e-100/3).^2/4000)','e');%发射机ON
ru=inline('0.5*(1+f0./f1)','e');
fun=inline('log2(ru.*r./lam).*f1.*exp(-(4-pi)/4.*r)*(4-pi)/4','r','e');
low=inline('lam./ru','e');
up=inline('1000','e');
numint=quad2dggen(fun,low,up,0,80,1e-6)