各位前辈，小白在此请教：
j的取值有1,2,3,4.想要实现KroneckerDelta[j,1],和KroneckerDelta[j,3]，这一步怎么处理呢。
\[Xi]j中的j在取1,2,3,4时的函数如代码所示。
\[Xi]1 = ((9*
#1)/(4 \[Pi]) Sin[(\[Pi]*#
2)/(2*
#1)] - #
1/(12 \[Pi])
Sin[(3*\[Pi]*
#2)/(2*#
1)] ) &;
\[Xi]2 = (-((9*#1)/(4 \[Pi])) Cos[(\[Pi]*
#2)/(2*#
1)] - (#1/(12 \[Pi]))
Cos[(3*\[Pi]*
#2)/(2*#
1)] ) &;
\[Xi]3 = (
#1^3/\[Pi]^3 Sin[(\[Pi]*#
2)/(2*
#1)] - #
1^3/(3 \[Pi]^3)
Sin[(3*\[Pi]*
#2)/(2*#
1)] ) &;
\[Xi]4 = (-(
#1^3/\[Pi]^3) Cos[(\[Pi]*#
2)/(2*
#1)] - #
1^3/(3 \[Pi]^3)
Cos[(3*\[Pi]*
#2)/(2*#
1)] ) &;
func\[Xi] = {\[Xi]1[a, 0], \[Xi]2[a, 0], \[Xi]3[a, 0], \[Xi]4[a, 0]}
ffunction1 =
kx0/D0*func\[Xi] -
KroneckerDelta[j, 1]*(2 - \[Upsilon]) ((n \[Pi])/b)^2 +
KroneckerDelta[j, 3]

Table[ffunction1, {j, 1, 4}]

需要funXi 和 j 的列表一一对应？
\[Xi]1 = ((9*
#1)/(4 \[Pi]) Sin[(\[Pi]*#
2)/(2*
#1)] - #
1/(12 \[Pi]) \
Sin[(3*\[Pi]*
#2)/(2*#
1)]) &;
\[Xi]2 = (-((9*#1)/(4 \[Pi])) Cos[(\[Pi]*
#2)/(2*#
1)] - (#1/(12 \
\[Pi])) Cos[(3*\[Pi]*
#2)/(2*#
1)]) &;
\[Xi]3 = (
#1^3/\[Pi]^3 Sin[(\[Pi]*#
2)/(2*
#1)] - #
1^3/(3 \[Pi]^3) \
Sin[(3*\[Pi]*
#2)/(2*#
1)]) &;
\[Xi]4 = (-(
#1^3/\[Pi]^3) Cos[(\[Pi]*#
2)/(2*
#1)] - #
1^3/(3 \[Pi]^3) \
Cos[(3*\[Pi]*
#2)/(2*#
1)]) &;
func\[Xi] = {\[Xi]1[a, 0], \[Xi]2[a, 0], \[Xi]3[a, 0], \[Xi]4[a, 0]}
j = Range[4]
MapThread[
kx0/D0*#1 - KroneckerDelta[#2, 1]*(2 - \[Upsilon]) ((n \[Pi])/b)^2 +
KroneckerDelta[#2, 3] &,
{func\[Xi], j}
]