Question

关于半导体费米积分的积分,求助matlab大神 本人matbal很弱，最近做一个关于费米积分的数值计算

整个方程是在一个费米积分的基础上，对费米积分再积分
我写的费米积分程序如下
function f=fermi(u,E)
syms x;
n=(u-E)/0.026;
y=x.^(-0.5)./(1+exp(x-n));
f=quad(@(x)eval(vectorize(y)),0,100);
end
能跑通。
但我的目的是以fermi(u,E)中的E为积分对象，再积分，就出问题了
比如会出
Undefined function or method 'isfinite' for input arguments of type 'sym'
这种错误。
求解我想把E在(1,2)上积分，我应该怎么写外部的函数呢，还是说我费米积分的本身表达就有问题呢？

Answer 1

可以这样直接写。

 u=1;E=0.8;n=(u-E)/0.026;

 F0=quadgk(@(x) x.^(-0.5)./(1+exp(x-n)),0,inf)

运行结果

如还有问题，可以进一步交流。