Wolfram mathematica 我怎样才能找出PN结中带有Kronecker符号的四个元件？

如何找出PN结中带有Kronecker符号的四个组件？

a=1；c=1.5

  kp1 := 1/a Sqrt[(\[Epsilon])^2/11.25 - ( 2 \[Epsilon] c)/
11.25 + (c)^2/11.25 + \[Epsilon]/11.25 - c/11.25 - ( ky)^2];`/*second region,first layer*/`

kp2 := 1/a Sqrt[(\[Epsilon])^2/11.25 - ( 2 \[Epsilon] c)/
11.25 + (c)^2/11.25 - \[Epsilon]/11.25 + c/11.25 - ( ky)^2];`/*second region,second layer*/`


k1 := 1/a Sqrt[(\[Epsilon])^2/11.25 + \[Epsilon]/11.25 - ( ky)^2];`/*First region,first layer*/`


k2 := 1/a Sqrt[(\[Epsilon])^2/11.25 - \[Epsilon]/11.25 - ( ky)^2];`/*First region,second layer*/`


p = {{1, 1, 0, 0}, {(  Sqrt[11.25] k1)/\[Epsilon], -(( Sqrt[11.25] k1)/\[Epsilon]), -((
I Sqrt[11.25] ky )/\[Epsilon]), -((I Sqrt[11.25] ky )/\[Epsilon])}, {0, 0, 1, 
1}, {-((I Sqrt[11.25] ky )/\[Epsilon]), -((I Sqrt[11.25] ky )/\[Epsilon]), (
Sqrt[11.25] k2)/\[Epsilon], -((Sqrt[11.25] k2)/\[Epsilon])}};`intermediary matrix`


\[Zeta][x_] = DiagonalMatrix[{Exp[I k1 x], Exp[-I k1 x], Exp[I k2 x], 
Exp[-I k2 x]}]; `Diagonal matrix`

p1 = {{1, 1, 0, 0}, {(
Sqrt[11.25] kp1)/(\[Epsilon] - 1.5), -((
Sqrt[11.25] kp1)/(\[Epsilon] - 1.5 )), -((
I Sqrt[11.25] ky )/(\[Epsilon] - 1.5 ) ), -((
I Sqrt[11.25] ky )/(\[Epsilon] - 1.5 ) )}, {0, 0, 1, 
1}, {-((I Sqrt[11.25] ky )/(\[Epsilon] - 1.5 ) ), -((
I Sqrt[11.25] ky )/(\[Epsilon] - 1.5 ) ), (
Sqrt[11.25] kp2)/(\[Epsilon] - 1.5 ) , -((
Sqrt[11.25] kp2)/(\[Epsilon] - 1.5 ))}};`intermediary matrix`


\[Zeta]1[x_] =  
DiagonalMatrix[{Exp[I kp1 x], Exp[-I kp1 x], Exp[I kp2 x], 
Exp[-I kp2 x]}]; `Diagonal matrix`


  M = Inverse[\[Zeta][0]].Inverse[p].p1 .\[Zeta]1[0];  `/*Transfer matrix*/`

{(t +)^l, 0, (t -)^l, 0}.M = {KroneckerDelta[1, l], (r +)^l, 
KroneckerDelta[-1, l], (r -)^l};`   The coefficient of the transmission and the reflection vectors`

**任何人谁能完成这个程序，找出这八个组成部分取决于克罗内克德尔塔** 其目的是绘制E（ky，E）的密度图