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