Need help with the following code. Trying to show
$\epsilon^{iklm} \epsilon_{prlm}
= - 2( \delta_p^i \delta^k_r - \delta^i_r \delta^k_p ) $
Tried the following code:
yields $2( \delta_p^i \delta^k_r - \delta^i_r \delta^k_p )$ without negative sign!
{i,k ,l,m,p,r}::Indices.
{i,k,l,m,p,r}::Integer(1..4).
\delta{#}::KroneckerDelta.
\epsilon^{i k l m}::EpsilonTensor(delta = \delta);
eq3 := \epsilon^{i k l m} \epsilon_{p r l m};
epsilon_to_delta(_);
expand_delta(_);
enter code here