Solution (source code)

= Solution

Use index contraction:
$$
(MN)_{ij}=(\delta_{ik}+2\epsilon_{ik\ell}a_\ell)a_kb_j=a_ib_j+2\epsilon_{ik\ell}a_ka_\ell b_j=N_{ij}.
$$
The last term vanishes because $a_ka_\ell$ is symmetric in the contracted indices while the <Levi-Civita symbol> is antisymmetric. Hence $\boxed{MN=N}$.