Без заголовка

Грамотного там увы не так много, их операции над матрицами "монолитны" и не дают возможность на высоком уровне получить промежуточные варианты вычислений (что они и продемонстрировали, поспешив "первыми задать" неудобный для них вопрос тест). А "быстрый цикл" порождает редукцию алгоритма к нагромождению ифов, форов и присваиваний. Гораздо проще и эффективнее сделать вставку кода на C(++).