[ykuai]

函数f=y(E)

c=1/(exp((E-Ef)/(k*Te(p)))+1);
d=(1/(16*pi^2*h2^2))*Fps*1/sqrt((h2/mp2)*(h2*w+E-Hwg)-(h2/ms2)*(E+Hwf));
Ree=0.01*((pi*k*Te(p)).^2+(E-Ef).^2)/(1+exp((E-Ef)/(k*Te(p))));
g=Ree.^2/((E-h2*w(j)).^2+Ree.^2);
y=@(E)c*d*g;

式中除了E是变量外，其余的都是具体的参数，怎么计算函数y在-infinity 到 +infinity上的积分呢？

非常感谢指点。