代码如下：
h = 6.62607015*10^-34 ;(*普朗克常数*)
c = 299792458;(*光速*)
k = 1.380649*10^(-23);(*玻尔兹曼常数*)
e = 1.6*10^-19;
b[x_] := (1 - x) (11.4 - 19.4 x);
EG[x_] := x 0.64 e + (1 - x) 3.43 e - b[x] x (1 - x) e;
Subscript[EG, 1] = 3.43 e;
Subscript[EG, 2] = EG[0.5];
\[CapitalDelta]EG = (Subscript[EG, 1] - Subscript[EG, 2]);
\[Phi][\[Epsilon]a_, \[Epsilon]b_, Ts_, \[Mu]_] := (2 \[Pi])/(h^3 c^2)
Integrate[\[Epsilon]^2/(
E^((\[Epsilon] - \[Mu])/(k Ts)) -
1), {\[Epsilon], \[Epsilon]b, \[Epsilon]a}] e;
\[Phi]2 = \[Phi][\[CapitalDelta]EG, Subscript[EG, 2], 300, \[Mu]1]
或者可以利用这个求个数值方程也可以不必要表达式

要求的方程中要对mu1求导。求教！谢谢啦

可能是数值计算误差，可以尝试直接把μ1具体的值代入，然后用数值法求导数，或者先保留所有的量进行符号计算，最后再代入具体数值

和LZ另一帖差不多是同一个问题：
https://tieba.baidu.com/p/7881092151?pid=144452920356&cid=#144452920356
用Rationalize或者SetPrecision强制变准确数就行了