Question: How to compute only the diagonal of a matrix product in Octave?