Алгоритм Хука - Дживса

Нужен алгоритм для решения задачи по оптимизационной математике по алгоритму Хука - Дживса. Дана функция a*x1^2+2*x1*x2+b*x2^2-2*x1-3*x2. В алгоритме должно быть реализовано: изменение начальной точки и точности. Должно выводится: номер итерации,значение функции и экстремум(в виде таблицы), так же траектория