我正在尝试评估此积分，但底部的减号会导致一些错误。
Kb=1.381e-23; %J/k
h=6.626e-34; %J*s
Ts=5760; %k
Tc=300; %k
q=1.60e-19; %J
c=3.0e8; %m/s
theta_s=atan(7e8/1.5e11) %rad
theta_e=pi/2;
theta_x=4.66e-3
syms E;
Eg=1.3*1.602e-19;
Il = int(E^2/((exp(E/Kb/Ts)-1)),E,Eg,Inf)
eval(Il)
当我尝试运行代码时出现以下错误。
Error using evalin
Undefined function 'polylog' for input arguments of type 'sym'.
Error in sym/eval (line 11)
s = evalin('caller',vectorize(map2mat(char(x))));
Error in homework3_stackoverflow_help (line 15)
eval(Il)
如果我底部有exp(E/Kb/Ts)-1，则会给我错误。但是，当我删除一个。 exp(E/Kb/Ts)工作得很好。积分可以在http://pveducation.org/pvcdrom/solar-cell-operation/detailed-balance