мда, один цикл заменили на 3, и после этого говорите о оптимизации по времени.
причем в первом варинте компилятор еще смог бы сам что-то подправить,
а в вашем будет тупо крутить 3 цикла.  причем с жуткими тормозами по памяти.

чем плоха книга Горелик - теоретическим походом. Читать надо все подряд и не только по фортрану.

Законные замечания, но.. не верьте своим глазам, не стоит всё буквально понимать.

Во-первых, выбор промежуточных переменных передан компилятору. Во-вторых, кто вам сказал, что эти три цикла будут выполнятся как три? В-третьих, это было первое промежуточное приближение. При такой записи лучше видно, что делать дальше. т.к. задачу приходилось поднимать снизу вверх с "низкоуровнего" С++. При записи же потом в одну строчку на полную катушку будет работать векторизация, а это ~x3.х при больших значениях.


  А лично я её и не открываю. Она считает, что в Fortran есть указатели. Это не компетентно, формальный перевод, расходящаяся с практикой теоретическая обрисовка.