【请指点】数值计算重积分总是报错(两个明确积分限、一个变限)

数学符号变成了代码：
[CapitalTheta] = m^2 X Cos[ [Sigma]]^2 Cot[ [Theta]]^2;
R[ [Epsilon]_, [Sigma]_, F_] := ( [Epsilon]^2 - 1) r^4 +
2 M r^3 + (a^2 ( [Epsilon]^2 - 1) - X Sin[ [Sigma]]^2 - F) r^2 +
2 M (F + (Sqrt[X] Sin[ [Sigma]] - a [Epsilon])^2) r - a^2 F;

[Sigma]m = 0; [Sigma]M = 2 [Pi];
[Epsilon]m = 1; [Epsilon]M = 100;
Fm = 0.01;
Fc[ [Epsilon]_, [Sigma]_] := (-377.54 + 585.47 [Epsilon]^2 +
2.61 Sin[ [Sigma]]^2 - 2.61 [Epsilon]^2 Sin[ [Sigma]]^2)/(
13.68 + 41.78 Sin[ [Sigma]]^2);

n = 100;
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}]

请大佬指点报错部分如何修正处理。
将报错点邻域挖去后会出现新的报错点。