Метод вращеющихся координат, метод Розенброка Опции

[setare]

Здравствуйте!!!! Очень нужна помощь в написании этого метода!!!! К сожалению, возникают очень много вопросов. Исходные данные: х0 (вектор), задается направление (S1,.....Sn)Б которое совпадает с направлением оси координат, задается точность.
Нужно написать этот метод и поитерационна вывести каждое полученное промежуточное значение до прихода в минимум.

Пожалуйста, очень нужно!!!! Завтра уже сдавать! Я с этим уже сижу неделю!

[setare]

есть такой алгоритм, по которому надо написать:


1) задаваемые значение (написаны ранее)
2)поиск экстремума вдоль каждого из направленй x(1)=x(k)+l(k)s(k), l1(k)=argminf(x(k)+l1s(k))....
Вопрос как это можно написать в коде? И как задавть направление?
3)получаем (хn(k)-x(k)-направление, новая точка x(k+1)=xn(k)
4) Проверяем разность модуль ( х(к+1)-х(к))<=E
Если выполняется то в качестве приближ точки экстр берем х"=х(к+1) иначе переход к шагу 5
5)проводим ортогонализацию
Вычисляем вспомогательный вектор
(А1(к)....Ан(к))=(S1(k+1).....sn(k+1))
Если l=0 то меняем направление на направления не предыдушеим шаге.k=k+1 и переход к шагу2.

Как проводить ортогонализацию? Каким образом можно все это запрограммировать?