pub fn kron(a: &Array2<Complex64>, b: &Array2<Complex64>) -> Array2<Complex64>
Calculate the kronecker product of two matrices