对5楼有点好奇,为什么不写成:
n = 100;
Epsilon_m = 0.1; Epsilon_M = 1.0; Sigma_m = 0.1; Sigma_M = 1.0;
para = {Beta -> 1, Theta -> Pi/4, R -> Function[{Epsilon, Sigma, F}, Epsilon + Sigma + F]};
core[Epsilon_, Sigma_, F_] := (Epsilon Exp[-Beta Epsilon])/(2 Sqrt[R[Epsilon, Sigma, F]] Cos[Theta]) /. para;
core1[Epsilon_, Sigma_] := NIntegrate[core[Epsilon, Sigma, F], {F, Fm, Fc[Epsilon, Sigma]}] /. para;
core1data = Flatten[Table[{Epsilon, Sigma, core1[Epsilon, Sigma]},
{Epsilon, Epsilon_m, Epsilon_M, (Epsilon_M - Epsilon_m)/n},
{Sigma, Sigma_m, Sigma_M, (Sigma_M - Sigma_m)/n}], 1];
corefunc = Interpolation[core1data];
NIntegrate[corefunc[Epsilon, Sigma], {Epsilon, Epsilon_m, Epsilon_M}, {Sigma, Sigma_m, Sigma_M}]