Будем считать, что матрица записана в памяти машины по строкам, так что адрес аы равен а + (k — 1) п + /. Как всегда, программирование начинаем с составления программы внутреннего цикла. Переадресацию переменных команд будем производить после их выполнения. Рабочая часть цикла, очевидно, состоит из трех операций пересылок чисел.
Какой вид должна иметь заготовка для сравнения в ячейке г2? Параметр /" изменяется до значения i — 1, переадресация команды
происходит после ее выполнения, т. е. мы должны прекращать выполнение цикла при / = i.
Чтобы закончить программирование внутреннего цикла, нужно в его подготовительной части осуществить формирование всех переменных команд для значения.
Во внешнем цикле меняется содержимое ячеек г2, г3, г4, г5. Добавим команды переадресации (поместим их перед командами формирования внутреннего цикла).
Сравнение на окончание цикла по i будем проводить по содержимому ячейки r5.
Выпишем готовую программу, дополнив все написанные команды формированием для i = 1 (переадресация команд происходит до их выполнения, начинать цикл нужно со значения i == 2).
особых затруднений не вызвала, программа была составлена обычным образом. Возникает вопрос, нельзя ли сократить написанную программу, применив восстановление команд? Оказывается, нет, — в этом случае зависимость верхней границы значения / от параметра i препятствует сокращению: Полное изменение переменных команд в цикле по / зависит от значения г, значит, должны зависеть от ь и «константы» восстановления. Их, в свою очередь, нужно переадресовывать во внешнем цикле. Число команд в программе применением восстановления сократить не удается.
Упражнения. 1. Переделать программу 4, используя восстановление по /.
2. Составить программу транспонирования квадратной матрицы так, чтобы во внутреннем цикле элементы одной из диагоналей, параллельных главной диагонали матрицы (т. е. элементы a [i, i + k] при i = 1, n — к и фиксированном k) переставлялись с симметричными им элементами, а во внешнем цикле осуществлялся перебор всех таких диагоналей (для & = 1,      п— 1).
3. Составить программу вычисления разделенных разностей по заданной таблице значений аргумента и функции. Аппарат разделенных разностей очень удобен для интерполяции функций одной переменной при неравноотстоящих значениях аргумента. Напомним определение разделенных разностей, одновременно дающее метод их вычисления.
где т может принимать значения т = 1, 2, п, a i для данного т — значения i — О, 1, 2,      п — га.